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December 15, 1998 (19981215) 
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Semiconductors, MOS) ; R098 (ELECTRONIC MATERIALS — Charge 
Transfer Elements, CCD & BBD) ; R131 (INFORMATION PROCESSING 
-- Microcomputers & Microprocessers ) 



ABSTRACT 

PROBLEM TO BE SOLVED: To properly control an actuator by providing a neuro 
computer to which an output of a second sensor different from a first 
sensor, and controlling the actuator based on an output of the neuro - 
computer . 

SOLUTION: A signal of a knocking sensor 8 is input to a sample hold 
circuit 11 via an amplifier 9. The input signal 10 is a time series 

signal . The sample hold circuit 11 sample holds the signal 10 at a 
constant time cycle, and it is input as a special multivariate signal , 
12a, 12b, 12c,..., 12n in input order, to an input phase of a neuro - 

computer 13.. According to a signal 15 from a control circuit 2 the 
circuit 11 is controlled. An output 14 of the neuro - computer 13 is 
output as a voltage proportional to the strength of knocking in an engine, 
and then it is converted into a digital signal by a D/A converter of the 
circuit 2 to control 'the ignition timing'. Accordingly/ a proper control is 
assured. 
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— Microcomputers & Microprocessers ) 

ABSTRACT 

PROBLEM TO BE SOLVED: To alternate the control patterns of a plurality of 
actuators together according to target value, by providing a neurocomputer 

having a plurality of input and output ports, and alternating the 
weighting of the neurocomputer or the like according to the difference 

between an output signal of a certain output port, and a set point of 
the same. 

SOLUTION: A neurocomputor comprises the neuroelements 30a-30c on an input 
layer, and the singals from a sample hold circuit are successively input 
thereto through the input ports 12a-12c. The outputs of the neuroelements 
30a-30c are respectively input to the neuroelements 31a-31c of an 
intermediate layer after the weighting factors are added thereto. The 
outputs of the neuroelements 31a-31c to which a threshold value .theta,(sub 
1) is added, are input to the neuro element 32 of an output layer to which 
the weighting 'fa'ctors are ' respectively added/* "arid ' are "output with a 
threshold value .theta.(sub 2) of the neuroelement 32. On this occasion, 
whether the same is correctly output to the output layer or not, is 
determined, and in the case of YES, the weighting coefficient is increased, 
and the threshold value is reduced, to adjust the weighting factor and the 
threshold value corresponding to the difference from the correct answer. 
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Cbrp'o'ratio'n) , JP (Japan)'" " <->■•-•.—• 
08-298807 [JP 96298807] 
November 11, 1996 (19961111) 

[6] H04N-007/16; G06F-015/18 ; H04L-029/14; H04N-007/12 
44.6 (COMMUNICATION — Television); 44.3 (COMMUNICATION — 
Telegraphy); 45.4 (INFORMATION PROCESSING — Computer 
Applications ) 

:R012 (OPTICAL FIBERS); R131 (INFORMATION PROCESSING -- 
Microcomputers & Microprocessers) 



ABSTRACT 

PROBLEM TO BE SOLVED: To improve transmission efficiency, to shorten 
nonoperating time of a system and to improve the transmission quality of 
the incoming line of the system without performing an extensive 
modification work on the transmission line. 



SOLUTION: In the two-way CATV (community antenna television system) system 
in which a center station 1 and plural local stations 2 are connected via a 
tree-form transmission line 4, a relation between noise characteristic 
amount and code error ratio in the incoming line of the two-way CATV system 
is measured before s'tarting communication service and' previously making a 
neurocomputer 18 learn this relation, and after the system operation is 



started, making the neurocomputer 18 predict types of code errors using 
the noise characteristic amount observed in a flowing noise observation 
device 16, and an error correction or the change of control device is 
directed to the local stations based on this code errors. 
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WEATHER FORECASTING SYSTEM AND NEURO - COMPUTER CONTROL SYSTEM 

PUB. NO.: 06-347563 [JP 6347563 A] 
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INVENTOR (s) : EMURA NORIO 
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APPLICANT (s) : HITACHI LTD [000510] (A Japanese Company or Corporation), JP 
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APPL. NO. : 05-136267 [ JP 93136267] 
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JAPIO CLASS: 46.1 (INSTRUMENTATION — Measurement); 45.4 (INFORMATION 

PROCESSING — Computer Applications)' 
JAPIO KEYWORD: R131 (INFORMATION PROCESSING -- Microcomputers & 

Microprocessers ) 

ABSTRACT 

PURPOSE: To estimate following weather information by making a neuro - 
computer learn change of past weather information, and then inputting 
previous weather information in a time series manner. 

CONSTITUTION: Weather information Fin, Fout for inputting, outputting 
observing spots L are allocated to neurons Nin, Nout of an input layer 20, 
an output layer 22 of a neuro - computer . Weather actual measurement 
information data of the spots L measured at different times T are input 
to the layer 20, the following weather information of the spot L is 
estimated by using learning information, and output 22 as weather 
estimation information Hout. The Fin observed at the time is converted 
into weather actual measurement information Hin, and input to the layer 22 
thereby to update the learning information. The Fin before a predetermined 
time and present time point sent from the spots L to an information 
processor is converted, input to the layer 22, following weather state is 
estimated by using the leaning information,' output 22 "as "Hout, and 
converted into Fout. Thus, weather information of the spot L after a 
predetermined time is estimated. 
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ABSTRACT 

PURPOSE: To decide convergence after learning at a high speed by outputting 
the convergence decision result of an output neuron to a general flag 
provided in the neuron at each time of learning and also outputting the 
OR coupling result of general flags of all output neurons to a control 
part . 



CONSTITUTION: A register for a gener 
An arithmetic operation element 
calculates the difference between 
allowable error and outputs state 
the general flag 40. Signals 
neurons are put together in one and 
part 10. The OR element 42 store 
with the output signal of a 
control part 13 executes a branch 
the summarizing register 43 when all 
the convergence. 



al flag 40 is placed in each neuron 30. 
33 when deciding convergence conditions 
an output value and tutor data, and an 
information on. the arithmetic result to 
from general flags 40 from respective 
sent to an OR element 42 in the control 
s the result of OR arithmetic operation 
summarizing register 43. An instruction 
instruction based upon the contents of 
learning patterns are learnt to decide 
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JP 



ABSTRACT 

PURPOSE: To generate a high precise multiple regression model by a method 
for generating the multiple regression model of an object by specifying the 
multiple regression constant matrix of the multiple regression model formed 
between input data indicated by plural input components inputted to the 
object and output data indicated by a single output component outputted 
from the object. 

CONSTITUTION: A neurocomputer 1 constituted of a single input unit for 
receiving fixed data, the input unit for receiving the input data indicated 
by the input components and a single output unit for calculating and 
outputting • the value of the sum o£ the .products of: the data received by the 
input unit and a weight value is prepared. At first, when the measured 

input data are inputted to the neurocomputer 1, the weight values are 
learned so that the measured output data are outputted from the 

neurocomputer 1. Then, the learned weight data are specified as the 
multiple regression constant matrix and the multiple regression model is 



generated. 
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ABSTRACT 

enable a state of a layer to be detected when a plurality of 
are distributed. 



CONSTITUTION: A title item is a method for detecting a layer state when 
upper and lower different liquids form layers in a target system, a 
measurement value which is obtained from the target system in time series 
or a value based on the measurement value is input to a neuro computer , 
a signal which is output' from the neurocomputer is compared with a 
signal which indicates the layer state which is known by actual sampling 
a this time / a weight indicating a strength of combination from an input 
to an output in the neurocomputer is adjusted, an input/output model with 
the adjusted weight as a coefficient is created, the measured value or a 
value based on the measured value is "input to the input/output model, and 
the state of the layer is recognized from the output at this time . 
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PURPOSE: To 



ABSTRACT 

recognize an operation having an individual difference such 



as* talking with hands by constituting an operation learning part and an 
operation recognition part by means of neuro computers consisting of 
input layers, hidden layers, output layers and context layers, and 
inplanting a learning result to the operation recognition part through a 
weight storage part. 

CONSTITUTION: An operation measurement part 10 measures the three- 
dimensional movement of hands and outputs it as time sequential data. A 
learning data generation part 14 generates learning data from the time 
sequential data and a teach signal . The neuro computer 16-1 has a 
recurrent function and learns a weight from the learning data. The weight 
storage part 18 stores the weight of network connection at every operation. 
The weight of the weight storage part 18 is inplanted to the neuro 
computer 16-2 in the operation recognition part 20. When time 
sequential data is inputted from the operation measurement part 10, the 
neuro computer 16-2 executes the algorithm of the recurrent type neural 
network and outputs a recognition result. 
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ABSTRACT 

PURPOSE: To reduce the number 1 of "times of* learning, and to improve 
recognition rate by inputting time series data after removing its DC 
component by differentiating the time series data to the input layer of a 
neuro - computer . 

CONSTITUTION: A differentiation unit 1 differentiates the inputted time 
series data, and removes its DC component, and inputs the time series 
data after removing the DC component to the input layer constituting the 
neuro - computer 2. Then, a teacher signal is inputted to an output 
layer in accordance with it, and after repeating the learning, the 
inference of the event of comming time , etc., is executed. Thus, by 
inputting the time series data after eliminating its variation due to 
the DC component by differentiating the time series data to the neuro — 
computer , the number of times of the learning can be reduced, and the 
recognition rate can be improved. 
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ABSTRACT 

PURPOSE: To enable attainment of correct judgement .on the necessity of re- 
inspection by interpreting measured data of a specimen inspecting 
apparatus in a neuro - computer and by detecting an error in measurement 
therein . 



CONSTITUTION: When measured data showing the result of inspection on a 
specimen are outputted to a neuro - computer 2 from a specimen inspecting 
apparatus 1, the computer 2 determines whether or not an error in 
measurement is expected from the measured data, and outputs to the 
apparatus 1 the result as to whether re-inspection is needed or not. The 
apparatus 1 receiving the result of decision of the computer 2 as to the 
error in measurement in this way can display 4 this result. When an 
operator P recognizes that the content of display on the device 4 is 
different from his own judgement, on the occasion, he modifies teaching 
data by operating a keyboard 5. When the content of display on the device 4 
is in accord with judgement of the operator P and shows that the 
re-inspection is needed, he can execute the re-inspection in accordance 
with the contents. Then, the measured data of the apparatus 1 and the 
result of decision of the re-inspection by the computer 2 are -stored 3.- 
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ABSTRACT 

PURPOSE: To optically execute an intellectual information processing of a 
complicated pattern recognition, etc. which necessitate a learning function 
by constituting the coupling strength of neurones between each layer, of a 
real time hologram, and constituting an intermediate layer, an output 
layer and an error correcting signal generating part, of a space optical 
modulator or a non-linear light optical element. 

CONSTITUTION: When a hologram 28 is irradiated by an optical input signal 
27, a diffracted light 29 is outputted, and this diffracted light is 
brought to threshold processing optically by an intermediate layer 6, and 
thereafter, radiated to a hologram 34, and this diffracted light 36 is 
brought to threshold processing optically by an output layer 10, and 
thereafter, outputted as an optical output signal 11. On the other hand, 



in* an error correcting signal generating part 13, an error signal is 
generated from a difference between the optical output signal 11 and a 
teacher signal 12 applied from the outside, and this error signal is 
fed back to the hologram 28 and recorded. That is, by error signals 14, 
16, weight of coupling between neurones between the intermediate layer 6 
and the output layer 10, between an input layer 2 and the intermediate 
layer 6, and plural intermediate layers 6 is < corrected, and learning is 
executed. In such a way, an optical neuro - computer for executing an 
intellectual information knowledge of a complicated pattern recognition, 
etc. which necessitate a learning function by utilizing an advantage of a 
parallel property and a high speed property, etc. can be realized. 
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NOVELTY - Several connectors (80A-80E) are operably coupled with 
corresponding oscillators (60A-60E) . A conductive medium (70) is 
operably coupled to the connectors, simultaneously applies oscillatory 

signal to each oscillator via the connector. A sourcing apparatus 
with rhythmic external forcing input (90) generating an oscillatory 
signal , is operably coupled with the medium. 

DETAILED DESCRIPTION - An INDEPENDENT CLAIM is also included for 
communication establishing method between two oscillator having 
different frequencies . i ^ 

USE - For simulating oscillatory nature of brain neurons. 
ADVANTAGE - Neuro - computer can act as a classical fully 
connected Hopfield network even when there are only interconnections. 

DESCRIPTION OF DRAWING (S) - The figure shows the schematic diagram 
of neural network having five neural processing elements. 
Oscillators (60A-60E) 
Conductive medium (70) 
Connectors (80A-80E) 
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Model comprises two interconnected opto-electronic operational 
modelling medium distribution units, and now has additional inverse 
coupling units, electron-optical indicators. Also, it is now proposed 
that conversion of optically continuous signals into discrete 
signals is carried out in the neurite-type units, and conversion of 
the optical * discrete signals v into 4 continuous"' optical signals* is' 
carried out in the synapse-type units. The optical outputs" of the 
opto-electronic feedback units electron-optical indicators on a planar 
photoresistor, modelling divergence of the optical signals , are 
connected optically to the optical inputs of a planar operational 
photoresistor modelling convergence of input signals .* This enables 
realisation of the function of adaptation and learning in neuron 
circuits . 

USE - Model concerns opto-electronic neuro-type computing 
structures and is for use as coprocessors in neuro - computer 
systems. 
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Time sharing log input signals and weight data, which are 
inputted sequentially via analog signal buses, are subjected to 
sum/product operation. An analog neuron processor (ANP) outputs a 
signal to the analog signal buses through a non-linear circuit. The 
neutral network is controlled by reading required data "from a control- 
pattern memory and required weight data from a weight memory under the 
control of a microsequencer, to obtain a practically operable 
neuro computer . 

In this neuro computer , a number of ANPs are connected by one 
analog bus. This enables the number of wires in the neutral network to 
be reduced greatly and the scale of the circuit to be minimised. 
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...SPECIFICATION 7A illustrates of a fourth embodiment of the present 

invention. The present embodiment is directed to a neuro - computer . In 
the figure, like reference characters are used to designate corresponding 
parts to those in Figure 6... 

. ..1 and the trays 2. 5a designates a clock generator which may be 
constructed from a crystal oscillator , and 5b designates a clock 
distributor which may be constructed-^: rom a buf fer- ei-rcuit .In addition, 
101. . . 

...end of learning on the basis of a plurality of values of output errors 
to bring the neuro - computer to a stop. 102 designates the whole of 
the neuro - computer . 

Figure 7B illustrates a neuron model which is a basic element of the 
neuro-computer of the. . . 

...the nonlinear function f. 

Figure 7C is a conceptual diagram of a hierarchical neural network 
constituting a neuro - computer of a three-layer structure comprised of 
an input layer, an intermediate layer and an output layer. . . 

...layer. The outputs are drawn to the outside. In this neural network, an 
error between a teacher signal corresponding to input pattern signals 
applied to the input layer at the time of learning and an output 
signal from the output layer is obtained. The weights between the 
intermediate layer and the output layer and. . . 

...weight values determined by the law**of . back propagation learning, if -an 
imperfect pattern which is slightly different from a pattern to be 
recognized is applied to the inputs of the first layer, an output signal 
corresponding to the pattern is output from the output layer. The 



signal is very similar to a teacher signal corresponding to the 
pattern applied at the time of learning. If the difference between 
the output signal and the teacher ^signal is t ve.ry ,small,_ the, imperfect 
pattern will become recognized. 

The operation of the neural network can be... Figure 7E is a flowchart 
of the learning process in the fourth embodiment. The learning in the 
neuro - computer serves to correct the weight of earh neuron until the 
network satisfies a desired input and output relationship. The method of 
learning is to prepare a plurality of pairs of an input signal vector 
and a teacher signal vector, namely, the number of the pairs 
corresponds to the number of teacher signals in a teacher signal set. 
Then, one pair is selected from among the pairs, its input signal Ip) ) 
is entered into the network as an object of learning and the output of 
the network for that input is compared with the correct output signal , 
that is, the teacher signal Op) ) corresponding to the input signal . 
The difference is referred to as an error. The weights of the neurons 
are corrected on the basis of the error and the values of the input and 
output signals at that time . This process is repeated for each element 
of the set of teacher signals until the learning converges... 
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Claims 
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English Abstract 

A neurocomputer (50) comprises n oscillating processing elements 
(60A, 60B, 60C, 60D and 60E) that communicate through a common medium 
(70) so that there are required only n connective junctions (80A, 80B, 
80C, 80D and 80E) . A rhythmic external forcing input (90) modulates the 
oscillatory frequency of the medium (7 0) which, in turn, is imparted to 
the n oscillators (60A, 60B, 60C, 60D and 60E) . Any two oscillators 
oscillating at different frequencies may communicate provided that 
input ! s power spectrum includes the frequency equal to the difference 
between the frequencies of the two oscillators in question. Thus, 
selective communication, or dynamic connectivity^ .bet.ween different 
neurocomputer oscillators occurs due to the frequency modulation of 
the medium (70) by external forcing. 



French Abstract 



' Cet ordinateur neuromimetique (50) comprend n elements de traitement 
oscillants- (6GA,.-60B, 60C, -60D &t 60R) qui communiquent par- * - 
1 1 intermediaire d ! un support commun (70) de sorte que seulement n 
jonctions de connexion (80A, 80B, 80C, 80D et 80E) sont necessaires. 
L 1 entree d ! une contrainte, exterieure, rythmique (90) module la frequence 
oscillatoire du support (70), laquelle est a son tour appliquee aux n 
oscillateurs (60A, 60B, 60C, 60D et 60E) . Deux oscillateurs quelconques, 
oscillant a des frequences differentes, peuvent communiquer pourvu que le 
spectre de la puissance d f entree comprenne la frequence egale a la 
difference entre les frequences des deux oscillateurs en question. Ainsi, 
il se produit une communication selective, ou une connectivite dynamique, 
entre differents oscillateurs de l 1 ordinateur neuromimetique, par suite 
de la modulation de frequence du support (70) au moyen d ! une contrainte 
exterieure . 

OSCILLATARY NEUROCOMPUTERS WITH DYNAMIC CONNECTIVITY 

Main International Patent Class: G06F-015/18 
International Patent Class: G06F-015/80 
Fulltext Availability: 

Detailed Description 

Claims 

English Abstract 

A neurocomputer (50) comprises n oscillating processing elements 
(60A, 60B, 60C, 60D and 60E) that communicate through a common medium 
(70) so that there are required only n connective junctions (80A, 80B, 
80C, 80D and 80E) . A rhythmic external forcing input (90) modulates the 
oscillatory frequency of the medium (70) which, in turn, is imparted to 
the n oscillators (60A, 60B, 60C, 60D and 60E) . Any two oscillators 
oscillating at different frequencies may communicate provided that 
input's power spectrum includes the frequency equal to the difference 
between the frequencies of the two oscillators in question. Thus, 
selective communication, or dynamic connectivity, between different 
neurocomputer oscillators occurs due to the frequency modulation of 
the medium (70) by external forcing. 

Detailed Description 

OSCILLATORY NEUROCOMPUTERS WITH DYNAMIC 
CONNECTIVITY 

The present application claims priority rights based on U.S. Provisional 
Application Serial No... 

...input I 0 vector produces the desired output vector. 

Because of their ability to simulate the apparently oscillatory nature 
of brain neurons, oscillatory neuro computers are among the more 
promising types of neurocomputers . Simply stated, the elements of an 
oscillatory neurocomputer consist of oscillators rather than 
amplifiers or switches. Oscillators are mechanical, chemical or 1 5 
electronic devices that are described by an oscillatory signal 
(periodic, quasi-periodic, almost periodic function, etc.). Usually the 
output is a scalar function of the form V((ot+@) where V is a fixed wave 

forin (sinusoid, saw-tooth or square wave ) , 0) is the frequency of 
oscillation , and @ is the phase deviation (lag or lead) . 

Recurrent neural networks have feedback paths from their outputs back to 
their inputs. As... 

...connective junctions for every n processing elements employed thereby. 

In a preferred embodiment of the invention, the neurocomputer comprises 
n. oscillating processing .-elements that can communicate through a common 
medium so that there are required only n connective junctions. A 
rhythmic external forcing input modulates the oscillatory frequency of 
the medium which, in turn, is imparted to the n oscillators . Any two 
oscillators oscillating at different frequencies may communicate 
provided that the input's power spectrum includes the frequency equal to 



the difference between the frequencies of the two oscillators in 

question. Thus, selective communication, or dynamic connectivity, between 

different neuro computer oscillators occurs due to frequency 
modulation of the medium by external forcing. 

3 

BRIEF DESCRIPTION OF THE 
DRAWING, .'.on the input. 

I 0 DESCRIPTION OF THE PREFERRED EMBODIMENT 

FIG. I schematically illustrates a conventional recurrent neurocomputer 
10 comprising n (in this case, n=5) neural processing elements 20. 
Elements 20 may comprise switches, amplifiers, oscillators or any other 
suitable neurocomputer element type known in the art. In order for each 
of elements 20 to communicate with the others 1 5 of elements 20, 
neurocomputer 10 necessarily includes n 2 (in this case, n 2 =25) 
connective junctions 30 to which conductors... 

.can be observed, where the number n of elements 20 grows large, the 
implementation of such a neurocomputer becomes prohibitively difficult, 
from both cost and practicability standpoints. 

FIG. 2 schematically illustrates a neurocomputer 50 according to 
principles of the present invention. Neurocomputer 50 comprises a 
finite number n (in this case, n=5) oscillatory neural processing 
elements 60A, 6013, 60C, 60D and 60E. Elements 60A, 6013, 60C, 60D and 
60E can comprise voltage-controlled oscillators , optical oscillators , 
lasers, microelectromechanical systems, Josephson junctions, 
macromolecules, or any other suitable" oscillator known in "the art. Each 
element 60A, 6013, 60C, 60D and 60E oscillates at a particular 
frequency that may or may not be the same frequency as that of 
5 

the others of elements 60A, 6013, 60C, 60D and 60E. In its most general 
sense, the neurocomputer 50 further comprises a medium 70 connected to 
each of elements 60A, 6013, 60C, 60D and 60E. . . 

.80B, 80C, 80D and 80E, respectively. Medium 70 may comprise a unitary 
body or multiple connected bodies . 

Neurocomputer 50 further comprises a rhythmic forcing signal source 
90 able to apply a modulated oscillatory frequency to medium 70 by 
means of a connection 100. 

Specifically, the medium 70 can be a... 

.70 enables elements 60B and 60E to communicate data to each other. 

Mathematical analysis of the said neurocomputer architecture, which is 
based on the theory developed by F. C. Hoppensteadt and E. M. Izhikevich 
( Oscillatory' neurocomputers with ^dynamic connectivity, Physical 
Review Letters 82(1999)29832986) shows that the neurocomputer dynamic is 
equivalent to. . . 

.mechanisms and architectures, Neural Networks 1(1988)17-61) to show that 
a network of n - 60 oscillators can memorize and successfully retrieve 
through associative recall three patterns corresponding to the images 
11011, "I 1% 112", as we illustrate in Figure 3. Thus, the neurocomputer 
can act as a classical fully connected Hopfield network despite the fact 
that it has only n...form (4) with cij = sij for all i and j . 

(vi) Initializing the Network 

To use the proposed neurocomputer architecture to implement the 
standard Hopf ieldGrossberg paradigm, as we illustrate in FIG. 3, we need 
a way. . . 

. field activity, " the former task requires some ingenuity since we do not 
have direct access to the oscillators . 

Suppose we are given a vector e E Rn to be recognized. Let us apply the 
external . . . 



Claim 
... of 

I 0 said connectors being operably coupled with a corresponding one of 
said elements , 

3 The neurocomputer of claim 1», wherein: - - . 
said forcing apparatus comprises a rhythmic input. 

4 The neurocomputer of claim 1, wherein: 
said elements comprise oscillators . 

1 5 5. A neurocomputer comprising: 

a plurality of n oscillating processing elements; 

a plurality of no more than n connectors, each of said connectors 

being operably coupled with a corresponding one of said elements; 

a conductive medium operably coupled with said connectors; and 

a rhythmic input operably coupled with said medium. 

6 A neurocomputer comprising: 

a plurality of n processing element means; 

a plurality of connectors operably coupled with said element 
means ; 

I 8 

means for simultaneously applying an oscillatory signal to each of 
said element means via said connectors; and 

means for generating said oscillatory signal operably coupled with 
said -means- for applying. - *. * - • . . * . 

7 The neurocomputer of claim 6, wherein: 
said plurality. . . 

. . . each of said 

connectors being operably coupled with a corresponding one of said 
element means . 

8 The neurocomputer of claim 6, wherein: 
said element means comprise oscillators . 

10 9. The neurocomputer of claim 6, wherein: 

said means for applying comprises a conductive medium. 

10 The neurocomputer of claim 6, wherein: 

said means for generating comprises a rhythmic input. 

II An oscillatory neurocomputer comprising: 
15a number n of oscillating elements 

a source of a rhythmic forcing input, 

a- medium i-nter connecting the source o*f * rhythmic 1 * forcing input^to 
each 

oscillating element, 
each oscillating element having an oscillating frequency, 
the oscillating frequency f, of at least one of the oscillating 
elements differing from the oscillating frequency f2 of at least one 
other of the oscillating 
elements , 
19 

the source of a rhythmic forcing input producing an input 
of a third frequency f3, establishing communication between the at least 
one oscillating element and the at least one other oscillating 
element . 

12 An oscillatory neurocomputer according to claim I 1, wherein f3 is 
substantially the difference between f, and f 2 . 

13 An oscillatory neurocomputer according to claim 1 1, ftirther 
comprising a number nj of connections of the source of a rhythmic 
forcing input to the oscillating 

elements, • wherein - • - t «~ . . . 



nj <n. I 0 14. An oscillatory neuro computer according to claim 12, 
further comprising a number n, of connections of the source of a 
rhythmic forcing input to the oscillating 

elements, wherein 
nj < n. 

15 An oscillatory neurocomputer according to claim I 1, wherein the 
oscillating elements are electronic oscillators , the source of a 

rhythmic forcing input is a function generator and the interconnecting 
medium is an electrically conductive medium electrically connecting the 
source of a rhythmic forcing input to the oscillators - 

16 An oscillatory neurocomputer according to claim 15, wherein the 
function generator provides a forcing signal having a carrier frequency 
and information content modulating the carrier frequency, the 
oscillators responding to the impression of the forcing signal onto the 
conductive medium to produce inforination content modulation 
substantially the same as that of the conductive medium. 

17 An oscillatory neurocomputer according to claim I 1, wherein the 
number 

n of 
20 

oscillating elements is greater than two, a first subset of the 
oscillating elements communicate at a frequency f3 of rhythmic forcing 
input from the source, and ,at le v ast ojie. second ..subset ,of the, oscillating 

elements communicate at least one further frequency f4 of rhythmic 
forcing input from the source. 

18 An oscillatory neurocomputer according to claim 15, wherein 
content 

varying one 

oscillator from its oscillating frequency is communicated to and 
varies from its oscillating frequency another oscillator in 
communication with the one oscillator . 

19 A neurocomputer including : 

I 0 (a) an array of oscillators , at least a plurality of said 
oscillators having 
differing frequencies, 

(b) a common conducting medium connected to each of the plurality 
of 

oscillators , 

(c) a source connected to the conducting medium to impart oscillator 
signals 

of various, frequencies to the conducting, mediuiji, t .the. .signals, of, various 
frequencies including frequencies effective to bring two or more of the 
oscillators into communication. 

20 An oscillatory neurocomputer according to claim 19, wherein the 
oscillators include feedback circuits connected with the medium. 

21 An oscillatory neurocomputer according to claim 20, wherein the 
oscillators are phase locked loops. 

21 

. A method of enabling communication of a characteristic between a first 
processing element oscillating at a first frequency and a second 
processing element oscillating at a second frequency different from the 
first frequency, the method comprising the steps of 
operably coupling. . . 

. to a medium; 

operably coupling the second element to said medium; 

operably coupling said medium to a rhythmic input; and 

causing said rhythmic input to oscillate said medium at a third 

frequency.. , ^ ....... . . . , . 
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Title: Epilepsy in a chaos neuro - computer model 
Author: Inoue, Masayoshi; Nakamoto, Kenji 
Corporate Source: Kagoshima Univ., Kagoshima, Jpn 
Conference Title: Chaos in Biology and Medicine 

Conference Location: San Diego, CA, USA Conference Date: 
19930712-19930713 , s ...... 

Sponsor: SPIE - Int Soc for Opt Engineering, Bellingham, WA USA 
E.I. Conference No.: 19894 

Source: Proceedings of SPIE - The International Society for Optical 
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Publication Year: 1993 

CODEN: PSISDG ISSN: 0277-786X ISBN: 0-8194-1285-6. 
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Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 
Journal Announcement: 9403W2 

Abstract: A network of chaos elements has been presented as an 
information processor where each element consists of two oscillators and 
it acts' as a neuron by making use of the synchronized state of the two 
oscillators . The model is considered as a dynamical model of the brain, 
and brain dynamics is metaphorically analyzed with the use of the model. 
The time sequences of Hopfield f s energy, which are generated by the network 
when it solves a traveling salesman problem, are investigated with the use 
of a fluctuation spectrum theory. The change of the energy reflects the 
active motion of neurons, and we consider that the time sequence 
corresponds to a brain wave. If the control parameters .of the .neuron are 
chosen properly, the model can efficiently find the solution where a low 
intermittent ?brain wave 1 is observed. On the other hand, the model will 
have epileptic fits if a certain control parameter takes a small value. 20 
Ref s . 

Descriptors: ^Bioelectric phenomena; Chaos theory; Neural networks; 
Mathematical models; Computer simulation 
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Journal Announcement: 9006 

Abstract: A brief discussion of methods for interfacing neurocomputers 
with host computers is presented, followed by a discussion of 
neurocomputer performance measures. Both of these issues are central to 
the development of neurocomputers to support neural network research as 
well as the development and deployment of practical applications of 
neurocomputing . In using the author's neurocomputer performance measures, 
a methodology that starts with operational requirements assessment is often 
a useful approach. The first step is to define the class of applications 
that are to be solved by the neurocomputer . From these, a list of neural 
networks of specified architectural types and sizes that will "have to run 
(along with duty-cycle information such as the expected mean and standard 
deviation of the time durations of the runs for each network, as well as 
the expected frequency of those runs as a function of the time - of day) . 
The end result of this analysis is a table expressing the neural network 
implementation requirements that must be met by the neurocomputer . 7 
Refs. 
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• Abstract: The authors' study of thalamo-cortical systems suggests a new 
architecture for a neurocomputer that consists of oscillators having 
different frequencies and that are connected weakly via a common medium 
forced by an external input. Even though such oscillators are all 
interconnected homogeneously, the external input imposes a dynamic 
connectivity. The authors use Kuramoto's model to illustrate the idea and 
to prove that such a neurocomputer has oscillatory associative 
properties. Then, they discuss a general case. The advantage of such a 
neurocomputer is that it can be built using voltage controlled 
oscillators , optical oscillators / lasers', * microelectromechanical 
systems, Josephson junctions, macromolecules , or oscillators of other 
kinds. (16 Refs) 
Subfile: ABC 
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Abstract: A unique architecture of a winner search hardware has been 
developed using a novel neuron-like high- functionality device called the 
neuron-MOS transistor (or nu MOS' for "short) as * a v key "circuit element. 'The 
circuits developed in this paper can find the location of the maximum (or 
minimum) signal among a number of input data on a continuous-time basis, 
thus enabling real-time winner tracking as well as fully-parallel sorting 
of multiple input data. We have developed two circuit schemes. One is an 
ensemble of self-loop-selecting nu MOS ring oscillators , finding the 
winner as an oscillating node. The other is an ensemble of nu MOS 
variable-threshold inverters, receiving a common ramp-voltage for 
competitive excitation, where data sorting is conducted through consecutive 
winner search actions. Test circuits were fabricated by a 
double-polysilicon CMOS process and their operation has been experimentally 
verified. (10 Refs) 
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Abstract: A network of 'coupled chaos oscillators 1 is presented as an 
information processor where the oscillators act as a neuron. The network 
runs on a deterministic rule, but it has the ability of stochastic search. 
Parallel synchronous computation can be carried out on the digital network 
and its power is illustrated. (7 Refs) 
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ABSTRACT: This paper describes a preliminary study for VLSI implementation 
of oscillator network model LEGION for image segmentation proposed by 
D.L. Wang, et al. Calculation precision needed for the proper 
operations is estimated to be 5bits by using software simulation. Since 
this model is described by differential equations with nonlinear 
functions, VLSI implementation using conventional neural network 
circuit architecture is difficult. A coupled oscillator circuit, 
which is the most important part of LEGION, is proposed by using a new 
pulse-width ^modulation ( PWM) method" The basic operation of this' 
proposed circuit is confirmed by circuit simulation, (author abst.) 
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ABSTRACT: With the aim of realizing the intellectual function of the brain, 
a research was carried out on a vibrational amplitude learning by pulse 
information processing and moderationism learning by the neural 
network. This paper shows that in pulse newral network the network 
structure unique to the input pulse column is acquired by learning. 
This paper shows that moderationism learning is a more robust learning 
method compared to conventional direct current moderationism and 
adaptively acquires the stable oscillation in the network. 
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ABSTRACT: By noticing dynamic nonlinear characteristics of the receptor 
(sensor) and the signal processing system (neural network) , it was 
aimed to ultimately realize a new type of information processing sysem 
similar to an organism on a silicon substrate. Especially, results of 
the following two studies are reported: extraction of gaseous species 



• information by the semiconductor gas sensing system using the nonlinear 
characteristic and information processing using the multiple stability 
mode of the network (dynamic neural network) of the nonlinear 
oscillator • 
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ABSTRACT: When N oscillators are coupled by one resistor, we can see 

N-phase oscillation / because the system tends to minimize the current 
through the coupling resistor. Moreover, when the hard oscillators 
are coupled, we can see N, N-l, , 3, 2-phase oscillation and get 
much more phase states. In this study, the two types of coupled 
oscillators networks with third and fifth-power nonlinear 
characteristics are proposed. One network has two-dimensional hexagonal 
structure and the other has two-dimensional lattice structure. In the 
hexagonal circuit, adjacent three oscillators are coupled by one 
coupling resistor. On the other hand, in the lattice circuit, four 
oscillators are coupled by one coupling resistor. In this paper we 
confirm the phenomena seen in the proposed networks by circuit 
experiments and numerical calculations. In the system with third-power 
nonlinear characteristics, we can see the phase patterns based on 
3-phase oscillation in the hexagonal circuit, and based on anti-phase 

oscillation in lattice circuit. In the system with fifth-power 
nonlinear characteristics, we can see the phase patterns based on 
3-phase and anti-phase oscillation in both hexagonal and lattice 
circuits. In particular, in these networks, we can see not only the 
synchronization based on .3-phase and anti-phase oscillation .but the 
synchronization which is not based on 3-phase and anti-phase 
oscillation - As a result, these networks are expected to generate 
various synchronization patterns. In these networks, each oscillator 
is connected to only its adjacent oscillators and various patterns 
are generated according to the initial condition. Therefore, we can 
consider that we can use these networks as a kind of cellular neural 
networks, (author abst.) 
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LANGUAGE: English COUNTRY OF PUBLICATION: Japan 

DOCUMENT TYPE: Journal 
ARTICLE TYPE: Commentary 
MEDIA TYPE: Printed Publication 

DESCRIPTORS: biomodel; information processing; neural network model; 

organic semiconductor; learning model; nonlinear vibration; thin film; 
neuro computer 

BROADER DESCRIPTORS: model; treatment; semiconductor; organic conductor; 

conductor; object; oscillation ; membrane and film; computer; hardware 
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Performance Characteristics of Analog MOS Type Neuro-Processor . 
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VOL. 77, NO. 11, PAGE. 2296-2305, FIG. 13, TBL.6, REF.13 
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UNIVERSAL DECIMAL CLASSIFICATION: 681.3:007.52 
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DOCUMENT TYPE: Journal 
ARTICLE TYPE: Original paper 
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DESCRIPTORS: MOS integrated circuit; analog integrated circuit; 
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DOCUMENT TYPE: Journal 
ARTICLE TYPE: Original paper 
MEDIA TYPE: Printed Publication 

ABSTRACT: We propose a fast computation method of stochastic relaxation for 



• the continuous-valued Markov random field (MRF) whose energy function 
is represented in the quadratic form. In the case of regularization in 
visual information processing, the probability density function of a 
state transition can be transformed to a Gaussian function, therefore, 
the probablistic state transition is realized with Gaussian random 
numbers whose mean value and variance are calculated based on the 
condition of the input data and the^ neighborhood.. Early visual . 
information processing can be represented with a coupled MRF model 
which consists of continuity and discontinuity processes. Each of the 
continuity or discontinuity processes represents a visual property, 
which is like an intensity pattern, or a discontinuity of the 
continuity process. Since most of the energy function for early visual 
information processing can be represented by the quadratic form in the 
continuity process, the probability density of local computation 
variables in the continuity process is equivalent to the Gaussian 
function. If we use this characteristic, it is not necessary for the 
discrimination function computation to calculate the summation of the 
probabilities corresponding to all possible states, therefore, the 
computation load for the state transition is drastically decreased. 
Furthermore, if the continuous-valued discontinuity process is 
introduced, the MRF model can directly represent the strength of 
discontinuity. Moreover, the discrimination function of this energy 
function in the discontinuity process, which is linear, can also be 
calculated without probability summation. In this paper, a fast method 
for calculating the state transition probability for the 
continuous-valued MRF on the visual information processing is 
theoretically explained. _ (abridged, ^author abst t ..).. , , , . _ x 
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JOURNAL NUMBER: Z0909AAV 
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DOCUMENT TYPE: Conference Proceeding 
ARTICLE TYPE: Original paper 
MEDIA TYPE: Printed Publication 

ABSTRACT: In an optical neuro-element using a thin film waveguide, 

refractive index distribution which controls weighted production was 
analyzed. It is clarified that differences in refractive index of a 
thin film wave guide are depend on incident light wave features, 
which represent neurons, such as impulse or beam, and frequencies , 
and it is shown that in case of beamed incident light waves , 
variation of refractive index distribution is concentrated at central 
part of a wave guide. SN ratios for a signal light interval and 
input/output light wave intensity distributions were illustrated and 
explained . 
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Identification of the Nonlinearity of Robot Manipulators by a Neural 
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TBL.3, REF.6 
JOURNAL NUMBER: F0045BAL ISSN NO: 0387-5024 
UNIVERSAL DECIMAL CLASSIFICATION: 007.52:681.51 
LANGUAGE: Japanese COUNTRY OF PUBLICATION: Japan 

DOCUMENT TYPE: Journal 
ARTICLE TYPE: Original paper 
MEDIA TYPE: Printed Publication 

ABSTRACT: In order to ensure the high-speed and high-precision control of 
trajectory tracking for robot manipulators, it is necessary to 
construct a.pr.eciae dynamic model by which the* manipulator is<- - 
controlled. In the motion equation of the manipulator, the driving 
force is a very important term from the standpoint of control. We 
assume the nonlinearities of the manipulator and driving system are 
expressed as nonlinear functions between the driving force and the 
control input of the driving motors. This paper proposes an 
identification method for the nonlinearities of a robot manipulator by 
means of a neural network. It was experimentally confirmed that the 
proposed method is effective in the case of trajectory tracking control 
for a 2-link manipulator, (author abst.) 

DESCRIPTORS: manipulator; neural network; neurocoinputer ; system 

identification; nonlinear vibration; linkage (machine element); control 
system (computer) ; positioning 

BROADER DESCRIPTORS: robot; network; computer; hardware; identification; 
recognition; oscillation ; motion mechanism; mechanism of machine 
element; mechanism; method 
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ARTICLE TYPE: Original paper 
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ABSTRACT: A novel harmonic retrieval algorithm is proposed in this paper 

based on Hopfield's neural network. Frequencies can be retrieved with 



■* high accuracy and high resolution under low signal to noise 
ratio (SNR). Amplitudes and phases in harmonic signals can also be 
estimated roughly by an energy constrained linear projection approach 
as proposed in the algorithm. Only no less than 2q neurons are 
necessary in order to detect harmonic signals with q different 
frequencies , where q denotes the number of different frequencies 
in harmonic signals . Experimental simulations show fast convergence 
and stable solution in spite of low signal to noise ratio can be 
obtained using the proposed algorithm, (author abst.) 

DESCRIPTORS: Hopfield model; neurpcomrmter ; SN ratio ;_optimizat ion 

problem';" projection ("mathematics ) ; resolving power; higher harmonic; 
parameter estimation; pattern recognition; signal detection 
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decision; statistical method; detection 
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ABSTRACT: A structured network model for robot control and its learning 
algorithm are presented. The model is divided into two sub-networks 
connected each other with a 'short term memory to process time dependent 
data. We have developed the learning algorithm, pseudoimpedance 
control, which has damped oscillation characteristic to avoid the 
local minimum problem. To evaluate the network model and the learning 
algorithm, small mobile robots controlled by neural networks have been 
developed. They were taught variety of habits to play a 

cops-and-robbers game. Through training, the robots learned habits such 

as capture and escape .( author abst.) 
DESCRIPTORS: neural network; robot; computing control; brain; network 
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Sensor data compression for smart mechanical systems by artificial neural 
networks 

( Datenkompression an intelligenten Strukturen mittels kuenstlicher 
neuronaler Netze) 

Lilienblum, T; Weihua Zhang; Michaelis, B 



Ottfo-von-Guericke-Univ. Magdeburg, D 

Smart Mechanical Systems - Adaptronics . Proc of the Second Scientific 
Conf., Preprint Nr. 2, Otto-von-Guericke-Univ. Magdeburg, 18th - 19th March 
19971997 

Document type: Conference paper Language: English 
Record type: Abstract 

ABSTRACT : 

Zur. stabilisderenden Wdrkung.>der a*if. ' i-ntelligente^ -mechanische Systeme- 
einwirkenden Kraf tbelastungen und Schwingungen werden erfolgreich 
sensorische Systeme in diese Bauelemente integriert. Notwendig ist 
allerdings der gleichzeigtige Einsatz zahlreicher sensorischer Aktoren und 
ihre zu Regelungsbef ehlen fuehrenden Abtastsignale . In Verbindung mit 
zahlreichen Prinzipskizzen und f ormelmaessigen mathematischen 
Erlaeuterungen wird dargelegt, dass mit einem modif izierten assoziativen 
Speicher zur 3-D-Bildverarbeitung, die Kompression der Sensordaten eines 
solchen mechanischen Systems moeglich ist. Die Sensordaten werden in Form 
von Neuron in eine erste Schicht eingespeist. In der zweiten Schicht 
geschieht dann die erwuenschte Datenkompression . Die Daten werden dann als 
a-priori-Kennwerte (Wichtung) in einem Neural-Netzwerk gespeichert. Die 
Berechnungsdauer laesst sich durch den Einsatz eines Neuro - Computers 
reduzieren. Die Zweckmaessigkeit dieses Datenkompressionsver f ahrens konnte 
durch Simulation nachgewiesen werden. 

DESCRIPTORS: ARTIFICIAL NEURAL NETWORKS; DATA COMPRESSION; ADAPTIVE SYSTEM; 

OSC I LLATI ON ATT ENUAT I ON 
IDENTIFIERS: AKTIVE DAEMPFUNG; ADAPTRONIK; intelligente Struktur; 
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Attentional neurocomputing deserves more attention 
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JOURNAL CODE: EET LANGUAGE: English 

RECORD TYPE: Fulltext 
SECTION HEADING: SR 
WORD COUNT: 1305 
TEXT : 

CHAIRMAN HNC INC. 

... tools and human/machine interfaces. 

The Grey and Singer work may also lead to new approaches to 
neurocomputer hardware design- for example, neural chips consisting of 
large pools of analog phase-locked loop oscillators that interact 
dynamically with incoming data, central* rhythm generators ancTeacn 'other. 
Attentional neurocomputers may be better able to exploit the natural 
strengths of analog VLSI technology (e.g., the implementation of dynamic 
system elements, such as oscillators ) than current neurocomputing design 
approaches, which mainly use analog circuits to carry out arithmetic 
operations-functions that. . . 
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WORD COUNT: 242 LINE COUNT: 00018 
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TEXT: 

CHAOTIC NEUROCOMPUTER DEVELOPED 09/16/91 TOKYO, JAPAN, 1991 SEP 16 
(NB) — Professor Masayoshi Inoue of Kagashima University has developed a 
new concept in the area of neurocomputing. He uses what's called a Chaos 
Oscillator which employs the Chaos concept of order within disorder. 

Professor Inoue's neurocomputer , equipped with the "Chaos 
Oscillator , " is said to be able to recognize a variety of vague patterns 
and, he contends, is much closer in operation and reasoning to the human 
brain. The Chaos Oscillator is based on the mathematical formula 
"AX(l-X)." The numbers 0 through 1 are input as... 
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New technology learns Wall Street f s mindset, (artificial intelligence) 
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ABSTRACT: More than 200 firms, including start-ups and such giants as IBM 
and DEC, are researching neurocomputing . A neural network is a computing 
architecture based on the way the brain is believed to draw upon experience 
to infer ways of recognizing similarities and ignoring differences among 
like patterns. One of the advantages of neurocomputers is that they 
reduce the cost and time of software development. In an expert system, a 
human expert must write an algorithm - a step-by-step procedure - that 
tells the computer how to perform a desired task. In contrast, a 
neurocomputer comes up with an algorithm of its own. There is a wide 
consensus that neurocomputing is in its infancy and will not duplicate the 
functioning of the human brain for a very long time , if ever. 
CAPTIONS: An application that predicts the Standard & Poor's 500 index, 
(graph) 

SPECIAL FEATURES : illustration; graph 

COMPANY NAMES: Hecht-Nielsen Neurocomputer Corp . --Innovations ; Nestor 
Inc. — Contracts 

DESCRIPTORS: Artificial Intelligence; Wall Street; Neural Network; 
Forecasting; Brokerage Industry; Cooperative Agreements; Computer Design; 
Expert Systems; Algorithm Complexity;* Cost of Programming; Future- of • 
Computing 

SIC CODES: 7372 Prepackaged software 

TRADE NAMES: NeuralWorks Professional (Expert system development 

software) — Usage 
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...ABSTRACT: way the brain is believed to draw upon experience to infer 
ways of recognizing similarities and ignoring differences among like 
patterns. One of the advantages of neurocomputers is that they reduce the 
cost and time of software development. In an expert system, a human 
expert must write an algorithm - a step-by-step procedure - that tells the 
computer how to perform a desired task. In contrast, a neurocomputer 
comes up with an algorithm of its own. There is a wide consensus that 
neurocomputing is in its infancy and will not duplicate the functioning of 
the human brain for a very long time , if ever. 
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Object recognition opens the eyes of machine -vision systems . 

Williams, Tom 

Computer Design, v27, n9, p69(10) 
May 1, 1988 

ISSN: 0010-4566 LANGUAGE: ENGLISH RECORD TYPE: FULLTEXT; ABSTRACT 

WORD COUNT: 4236 LINE COUNT: 00331 

ABSTRACT: Automatic object recognition demands high computing power 
because the process is not a linear one easily defined by algorithms, 
global image association is required, and all received data needs to be 
simultaneously processed. Object recognition system vendors, though, are 
using new architectures and image-processing techniques to maximize the 
throughput and accuracy in image capture, image conditioning for 
recognition, extraction of relevant image data, and decision on image 
content. Evolving feature extraction systems include true gray-scale 
processing, edge and corner identification, and pattern matching. Modular 
systems provide adaptability to specific recognition applications. Neural 
networks offer great promise. -for fast object recognition/ as the * - 
architecture is oriented to the global processing and resolution required. 
Various products are described. 

CAPTIONS: Diagram of Hecht-Nielsen Neurocomputers' face recognizer, 
(chart) 

SPECIAL FEATURES: illustration; photograph; chart 

DESCRIPTORS: Image Processing; Object Recognition; Trends; New Technique; 
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... addresses, processing routines and lookup tables. 

One neural-network system is tthe Anza Plus from Hecht-Nielsen 
Neurocoinputers (San Diego, CA) . Anza Plus is an aT-compatible board that 
contains a reduced-instruction-set-type... 

...fast enough to seem connected in real time. The processing rate of such 
a neural network is measured in interconnects per second (IPS) , and the 
Anza Plus has a sustained rate of 6 million IPS... 
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CHAOTIC NEUROCOMPUTER DEVELOPED 09/16/91 
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Word Count: 226 
TEXT: 

TOKYO, JAPAN, 1991 SEP 16 (NB) ~ Professor Masayoshi Inoue of Kagashima 
University has developed a new concept in the area of neurocomputing . He 
uses what's called a Chaos Oscillator which employs the Chaos concept of 
order within disorder. 

Professor Inoue ! s neurocomputer / equipped with the "Chaos 
Oscillator , " is said to be able to recognize a variety of vague patterns 
and, he contends, is much closer in operation and reasoning to the human 
brain. The Chaos Oscillator is based on the mathematical formula 
M AX(1-X)." The numbers 0 through 1 are input as "X." Professor Inoue's 
model continuously compares two formulas to reach a conclusion, which is 
different from traditional, Von-Neuman-type computers, which deal with only 
1 and 0. 

Under Prof. Inoue's system, it is said more complicated patterns can 
be processed. Currently, Prof. Inoue is running the system with software, 
but he wants to incorporate the features onto a computer chip. The chip 
could play an important role in the development of programs based on 
advanced artificial intelligence. 

It is said Prof. Inoue has already applied for copyright in Japan, and 
is planning to apply for the same in the U.S. and the U.K. soon. Details of 
the study are expected to be introduced in the "Physics Letters A" of the 
Physics International Research in Europe by the end of this month. 
(Masayuki Miyazawa/ 19910916) 
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Professor Inoue's neurocomputer , equipped with the "Chaos 
Oscillator , " is said to be able to recognize a variety of vague patterns 
and, he contends, is much closer in operation and reasoning to the human 
brain. The Chaos Oscillator is based on the mathematical formula 
"AX(l-X)." The numbers 0 through 1 are input as... 
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from Silicon Graphics. The system, to be installed for testing in an F-18 
and a CH-53, will give pilots\real- time solid color graphics for mission 
presentations ... .Unisys, McLean, Va., is part of a team chosen by NASA to 
provide hardware, software ancksystems integration for the space station 
Freedom control center. Unisys ud.ll support p£i>me contractor Ford 
Aerospace. -Database management \sof tware from Informix will be included in 
three projects awarded recently byvthe Navy. /They are the Desktop Tactical 
Computer awarded to systems integrator C3 Inrc., the Tomahawk Planning 
System awarded to McDonnell Douglas >and LorAl, and the Metrology Automated 
Support for Recall and reporting (MEASURE) / won by Honeywell. 
COPYRIGHT 1989 BY PASHA PUBLICATIONS AnC. 

COPYRIGHT 1989 Pasha Publications \ncl 

COPYRIGHT 1999 Gale Group Xf 
PUBLISHER NAME: Pasha Publications, Inc/\ 

INDUSTRY NAMES: AERO (Aerospace and Defense); BUSN (Any type of business 
) ; CMPT (Computers and Office Automation) 

(USE FORMAT 7 ' FOR FULLTEXT) / \ 

TEXT: / \ 

...Defense Data Network. In the environment of placket-switching 
communication, Blacker permits controlled access\to information at 
different security levels based/on the user's need to know. .The Defense 
Advanced Research Projects Agency has given Hecht-Nie^son Neuro - 
computers , San Diego, a $250,000 contract for the second phase of 
development for the next-generation VLSI . . . 

...to be installed for testing in an F-18 and a CH-53, will give pilots 
real- time solid color graphics for mission presentations .... Unisys , 
McLean, Va . , is part of a team chosen by NASA... 
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ABSTRACT : 

Suntory and researchers from the Tokyo Inst of Technology have developed an 
automated quality control system for whisky, with equipment that can 
distinguish as many as 5 different brands. The aromatic molecule 
discrimination system, scheduled for commercial introduction by end-FY91, 
contains a heated air feeding unit that delivers samples to 8 sensors 
containing 8 types of membranes, including lipid, cholesterol, and 
cellulose, and quartz oscillators . The signals from the sensors are 
sent to a neurocomputer , which studies the output for 10,000 cycles, and 
having learned the patterns, can discern aromas of up to 5 brands of 
whiskey, comparable or superior to human performance. The Ministry of 
Education subsidized the 2-yr R&D project. 
PUBLISHER NAME: Japan External Trade Organization 
COMPANY NAMES: *Suntory Ltd. 

EVENT NAMES: *330 ' (Product information) ' ' 
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ABSTRACT: 

...developed an automated quality control system for whisky, with equipment 
that can distinguish as many as 5 different brands. The aromatic 
molecule discrimination system, scheduled for commercial introduction by 
end-FY9L, contains a heated air... 

...delivers samples to 8 sensors containing 8 types of membranes , ' including 
lipid, cholesterol, and cellulose, and quartz oscillators . The signals 
from the sensors are sent to a neuro computer , which studies the output 
for 10,000 cycles, and having learned the patterns, can discern aromas of 
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AISTs Electrotechnical laboratory Develops Image Recognition System 

Comline Computers Deceiver 7, 1988 p. 1 

An information scienoe R&D group at the Xgency of Industrial Science 
and Technology's (AIST) Electrotechnical Laboratory has developed an image 
recognition system which c\n learn to recognize its target images in one 
tenth the time of existing rieuro - computer systems. The system employs 
multivariate analysis instead \>f a neuyo-network system, and includes a 
video camera, application specific I'Cs /(ASICs) / and a personal * computer . 
The system recognizes an image inputs, from /the camera by comparing the image 
against a series of 25 optical mask pfetteins and analyzing the results, and 
on a test using 2,000 samples of different characters the system 

achieved 98.7% accuracy. \ 

The system can be constructed wi\h off-the-shelf LSIs, and does not 
require special devices, such as Y the optical elements used in 
neuro-computer systems, making the system n^pre commercially viable. 

COMLINE NEWS SERVICE, Sugetsu Building, \- 12-7 Kita-Aoyama, Minato-Ku, 
Tokyo 107, Japan. Telex 2428134 COMLN \ 
PRODUCT: *Image Analysis Equip (3662653) \ 
EVENT: ^Product Design & Development: (33) \ 
COUNTRY: * Japan (9JPN) / \ 

. . . developed an image recognition system which>can learn to recognize 
its target images in one tenth/ the time of exisrsing neuro - computer 
systems. The system employs/ multivariate analysis instead of a 
neuro-network system, and includes a video camera... \ 

. . . optical ' mask patterns and/ analyzing the results,' and on a test using 
2,000 samples of 37 different characters the systefo achieved 98.7% 
accuracy. ' • \ 

The system can be constructed with off-the-shelf LSIs... 
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HNC RELEASES ADVANCED 
RECOGNITION SYSTEMS 

News Release October 30, 

Hecht-Nielsen Neurocomputer Corp. (HNC) introduces AR/NET, the first 
fully- functional version of the Adaptive Resonance Network, a neural 
network architecture that is especially / adept at the unsupervised 
recognition and classification of arbitrary/ patterns,. In combination with 
the HNC ANZA Neurocomputing Coprocessor System (whicn^can be installed in 
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OS C I LLATARY NEUROCOMPUTERS WITH DYNAMIC CONNECTIVITY 
ORDINATEUR NEUROMIMETIQUE OSCILLATOIRE A CONNECTIVITE DYNAMIQUE 

Patent Applicant/Assignee: 

ARIZONA BOARD OF REGENTS a body corporate acting; on behalf of ARIZONA 

STATE UNIVERSITY, 
HOPPENSTEADT Frank C, 
IZHIKEVICH Eugene, 
Inventor (s) : 
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Claims 
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English Abstract 

A neurocomputer (50) comprises n oscillating processing elements (60A, 



-60B, 60C, 60D and 60E) that communicate through a common medium (7 0) so 
that there are required only n connective junctions (80A, 80B, 80C, 80D 
and 80E) . A rhythmic external forcing input (90) modulates the 
oscillatory frequency of the medium (70) which, in turn, is imparted to 
the n oscillators (60A, 60B, 60C, 60D and 60E) . Any two oscillators 
oscillating at different frequencies may communicate provided that 
input's power spectrum includes the frequency equal to the difference 
between the frequencies of the two oscillators in question. Thus, 
selective communication, or dynamic connectivity, between different 
neurocomputer oscillators occurs due to the frequency modulation of the 
medium (70) by external forcing. 

French Abstract 

Cet ordinateur neuromimetique (50) comprend n elements de traitement 
oscillants (60A, 60B, 60C, 60D et 60E) qui communiquent par 
1 ' intermediaire d'un support commun (7 0) de sorte que seulement n 
jonctions de connexion (80A, 80B, 80C, 80D et 80E) sont necessaires. 
L' entree d'une contrainte, exterieure, rythmique (90) module la frequence 
oscillatoire du support (70), laquelle est a son tour appliquee aux n 
oscillateurs (60A, 60B, 60C, 60D et 60E) . Deux oscillateurs quelconques, 
oscillant a des frequences differentes, peuvent communiquer pourvu que le 
spectre de la puissance d 1 entree comprenne la frequence egale a la 
difference entre les frequences des deux oscillateurs en question. Ainsi, 
il se produit une communication selective, ou une connectivite dynamique, 
entre differents oscillateurs de 1 1 ordinateur neuromimetique, par suite 
de la modulation de frequence du support (70) au moyen d l une contrainte 
exterieure . 
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Oscillatory neuro -computer for simulating oscillatory nature of brain 
neurons, has conductive medium coupled to connectors, which applies 
oscillatory signal to each oscillator via corresponding connector 

Patent Assignee: UNIV ARIZONA STATE (UYAR-N) 

Inventor: HOPPENSTEADT F C ; IZHIKEVICH E 

Number of Countries: 022 Number of Patents: 001 
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Abstract (Basic) : WO 200029970 Al 

NOVELTY - Several connectors (8.0A-80E) are. operably coupled with 
corresponding oscillators (60A-60E) . A conductive medium (70) is 
operably coupled to the connectors, simultaneously applies oscillatory 
signal to each oscillator via the connector. A'sourcing apparatus with 
rhythmic external forcing input (90) generating an oscillatory signal, 
is operably coupled with the medium. 

DETAILED DESCRIPTION - An INDEPENDENT CLAIM is also included for 
communication establishing method between two oscillator having 
different frequencies. 

USE - For simulating oscillatory nature of brain neurons. 

ADVANTAGE - Neuro-computer can act as a classical fully connected 
Hopfield network even when there are only interconnections. 

DESCRIPTION OF DRAWING (S) - The figure shows the schematic diagram 



of neural network having five neural processing elements. 
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Connectors (80A-80E) 
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Title: Diffusion bay simulation and its impact on the overall FAB 

performance : a simplified example 

Author(s): Collins, D.W.; Flores-Godoy, J. -J.; Tsakalis, K.S.; 
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Conference Title: 2003 IEEE International Symposium on Semiconductor 
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Conference Date: 30 Sept. -2 Oct. 2003 Conference Location: San Jose, 

CA, USA 

Medium: Also available on CD-ROM in PDF format 

Language: English Document Type: Conference Paper (PA) 

Treatment: Practical (P) ; Theoretical (T) 

Abstract: In this presentation the importance of simulation as a tool to 
understand the behavior of systems and their interactions is discussed. As 
an example, with the help of simulations of a diffusion bay it was possible 
'to compute the probability of failure of the equipment for different 
controllers. This information is then used in a simplified FAB to study the 
performance that can be expected for different scheduling policies. (16 
Refs) 
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Descriptors: failure analysis; integrated circuit manufacture; 
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Identifiers: diffusion bay simulation; overall FAB performance; 
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Journal: Neural Computation vol.15, no. 10 p. 2359-77 
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Publication Date: Jan. 2003 Country of Publication: USA 
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Material Identity Number: N733-2003-010 
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Abstract: A procedure is developed to probe the changes in the functional 
interactions among neurons in primary motor cortex of the monkey brain 
during adaptation. A monkey is trained to learn a new skill, moving its arm 



to reach a target under the influence of external perturbations. The spike 
trains of multiple neurons in the primary motor cortex are recorded 
simultaneously. We utilize the methodology of directed transfer function, 
derived from a class of linear stochastic models, to quantify the causal 
interactions between the neurons. We find that the coupling between the 
motor neurons tends to increase during the adaptation but return to the 
original level after the adaptation. Furthermore, there is evidence that 
adaptation tends to affect the topology of the neural network, despite the 
approximate conservation of the average coupling strength in the network 
before and after the adaptation. (40 Refs) 
Subfile: C 
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Title: Capacity of oscillatory associative -memory networks with error-free 
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Language: English Document Type: Journal Paper (JP) 
Treatment: Theoretical (T) 

Abstract: Networks of coupled periodic oscillators (similar to the 
Kuramoto model) have been proposed as models of associative memory. 
However, error-free retrieval states of such oscillatory networks are 
typically unstable, resulting in a - near zero: capacity. -This puts • the 
networks at disadvantage as compared with the classical Hopfield network. 
Here we propose a simple remedy for this undesirable property and show 
rigorously that the error-free capacity of our oscillatory, 
associative-memory networks can be made as high as that of the Hopfield 
network. They can thus not only provide insights into the origin of 
biological memory, but can also be potentially useful for applications in 
information science and engineering. (23 Refs) 
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Abstract: We consider a system of N phase oscillators having randomly 
distributed natural frequencies and diagonalizable interactions among the 
oscillators. We show that, in the limit of N to infinity , all solutions of 
such a system are incoherent with probability one for any strength of 
coupling, which implies that there is no sharp transition from incoherence 
to coherence as the coupling strength is increased, in striking contrast to 
Kuramoto's (special) oscillator system." (20 Refs) 
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Abstract: Small-world and scale-free networks are known to be more easily 
synchronized than regular lattices, which is usually attributed to the 
smaller network distance between oscillators. Surprisingly, we find that 
networks with a homogeneous distribution of connectivity are more 
synchronizable than heterogeneous ones, even though the average network 



distance is larger. We present numerical computations and analytical 
estimates on synchronizability of the network in terms of its heterogeneity 
parameters. Our results suggest that some degree of homogeneity is expected 
in naturally evolved structures, such as neural networks, where 
synchronizability is desirable. (42 Refs) 
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Abstract: Motivated by the practical consideration of the measurement of 
chaotic signals in experiments or the transmission of these signals through 
a physical medium, we investigate the effect of filtering on chaotic 
symbolic dynamics. We focus on the linear, time-invariant filters that are 
used frequently in many applications, and on the two quantities 
characterizing chaotic symbolic dynamics: topological entropy and bit-error 
rate. Theoretical consideration suggests that the topological entropy is 
invariant under filtering. Since computation of this entropy requires that 
the generating partition for defining the symbolic dynamics be known, in 
practical situations the computed entropy may change as a filtering 
parameter is changed. We find, through numerical computations and 
experiments with a chaotic electronic circuit, that with reasonable care 
the computed or measured entropy values can be preserved for a wide range 
of the filtering parameter. (32 Refs) - .■-*■* 
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Abstract: Efficiency in passage times is an important issue in designing 
networks, such as transportation or computer networks. The small-world 
networks have structures that yield high efficiency, while keeping the 
network highly clustered. We show that among all networks with the 
small-world structure, the most efficient ones have a "single center" node, 
from which all shortcuts are connected to uniformly distributed nodes over 
the network. The networks with several centers and a connected subnetwork 
of shortcuts are shown to be "almost" as efficient. Genetic-algorithm 
simulations further support our results. (14 Refs) 
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Abstract: The transition to phase synchronization in systems consisting * 
of a large number (N) of coupled nonlinear oscillators via the route of 



phase clustering (phase synchronization among subsets of oscillators) is 
investigated. We elucidate the mechanism for the merger of phase clusters 
and find an algebraic scaling between the critical coupling parameter 
required for phase synchronization and N. Our result implies that, in 
realistic situations, phase clustering may be more prevalent than full 
phase synchronization. (19 Refs) 
Subfile: A 

Descriptors: chaos; synchronisation 

Identifiers: phase clustering; transition to phase synchronization; 
coupled nonlinear oscillators; pha'se clusters; algebraic scaling; 'critical 
coupling parameter; chaos 

Class Codes: A0545 (Theory and models of chaotic systems) 

Copyright 2001, I EE 



1/5/9 (Item 9 from file: 2) 

DIALOG (R) File 2:INSPEC 

(c) 2004 Institution of Electrical Engineers. All rts . reserv. 

6969312 INSPEC Abstract Number: C2001-08-1290L-017 
Title: Oscillatory model of novelty detection 
Author (s): Borisyuk, R. ; Denham, M. ; Hoppensteadt , F. ; Kazanovich, Y. ; 
Vinogradova, O. 

Author Affiliation: Centre for Neural & Adaptive Syst, Plymouth, UK 
Journal: Network: Computation in Neural Systems vol.12, no . 1 p. 1-20 
Publisher: IOP Publishing, 

Publication Date: Feb. 2001 Country of Publication: UK 

CODEN: NTWREZ ISSN: 0954-898X 

SICI: 0954-898X (200102) 12 : 1L . 1 : OMND; 1-T 

Material Identity Number: N636-2001-001 

U.S. Copyright Clearance Center Code: 0954-898X/2001/ 010001+2 0$30 . 00 
Document Number: S0954-898X ( 01) 19653-3 

Language: English Document Type: Journal Paper (JP) 
Treatment: Theoretical (T) 

Abstract: A model of novelty detection is developed which is based on an 
oscillatory mechanism of memory formation and information processing. The 
frequency encoding of the input information and adaptation of natural 
frequencies of network oscillators to the frequency of the input signal are 
used as the mechanism of information storage. The resonance amplification 
of network activity is used as a recognition principle for familiar 
stimuli. Application of the model to novelty detection in the hippocampus 
is discussed. (34 Refs) 
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Abstract: We combine here two well-known and established concepts: 
microelectromechanical systems (MEMS) and neurocomputing . First, we 
consider MEMS oscillators having low amplitude activity and we derive a 
simple mathematical model that describes nonlinear phase-locking dynamics 
in them. Then, we investigate * a. theoretical possibility *of using MEMS 
oscillators to build an oscillatory neurocomputer having autocorrelative 
associative memory. The neurocomputer stores and retrieves complex 
oscillatory patterns in the form of synchronized states with appropriate 
phase relations between the oscillators. Thus, we show that MEMS alone can 
be used to build a sophisticated information processing system (U.S. 
provisional patent 60/178,654). (19 Refs) 
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Abstract: We describe here how a network of oscillators can perform 
neural computations. In particular, it shown how the connectivity within 
the network can be created to memorize data in terms of phase relations 
between synchronized states. The memorized states are extracted through 
correlation calculations. The influence of noise on the system is 
discussed. (4. Refs) 
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Abstract: We investigate here possible neurocomputational features of 
networks of laser oscillators. Our approach is similar to classical optical 
neurocomputing where artificial neurons are lasers, .and ..connection matrices 
are holographic media. However, we consider oscillatory neurons 
communicating via phases rather than amplitudes. Memorized patterns 
correspond to synchronized states where the neurons oscillate with equal 
frequencies and with prescribed phase relations. The mechanism of 
recognition is related to phase locking. (19 Refs) 
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Abstract: We propose a novel architecture of an oscillatory neural 
network that consists of phase-locked loop (PLL) circuits. It stores and 
retrieves complex oscillatory patterns as synchronized states with 
appropriate phase relations between neurons. (16 Refs) 
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Abstract: We describe a biologically inspired oscillatory neural network 
for memorizing temporal sequences of neural activity patterns. The neural 
network consists of interactive neural oscillators with all-to-all 
excitatory connections forced by a slow T-periodic signal. The dynamics of 
the network' aire viewed through a time" window with 'duration T . Tne * network 
memorizes binary patterns in terms of low and high activity of the 
corresponding oscillators. The learning rule is temporally asymmetric, and 
it takes into account the activity level of pre- and post- "synaptic" 
oscillators in two contiguous time windows. Recall by the network is fast: 
all memorized patterns of sequences are reproduced in the correct order 
during the same time window, but with a short time delay. The applicability 
of these results to studies of the hippocampus is discussed. (10 Refs) 
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Abstract: Spatial patterns of theta-rhythm activity in oscillatory models 
of the hippocampus are studied here using canonical models for both 
Hodgkin's class-1 and class-2 excitable neuronal systems. Dynamics of these 
models are studied in both the frequency domain, to determine phase-locking 
patterns, and in the time domain, to determine the amplitude responses 
resulting from phase-locking patterns. Computer simulations presented here 
demonstrate that phase deviations (timings) between inputs from the medial 
septum and the entorhinal cortex can create spatial patterns of 
theta-rhythm phase-locking. In this way, the authors show that the timing 
of inputs (not only their frequencies alone) can encode specific patterns 
of theta-rhythm activity. This study suggests new experiments to determine 
temporal and spatial synchronization. (45 Refs) 
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Abstract: The authors 1 study of thalamo-cortical systems suggests a new 
architecture for a neurocomputer that consists of oscillators having 



different frequencies and that are connected weakly via a common medium 
forced by an external input. Even though such oscillators are all 
interconnected homogeneously, the external input imposes a dynamic 
connectivity. The authors use Kuramoto ! s model to illustrate the idea and 
to prove that such a neurocomputer has oscillatory associative properties. 
Then, they discuss a general case. The advantage of such a neurocomputer is 
that it can be built using voltage controlled oscillators, optical 
oscillators, lasers, microelectromechanical systems, Josephson junctions, 
macromolecules , or oscillators of other kinds. (16 Refs) 
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Abstract: This Multiscale Integration of Manufacturing and Assembly 
Processes (MIMAP) demonstration project investigates the integration of two 
or more Thrust Area Groups (TAGs) by creating a flow of information and 
processes from two areas of research. The control theory research of two 
different controllers for diffusion furnaces used in semiconductor 
manufacturing which predict a specific window of machine failure, the 
mean-time-bef ore-failure (MTBF) . The objective is to increase yield, 
decrease cycle time, work-in-progress (WIP) and production costs. A global 
factory Minimum Inventory Variability Scheduling, policy (MIVP ( R) ). , - used to 
decrease cycle time and cycle time variance when compared to a 
first-in-first-out (FIFO) scheduling policy, was used to make' the 
comparisons between the two inner and outer loop controllers. The project's 
Cross Cutting Methodologies (CCMs) is ensured by the participation of 
faculty from three different colleges (LAS, CEAS, and CTAS) and two 
electrical engineering doctoral students. (20 Refs) 
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Treatment: Theoretical (T) 

Abstract: This paper describes some problems and investigations 
encountered when implementing new resource scheduling policies in a large 
semiconductor manufacturing- facility^ FAB) . The FAB '-described here uses a 
product release policy based on customer orders and a work-in-progress 
(WIP) chart. The scheduling of resource tools is done on a first in, first 
out (FIFO) basis on high speed tools and due date first (DDF) at bottleneck 
tools, except for high priority LOTS, called MAXIs. This presentation 
describes briefly the theory behind minimum inventory variability 
scheduling policies. A heuristic explanation of the minimum inventory 
variability for resource scheduling policies is given here. Finally a large 
semiconductor manufacturing facility is discussed in generic terms, 
including (sanitized) data collection. The results of the baseline output 
and historical data are compared to MIVSP. (13 Refs) 
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Abstract: The models and simulations here demonstrate that steady 
progressing waves are possible ^ in ^networks whose ^ elements % a f e near 
saddle-node on limit cycle bifurcations. Such networks arise in studies of 
the neocortex, and cortical waves of this kind have been observed. There 
are many areas of the brain whose function is based on propagation of 
activity, such as sound location in the nucleus laminaris. It is also 
possible that propagation of neural activity in the hippocampus is 
important for navigation and memorization. The authors have shown here that 
striated structures can support waves having different speeds which is not 
possible in single layer structures. These waves in striated structures 
moving at different speeds can create dynamic patterns with brief intervals 
of coincidence that are believed to play important roles in brain function. 
(15 Refs) 
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Abstract: This paper describes an implementation of the 1-Step Ahead 
Minimum Inventory Variability Resource Scheduling Policy (R) , in a large 
semiconductor facility (FAB) over the period from May, 1996, through 
January, 1997. The FAB described here uses a product release policy based 



on customer orders and a work-in-progress (WIP) chart. The scheduling of 
resource tools was done on a first in, first out (FIFO) basis on high speed 
tools and due date first (DDF) at bottleneck tools, except for high 
priority lots, called MAXIs. The FAB is discussed in generic terms 
(sanitized) because of the proprietary nature of the devices manufactured. 
Percentages of change in cycle time and output are presented. (18 Refs) 
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Abstract: It is a well-known fact that oscillatory networks can operate 
as Hopf ield-like neural networks, the only difference being that their 
attractors are limit cycles: one for each memorized pattern. The neuron 
activities are synchronized on the limit cycles, and neurons oscillate with 
fixed phase differences (time delays) . We prove that this property is a 
natural attribute of general weakly connected neural networks, and it is 
relatively independent of the equations that describe the network activity. 
In particular, we prove an analogue of the Cohen-Grossberg convergence 
theorem for oscillatory neural networks. (0 Refs) 
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Abstract: In this paper we do not discuss what a thalamo-cortical system 
modeled by a weakly connected oscillators can do, but we rather discuss 
what it cannot do. Interactions between any two cortical columns having 
oscillatory dynamics crucially depend on their frequencies. When the 
frequencies are different, the interactions are functionally insignificant 
(i.e., they average to zero) even when there are synaptic connections 
between the cortical columns. We say that there is a frequency gap that 
prevents interactions. When the frequencies are equal (or close) the 
oscillators interact via phase deviations. By adjusting the frequency of 
oscillations, each cortical column can turn on or off its connections with 
other columns. This mechanism resembles that of selective tuning in 
frequency modulated (FM) radios. A weak non-constant thalamic input can 
remove, the frequency gap and ..link .any two. oscillators -.provided the input is 
chosen appropriately. In the case of many cortical columns with 
incommensurable frequency gaps the thalamic forcing will be chaotic. By 
adjusting its temporal activity, the thalamus has complete control over 
information processing taking place in important parts of the cortex. (4 
Refs) 
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Abstract: A major drawback to most mathematical models in neuroscience is 
that they are either far away from reality or the results depend on the 
specific model. A promising alternative approach takes advantage of the 
fact that many complicated systems behave similarly when they operate near 
critical regimes, such as bifurcations. Using nonlinear dynamical system 
theory it "is -pos'sible to- prove • th'af all systems ''near certain critical 
regimes are governed by the same model, namely a canonical model. Briefly, 
a model is canonical if there is a continuous change of variables that 
transforms any other model that is near the same critical regime to this 
one. Thus, the question of plausibility of a mathematical model is replaced 
by the question of plausibility of the critical regime. Another advantage 
of the canonical model approach to neuroscience is that rigorous derivation 
of the models is possible even when only partial information is known about 
anatomy and physiology of brain structures. Then, studying canonical models 
can reveal some general laws and restrictions. In particular, one can 
determine what certain brain structures cannot accomplish regardless of 
their mathematical model. Since the existence of such canonical models 
might sound too good to be true, we present a list of some of them for 
weakly connected neural networks. Studying such canonical models provides 
information about all weakly connected neural networks, even those that 
have not been discovered yet. (7 Refs) 
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Abstract: By a weakly connected neural network we mean a dynamical system 
of the form x/sub i/=f/sub il (x, lambda )+ epsilon g/sub i/(x/sub 1/, 
x/sub n/, tau , epsilon ), Kor=y<or=n, (1) where the multidimensional 
variables x/sub i/, lambda and r are the activity of i/sup th/ neuron, an 
internal tuning (bifurcation) parameter and an external input parameter, 
respectively; 0<or= epsilon «1 is a small parameter, f/sub il describes 



the dynamics of the i/sup th/ neuron and g/sub i/ describes network 
properties. If (1) is near a hyperbolic equilibrium, then its dynamics are 
qualitatively like a linear system. In order to exhibit any interesting 
nonlinear behavior, the neurons must be near their thresholds. In this 
case, (1) can exhibit multiple bifurcations, and it can be significantly 
simplified by reduction to a canonical jnodel, for example to x/sub i/=r/sub 
i/+b/sub i/x/sub i/+x/sub i//sup 2/+ Sigma /sub j=l//sup k/c/sub ij/x/sub 
j/, x/sub i/ in R for a multiple fold bifurcation, or to x/sub i/=r/sub 
i/+b/sub i/x/sub i/+or-x/sub i//sup 3/+ Sigma /sub j=l//sup k/c/sub 
ij/x/sub j/, x/sub i/ in R for a multiple cusp singularity. (0 Refs) 
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Abstract: For pt . I see ibid., vol. 75, no. 2, p. 117-27 (1996). This is 
the second of two articles devoted to analyzing the relationship between 
synaptic organizations (anatomy) and dynamical properties (function) of 
networks of neural oscillators near multiple supercritical Andronov-Hopf 
bifurcation points. Here the authors analyze learning processes in such 
networks. Regarding learning dynamics, the authors assume: (1) learning is 
local (i.e. synaptic modification depends on pre- and postsynaptic neurons 
but not on others); (2) synapses modify slowly relative to characteristic 
neuron response times; and (3) in the absence ^ of either^ pre- or 
postsynaptic activity, the synapse weakens (forgets). The authors 1 major 
goal is to analyze all synaptic organizations of oscillatory neural 
networks that can memorize and retrieve phase information or time delays. 
They show that such networks have the following attributes: (1) the rate of 
synaptic plasticity connected with learning is determined locally by the 
presynaptic neurons; (2) the excitatory neurons must be long-axon relay 
neurons capable of forming distant connections with other excitatory and 
inhibitory neurons; and (3) if inhibitory neurons have long axons, then the 
network can learn, passively forget and actively unlearn information by 
adjusting synaptic plasticity rates. (8 Refs) 
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Abstract:'" The authors 'study weakly connected networks of neural 
oscillators near multiple Andronov-Hopf bifurcation points. They analyze 
relationships between synaptic organizations (anatomy) of the networks and 
their dynamical properties (function). The authors' principal assumptions 
are: (1) each neural oscillator comprises two populations of neurons: 
excitatory and inhibitory ones; (2) activity of each population of neurons 
is described by a scalar (one-dimensional) variable; (3) each neural 
oscillator is near a nondegenerate supercritical Andronov-Hopf bifurcation 
point; (4) the synaptic connections between the neural oscillators are 
weak. All neural networks satisfying these hypotheses are governed by the 
same dynamical system, which the authors call the canonical model. Studying 
the canonical model shows that: (1) A neural oscillator can communicate 
only with those oscillators which have roughly the same natural frequency. 
That is, synaptic connections between a pair of oscillators having 
different natural frequencies are functionally insignificant. (2) Two 
neural oscillators having the same natural frequencies might not 
communicate if the connections between them are from among a class of 
pathological synaptic configurations. In both cases the anatomical presence 
of synaptic connections between neural oscillators does not necessarily 
guarantee that the connections are functionally significant. (3) There can 
be substantial ' pha'se differences "(time delay's') *' between tne* neural 
oscillators, which result from the synaptic organization of the network, 
not from the transmission delays. Using the canonical model the authors can 
illustrate self-ignition and autonomous quiescence (oscillator death) 
phenomena. That is, a network of passive elements can exhibit active 
properties and vice versa. The authors also study how Dale's principle 
affects dynamics of the networks, in particular, the phase differences that 
the network can reproduce. The authors present a complete classification of 
all possible synaptic organizations from this point of view. The theory 
developed here casts some light on relations between synaptic organization 
and functional properties of oscillatory networks. The major advantage of 
the authors 1 approach is that they obtain results about all networks of 
neural oscillators, including the real brain. The major drawback is that 
the authors ' findings are valid only when the brain operates near a 
critical regime, viz. for a multiple Andronov-Hopf bifurcation. (32 Refs) 
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Abstract: A network of voltage-controlled oscillator neuron (VCON) models 
is formed by a gradient of negative feedbacks, similar to a gradient of 
inhibitory connections in a neural system. A network of this kind is called 
a prism. Computer simulation notation to describe all the variables in the 
network is introduced. Computer simulation of the network shows that a 
smooth gradient of output frequencies results from monochromatic 
stimulation. (2 Refs) 
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Abstract: The authors show how the form of the perturbation approximation 
for the solution of stiff systems of ordinary differential equations with 
an identifiable , small parameter can bemused to generate -associated nonstiff 
or relaxed equations. Solutions of these relaxed equations are easily 
calculated and appropriate combinations of these solutions furnish 



numerical approximations to the original stiff problem. Variations of this 
method are applied to two classes of initial value problems: those with 
highly oscillatory solutions and those with rapidly equilibrating 
solutions. (7 Refs) 
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Abstract:- A mult-L-ti-me perturbation algorithm is* derived to -study ■ certain 
synchronous control systems and nonlinear renewal processes. The new 
methods presented study singular perturbation problems or nonlinear 
Volterra integro-dif f erential equations whose kernels are close to the 
Dirac delta function. The algorithm provides useful approximate solutions 
to weakly filtered phase-locked loop circuits and to nonlinear renewal 
equations describing certain population dynamics. (11 Refs) 

Subfile: B C 

Descriptors: control theory; demography; ecology; feedback; 
integro-dif f erential equations; nonlinear differential equations; 
perturbation techniques; phase-locked loops 

Identifiers: multitime perturbation algorithm; nonlinear Volterra 
integrodiff erential equations; algorithm; weakly filtered synchronous 
control systems; nonlinear renewal processes; certain synchronous control 
systems; nonlinear renewal processes; Dirac delta function; filtered 
phase-locked loop circuits; nonlinear renewal equations; population 
dynamics 

Class Codes: B0220 (Analysis); B1160 (Nonlinear network analysis and 
design); B1250 (Modulators, demodulators, discriminators and mixers); 
C1120 (Analysis); C1290B (Natural resources and ecology); C1300 (Control 
theory) • - • ■ - ■ - , < ..... . . . „ . 



1/5/30 (Item 30 from file: 2) 

DIALOG (R) File 2: INSPEC 

(c) 2004 Institution of Electrical Engineers. All rts. reserv. 

01723009 INSPEC Abstract Number: C81025330 
Title: Threshold analysis of a drug use epidemic model 
Author(s): Hoppensteadt, F.C. ; Murray, J.D. 

Author Affiliation: Dept. of Math., Univ. of Utah, Salt Lake City, UT, 
USA 

Journal: Mathematical Biosciences vol.53, no. 1-2 p. 79-87 
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Treatment: Theoretical (T) 

Abstract: A model of an individual's response to a drug is formulated. 
This is based on interaction between applied dosage of the drug and active 
and inactive receptor sites in the body. A population of such individuals 
is considered. These are divided among nonusers ( susceptibles ) , active 
users, and individuals removed through treatment and cure, disenchantment, 
etc. Recruitment of susceptibles to active drug use is assumed to occur 
through contact with active users, and the effectiveness of recruitment (or 
infectiousness) depends on the age of a user. The model is based solely on 
the user's response to the drug, and it is shown that when a certain 
combination., of \ susceptible population^size, individual- -susceptibility, , and 
infectiousness does not exceed a critical threshold value, there will be 
only few users. But when the threshold value is exceeded, an epidemic of 
drug use ensues. (2 Refs) 
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Abstract: The sine-Gordon equation and systems of discrete approximations 
to it which are respectively a model of the Josephson junction and models 
of coupled-point Josephson junctions are studied. This is done in the 
so-called current-driven case. The voltage response of these devices is the 
average of the time derivative of the solution of these equations and the 
authors demonstrate the existence of running periodic solutions for which 
the average exists. Static solutions are also studied. These together with 
the running solutions explain the multiple— valuedness in the response of a 
Josephson junction in several cases. The stability of the various solutions 
is described in some of the cases. Numerical results are displayed which 
give details of structure of solution types in the case of a single point 
junction and of the one-dimensional distributed junction. (10 Refs) 
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Abstract: A van der Pol relaxation oscillator that is subjected to 
external sinusoidal forcing can exhibit stable and unstable periodic and 
almost periodic responses. For some forcing amplitudes it even happens that 
two stable subharmonics having different periods may coexist. The stable 
responses of such forced oscillators are investigated here. By numerically 
computing the rotation number of stable oscillations for various values of 
the forcing amplitude and oscillator tuning, descriptions are obtained of 
regions of* phase locking, successive bifurcation bf" stable subharmonic ' and 
almost periodic oscillations, and overlap regions where two distinct stable 
oscillations can coexist. (16 Refs) 
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Abstract: Exploitation of a fish population distributed in a habitat 
bounded on one side by a breeding ground and on the other by an unfavorable 
environment is studied. The population's dynamics on the breeding grounds 
are assumed to be described by a simple depensatory function. The effects 
on this population of harvesting by a fishing fleet are determined as 
maximum effort* arid harvesting qubtas^are varied/ In "'particular , * threshold 
values for these parameters are derived beyond which an open access fishery 
collapses. Competition with an external fleet and dynamic optimisation are 
discussed briefly. (2 Refs) 
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Abstract: It is shown here that certain systems of nonlinear (parabolic) 
reaction-diffusion equations have solutions which are approximated by 
oscillatory functions in the form R( xi -c tau )P(t*) where P(t*) 
represents a sinusoidal oscillation on a fast time scale t* and R( xi -c 
tau ) represents a slowly-varying modulating amplitude on slow space ( xi ) 
and slow time ( tau ) scales. Such solutions describe phenomena in chemical 
reactors, chemical and biological reactions, and in other media where a 
stable oscillation at each point (or site) undergoes a slow amplitude 
change due to diffusion. (19 Refs) 
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Abstract: Studies the large time behavior of solutions x(t) to the 
nonlinear equation x(t)= gamma integral /sub t-l//sup t/x( eta )(l-x( eta 
,))d eta which, satisfy 0<pr=x ( t>) <pr=J . . It is. Jcnown- -that .such, solutions 
approach constants as t to infinity . The authors investigate the way in 
which the limit is attained. For 0<or= gamma <or=3, the convergence is 
dominated by iterates of a quadratic polynomial. For gamma »1, a multitime 
perturbation expansion is derived for the solution. The leading coefficient 
satisfies Burgers 1 equation with periodic boundary conditions and the 
higher order coefficients are determined as solutions of a linear heat 
equation with periodic boundary conditions. (3 Refs) 
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Abstract: The model presented describes the ■ spread of an infection in a 
population by keeping track of the chronological ages of the participants 
as well as their 1 class ages 1 , (i.e. the length of time since entering 
their present state) . The reasoning behind this model is similar to that 
used in the equation of age dependent population growth. (7 Refs) 
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Conference Title: Proceedings of the 8th annual Allerton conference on 
circuit and system theory p. 303-9 
Publisher: IEEE, New York, NY, USA 

Publication Date: 1970 Country of Publication: USA xii+657 pp. 

Conference Sponsor: Univ. Illinois, Dept. Elect. Engng . , Co-ordinated 
Sci. Lab., IEEE, Groups on Circuit Theory and Automatic Control 

Conference Date: 7-9 Oct. 1970 Conference Location: Monticello, IL, 
USA 

Language: English Document Type: Conference Paper (PA) 
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Abstract: Systems of the form dx/dt=f ( t, x, y, epsilon ), epsilon 
dy/dt=g ( t , x, y, epsilon ) are considered. This paper studies conditions 
under which the solution of the initial value problem associated with this 
system can be written as (x( tau , epsilon ), y*(t, epsilon )) + (X( tau , 
epsilon ), Y( tau , epsilon j) where (x*,y*) is a solution of the system 
which is smooth at epsilon =0 and (X,Y) is an exponentially decaying 
function of tau =t . epsilon as tau to infinity . 
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sinusoidal forcing can exhibit stable and unstable periodic and almost 
periodic responses. For some forcing amplitudes it- even happens that two 
stable subharmonics having different periods may coexist. The stable 
responses of such forced oscillators are investigated. By numerically 
computing the rotation number of stable oscillations for various values of 
the forcing amplitude and oscillator tuning, descriptions are obtained of 
regions of phase locking, successive bifurcation of stable subharmonic and 
almost periodic oscillations, and overlap regions where two distinct stable 
oscillations can coexist. 
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asymptotically .stable, .as t neajrs. iafinity, and a- domain of stability is 
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Title: Pattern recognition via synchronization in phase-locked loop 
neural networks 

Author: Hoppensteadt, Frank C. ; Izhikevich, Eugene M. 

Corporate Source: Arizona State Univ, Tempe, AZ, USA 

Source: IEEE Transactions on Neural Networks v 11 n 3 May 2000. p 734-738 



Publication Year: 2000 
' CODEN: ITNNEP* ISSN: 1045-9227 1 * ' i-. ..... 

Language: English 

Document Type: JA; (Journal Article) Treatment: A; (Applications); X; 
(Experimental) 

Journal Announcement: 0010W2 

Abstract: We propose a novel architecture of an oscillatory neural 
network that consists of phase-locked loop (PLL) circuits. It stores and 
retrieves complex oscillatory patterns as synchronized states with 
appropriate phase relations between neurons. (Author abstract) 15 Ref s . 

Descriptors: *Neural networks; Phase locked loops; Variable frequency 
oscillators; Pattern recognition; Synchronization; Phase shifters; Natural 
frequencies 

Identifiers: Phase locked loop neural networks; Brain rhythms; 
Oscillatory associative memory 
Classification Codes: 

723.4 (Artificial Intelligence); 713.5 (Other Electronic Circuits); 
713.2 (Oscillators); 723.5 (Computer Applications) 
723 (Computer Software) ; 713 (Electronic Circuits) 
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Title: Oscillatory model of the hippocampal memory 
Author: Boris yuk, Roman; Hoppensteadt , Frank 

Corporate Source: Univ of Plymouth, Plymouth, UK 

Conference Title: International Joint Conference on Neural Networks 
(IJCNN* 99) 

Conference Location: Washington, DC, USA Conference Date: 
19990710-19990716 

Source: Proceedings of the International Joint Conference on Neural 
Networks v 1 1999. IEEE, USA. p 42-45 
Publication Year: 1999 
CODEN: 850FAE 
Language: English 

Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 
Journal Announcement: 0007W1 

Abstractive describe "a bidlogically'lnspired oscillatory neural network 
for memorizing temporal sequences of neural activity patterns. The neural 
network consists of interactive neural oscillators with all-to-all' 
excitatory connections forced by a slow T-periodic signal. The dynamics of 
the network are viewed through a time window with duration T. The network 
memorizes binary patterns in terms of low and high activity of the 
corresponding oscillators. The learning rule is temporally asymmetric, and 
it takes into account the activity level of pre- and post- % synaptic 1 
oscillators in two contiguous time windows. Recall by the network is fast: 
All memorized patterns of sequences are reproduced in the correct order 
during the same time window, but with a short time delay. The applicability 
of these results to studies of the hippocampus is discussed. (Author 
abstract) 10 Refs. 

Descriptors: *Neural networks; Pattern recognition; Learning systems; 
Knowledge based systems; Mathematical models 

Identifiers: Hippocampal memory; Oscillatory neural networks; Excitatory 
connections 

Classification Codes: 

723.4.1 (Expert Systems) 

461.1 (Biomedical Engineering); 723.4 (Artificial Intelligence); 461.4 
(Human Engineering) 

461 (Biotechnology); 723 (Computer Software); 921 (Applied 
Mathematics) 
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Title: Implementation of Minimum Inventory Variability Scheduling 1-Step 
Ahead Policy in a large semiconductor manufacturing facility 

Author: Collins, Donald W.; Williams, Ken; Hoppensteadt , Frank C. 

Corporate Source: Arizona State Univ East, Mesa, AZ, USA 

Conference Title: Proceedings of the 1997 IEEE 6th International 
Conference on Emerging Technologies and Factory Automation, ETFA 1 97 

Conference Location: Los Angeles, CA, USA Conference Date: 
19970909-19970912 

Sponsor: IEEE 

E.I. Conference No.: 47152 

Source: IEEE Symposium on Emerging Technologies & Factory Automation, 
ETFA 1997. IEEE, Piscataway, NJ, USA, 97TH8314 . p 497-504 
Publication Year: 1997 
CODEN: 8 5 ROAM 

Language :• -English - ■ * . 

Document Type: CA; (Conference Article) Treatment: G; (General Review); 
T; (Theoretical) 

Journal Announcement: 9712W4 

Abstract: This paper describes an implementation of the 1-Step Ahead 
Minimum Inventory Variability Resource Scheduling Policy, in a large 
semiconductor facility (FAB) over the period from May, 1996, through 
January, 1997. The FAB described here uses a product release policy based 
on customer orders and a work-in-progress (WIP) chart. The scheduling of 
resource tools was done on a first in, first out (FIFO) basis on high speed 
tools and due date first (DDF) at bottleneck tools, except for high 
priority lots, called MAXI 1 s . The FAB is discussed in generic terms 
(sanitized) because of the proprietary nature of the devices manufactured. 
Percentages of change in cycle time and output are presented. (Author 
abstract) 18 Ref s . 

Descriptors: ^Production control; Inventory control; Semiconductor device 
manufacture; Scheduling 

Identifiers: Inventory variability scheduling 
Classification Codes: 

913.2 (Production Control); 911.3 (Inventory Control); 714.2 
(Semiconductor Devi-ces* & Integrated Circuits) <■«■■- v--- - « - 

913 (Production Planning & Control); 911 (Industrial Economics); 714 
(Electronic Components) 
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Title: Investigation of minimum inventory variability scheduling policies 
in a large semiconductor manufacturing facility 

Author: Collins, Donald W.; Hoppensteadt, Frank C. 

Corporate Source: Arizona State Univ East, Mesa, AZ, USA 

Conference Title: Proceedings of the 1997 American Control Conference. 
Part 3 (of 6) 

Conference Location: Albuquerque, NM, USA Conference Date: 
19970604-19970606 

Sponsor: IEEE 

I. -Conference No. 4-7 028 - ■* < ■ 

Source: Proceedings of the American Control Conference v 3 1997. IEEE, 
Piscataway, NJ, USA, 97CH3604 1 . p 1924-1928 

Publication Year: 1997 

CODEN: PRACEO ISSN: 0743-1619 

Language: English 

Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 



-Journal Announcement: 9711W3 

Abstract: This paper describes some problems and investigations 
encountered when implementing new resource scheduling policies in a large 
semiconductor manufacturing facility (FAB) . The FAB described here uses a 
product release policy based on customer orders and a work-in-progress 
(WIP) chart. The scheduling of resource tools is done on a first in, first 
out (FIFO) basis on high speed tools and due date first (DDF) at bottleneck 
tools, except for high priority LOTS, called MAXI's. Minimum Inventory 
Variability Scheduling Policies (MIVSP) were introduced in February 1996. 
The development of a simulation model representing the FAB and a partial 
implementation of MIVSP were completed over the period from May, 1996, 
through January, 1997. This presentation describes briefly the theory 
behind MIVSP. A heuristic explanation of the minimum inventory variability 
for resource scheduling policies is given here. Finally a large 
semiconductor manufacturing facility is discussed in generic terms, 
including (sanitized) data collection. The results of the baseline output 
and historical data are compared to MIVSP. (Author abstract) 13 Refs. 

Descriptors: ^Semiconductor device manufacture; Production control; 
Resource allocation; Computer simulation; Inventory control; Scheduling 

Identifiers: minimum, inventory variability scheduling, ^policies . ^MIVSP ) ; 
Work in progress (WIP) chart 

Classification Codes: 

714.2 (Semiconductor Devices & Integrated Circuits); 913.2 (Production 
Control); 912.2 (Management); 723.5 (Computer Applications) 

714 (Electronic Components); 913 (Production Planning & Control); 912 
(Industrial Engineering & Management); 723 (Computer Software) 
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Title: Associative memory of weakly connected oscillators 
Author: Hoppensteadt , Frank C. ; Izhikevich, Eugene M. 
Corporate Source: Arizona State Univ, Tempe, AZ, USA 

Conference Title: Proceedings of the 1997 IEEE International Conference 
.on Neural Networks.. Pa,rt 2 (o.f 4) „ t ^ ( ....... . . „ 

Conference Location: Houston, TX, USA ' Conference Date: 

19970609-19970612 

Sponsor: IEEE 

E.I. Conference No.: 46924 

Source: IEEE International Conference on Neural Networks - Conference 
Proceedings v 2 1997. IEEE, Piscataway, NJ, USA, 97CB36109 . p 1135-1138 
Publication Year: 1997 
CODEN: ICNNF9 
Language: English 

Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 
Journal Announcement: 9710W4 

Abstract: It is a well-known fact that oscillatory networks can operate 
as Hopf ield-like neural networks, the only difference being that their 
attractors are limit cycles: one for each memorized pattern. The neuron 
activities are synchronized on the limit cycles, and neurons oscillate with 
fixed phase differences (time delays) . We prove that this property is a 
natural attribute of general weakly connected neural networks, and it is 
relatively independent of the equations that describe the network activity. 
In particular, we prove an analogue of the Cohen-Grossberg convergence 
.theorem for. .oscillatory neural networks,. .(Author , ajDstr.act) 8Befs. . 

Descriptors: *Neural networks; Associative storage; Bifurcation 
(mathematics); Convergence of numerical methods; Vectors 

Identifiers: Weakly connected oscillators; Hopfield neural networks 

Classification Codes: 

723.4 (Artificial Intelligence); 722.1 (Data Storage, Equipment & 
Techniques); 921.6 (Numerical Methods); 921.1 (Algebra) 

723 (Computer Software); 722 (Computer Hardware); 921 (Applied 
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Title: Thai amo -cortical interactions modeled by forced weakly connected 
oscillatory networks 

Author: Hoppensteadt, Frank C. ; Izhikevich, Eugene M. 

Corporate Source: Arizona State Univ, Tempe, AZ, USA 

Conference, Title.: ^Proceedings .pf t thje 1997 IE£E t International . Conference 
on Neural Networks. Part 1 (of 4) 

Conference Location: Houston, TX, USA Conference Date: 

19970609-19970612 

Sponsor: IEEE 

E.I. Conference No.: 46924 

Source: IEEE International Conference on Neural Networks - Conference 
Proceedings v 1 1997. IEEE, Piscataway, NJ, USA, 97CB36109 . p 328-331 
Publication Year: 1997 
CODEN: ICNNF9 
Language: English 

Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 
Journal Announcement: 9710W4 

Abstract: In this paper we do not discuss what a thalamo-cortical system 
modeled by a weakly connected oscillators can do, but we rather discuss 
what it cannot do. Interactions between any two cortical columns having 
oscillatory dynamics crucially depend on their frequencies. When the 
frequencies are different, the interactions are functionally insignificant 
(i.e., they average to zero) even when there are synaptic connections 
between the cortical columns. We say that there is a frequency gap that 
prevents in.t.eractip.ns When the frequencies are equal v (.or close) the 
oscillators interact via phase deviations. By adjusting the frequency of 
oscillations, each cortical column can turn on or off its connections with 
other columns. This mechanism resembles that of selective tuning in 
Frequency Modulated (FM) radios. A weak non-constant thalamic input can 
remove the frequency gap and link any two oscillators provided the input is 
chosen appropriately. In the case of many cortical columns with 
incommensurable frequency gaps the thalamic forcing will be chaotic. By 
adjusting its temporal activity, the thalamus has complete control over 
information processing taking place in important parts of the cortex. 
(Author abstract) 4 Refs. 

Descriptors: *Neural networks; Mathematical models; Oscillators 
(electronic); Natural frequencies; Oscillations; Chaos theory; Data 
processing; Frequency modulation 

Identifiers: Thalamo cortical interactions; Weakly connected neural 
networks (WCNN) 

Classification Codes: 
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Title: Canonical models for mathematical neuroscience 
Author: Hoppensteadt, Frank C. ; Izhikevich, Eugene M. 
Corporate Source: Arizona State Univ, Tempe, AZ, USA 

Conference Title: Proceedings of the 1997 IEEE International Conference 



on Neural Networks. Part 1 (of 4) 

Conference Location: Houston, TX, USA Conference Date: 

19970609-19970612 

Sponsor: IEEE 

E.I. Conference No.: 46924 

Source: IEEE International Conference on Neural Networks - Conference 
Proceedings v 1 1997. IEEE, Piscataway, NJ, USA, 97CB36109 . p 324-327 
Publication Year: 1997 
CODEN: ICNNF9 
Language: English 

Document Type: CA; (Conference Article) Treatment: T; (Theoretical) 
Journal Announcement: 9710W4 

Abstract :* A majo'r drawback 'to most ' mathematical mode' Is in neuroscience is 
that they are either far away from reality or the results depend on the 
specific model. A promising alternative approach takes advantage of the 
fact that many complicated systems behave similarly when they operate near 
critical regimes, such as bifurcations. Using non-linear dynamical system 
theory it is possible to prove that all systems near certain critical 
regimes are governed by the same model, namely a canonical model. Briefly, 
a model is canonical if there is a continuous change of variables that 
transforms any other model that is near the same critical regime to this 
one. Thus, the question of plausibility of a mathematical model is replaced 
by the question of plausibility of the critical regime. Another advantage 
of the canonical model approach to neuroscience is that rigorous derivation 
of the models is possible even when only partial information is known about 
anatomy and physiology of brain structures. Then, studying canonical models 
can reveal some general laws and restrictions. In particular, one can 
determine what certain brain structures cannot accomplish regardless of 
their mathematical model. Since the existence of such canonical models 
might sound too good to be true, we present a list of some of them for 
weakly connected neural networks. Studying such canonical models provides 
information about all weakly connected neural networks, even those that 
have not been discovered yet." (Author abstract) 1 "ReEsV ' 

Descriptors: *Neural networks; Mathematical models; System theory; 
Neurophysiology; Bifurcation (mathematics); Brain models 

Identifiers: Mathematical neuroscience; Canonical models 
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Title: Singular perturbation solutions of noisy systems 
Author: Hoppensteadt , Frank C. 

Corporate' "Source':' Mi'chigan "State Univ, East Lansing,'* Ml, USA 

Source: SIAM Journal on Applied Mathematics v 55 n 2 Apr 1995. p 544-551 

Publication Year: 1995 

CODEN: SMJMAP ISSN: 0036-1399 

Language: English 

Document Type: JA; (Journal Article) Treatment: T; (Theoretical) 
Journal Announcement: 9506W3 

Abstract: Recent work on singular perturbation solutions that persist in 
the presence of noise is described. Two different settings are considered: 
small deviation theory in quasi-static problems, where there are small 
amplitude but highly irregular perturbations, and averaging problems where 
there are ergodic stochastic perturbations. In the first case, it is shown 
that quasi-static approximations can be valid when the underlying problem 
experiences small deviation perturbations in problems that are stable under 
persistent disturbances. In the second, averaging principles are described 



for certain dynamical systems in Hilbert spaces that include applications 
to a wide variety of initial-boundary value problems for partial 
differential equations and for Volterra integral equations. These methods 
are applied here to four problems arising in applications. (Author 
abstract) 7 Refs. 

Descriptors: ^Perturbation techniques; Integral equations; Approximation 
theory; Problem solving; Convergence of numerical methods; Partial 
differential equations; Boundary value problems; Vectors 

Identifiers: Singular perturbation methods; Stochastic integral equations 
; Noisy systems; Averaging problems; Hilbert spaces 
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Title: Particle method for population waves 
Author: Chiu, Chichia; Hoppensteadt, Frank C. 

Corporate Source: Michigan State Univ, East Lansing, MI, USA 

Source: SIAM Journal on Applied Mathematics v 54 n 2 Apr 1994. p 466-477 

Publication Year: 1994 

CODEN: SMJMAP ISSN: 0036-1399 

Language: English 

Document Type: JA; (Journal Article) Treatment: T; (Theoretical) 
Journal Announcement: 94 09W1 

Abstract: Phase models are useful for studying synchronization of 
bacterial cell culture growth and other* biological* everrts associated with 
cell cycles. This paper considers a model that allows the growth rate of 
cells to change at phases of cell cycle. In this paper, a particle method 
is derived for solving the weak formulation of this model, convergence of 
the particle method is proved. 

Descriptors: *Cell culture; Growth kinetics; Mathematical models; Error 
analysis 

Identifiers: Cell populations; Vortex methods; Particle methods; Cell 
cycle 
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Title: Computer simulation > of a neural prism. - • * - 

Author: Hoppensteadt, F. C. 

Corporate Source: Dept of Math, Michigan State Univ, E Lansing, MI, USA 
Conference Title: Proceedings of the 32nd Midwest Symposium on Circuits 
and Systems 

Conference Location: Champaign, IL, USA Conference Date: 19890814 
Sponsor: IEEE Circuits and Systems Soc; IEEE Acoustics, Speech, and 
Signal Processing Soc; IEEE Control Soc; IEEE Education Soc; IEEE 
Industrial Electronics Soc; IEEE Central Illinois Section 
E.I. Conference No.: 13999 

Source: Midwest Symposium on Circuits and Systems. Publ by IEEE, IEEE 
Service Center, Piscataway, NJ, USA (IEEE cat n 89CH2785-4). p 238-239 
Publication Year: 1990 
CODEN: MSCSDL 
Language: English 

Document Type: PA; (Conference Paper) Treatment: X; (Experimental) 



Journal Announcement: 9102 

Abstract: A network of voltage-controlled oscillator neuron (VCON) models 
is formed by a gradient of negative feedbacks, similar to a gradient of 
inhibitory connections in a neural system. A network of this kind is called 
a prism. Computer simulation notation to describe all the variables in the 
network is introduced. Computer simulation of the network shows that a 
smooth gradient of output frequencies results from monochromatic 
stimulation. 2 Ref s . 

Descriptors: *NEURAL NETWORKS--* Computer Simulation; VOLTAGE REGULATORS 

Identifiers: NEURAL PRISMS; NEGATIVE FEEDBACKS 
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Title: ALGORITHM FOR APPROXIMATE SOLUTIONS TO WEAKLY FILTERED SYNCHRONOUS 

CONTROL SYSTEMS AND NONLINEAR RENEWAL PROCESSES. 
Author: Hoppensteadt , F. C. 

Corporate Source: Univ of Utah, Dep of Mathematics, Salt Lake City, Utah, 
USA 

Source: SIAM Journal on Applied Mathematics v 43 n 4 Aug 1983 p 834-843 

Publication Year: 1983 

CODEN: SMJMAP ISSN: 0036-1399 

Language : . ENGLISH ^ t . 

Journal Announcement: 8312 

Abstract: A multi-time perturbation algorithm is derived to study certain 
synchronous control systems and nonlinear renewal processes. The new 
methods presented here study singular perturbation problems for nonlinear 
Volterra integro-dif f erential equations whose kernels are close to the 
Dirac delta function. The algorithm provides useful approximate solutions 
to weakly filtered phase-locked loop circuits and to nonlinear renewal 
equations describing certain population dynamics. 11 refs. 

Descriptors: ^CONTROL SYSTEMS; MATHEMATICAL TECHNIQUES-- 
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Title: INTEGRATE -AND -FIRE MODELS OF NERVE MEMBRANE RESPONSE TO 
OSCILLATORY INPUT. 

Author: Keener, J. P.; Hoppensteadt, F. C. ; Rinzel, J. 

Corporate Source: Univ of Utah, Salt Lake City, USA 

Source: SIAM Journal on Applied Mathematics v 41 n 3 Dec 1981 p 503-517 

Publication Year: 1981 

CODEN: SMJMAP ISSN: 0036-1399 

Language: ENGLISH 

Journal Announcement: 8205 

Abstract: Nerve membranes exhibit curious responses to alternating 
current stimulation, among which are phase locking, as well as reponses 
without apparent periodic pattern. These phenomena are investigated by 
presenting a complete analysis of the response to periodic input of an 
integrate-and-f ire model, which is a simplified version of the 
Hodgkin-Huxley theory for space clamped nerves. 13 refs. 
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Author: Hoppensteadt , F. C. ; Murray, J. D. 
Corporate Source: Univ of Utah, Salt Lake City 
Source: Mathematical Biosciences v 53 n 1-2 Feb 1981 p 79-87 
Publication Year: 1981 
CODEN: MABIAR ISSN: 0025-5564 
Language: ENGLISH 
Journal Announcement: 8107 

Abstract: A model of an individual's response to a drug is formulated 
first. This is based on interaction between applied dosage of the drug and 
active and inactive receptor sites in the body. Next, a population of such 
individuals is considered. These are divided among nonusers ( susceptibles ) , 
active users, and individuals removed through treatment and cure, 
disenchantment, etc. Recruitment of susceptibilities to active drug use is 
assumed to occur through contact with active users, and the effectiveness 
of recruitment (or infectiousness) depends on the age of a user. The model 
is based solely on the user's response to the drug, and it is shown that 
when a certain combination of susceptible population size, individual 
susceptibility, and infectiousness does not exceed a critical threshold 
value, there will be only few users. But when the threshold value is 
exceeded, an" epidemic *of drug use ensues 2 refsl ~ * * • 
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Author: Levi, M. ; Hoppensteadt , F. C. ; Miranker, W. L. 
Corporate Source: New York Univ, Courant Inst, NY 

Source: Quarterly of Applied Mathematics v 36 n 2 Jul 1978 p 167-198 

Publication Year: 1978 

CODEN: QAMAAY ISSN: 0033-569X 

Language: ENGLISH 

Journal Announcement: 7902' " * ^ " .... 

Abstract: A study is carried out of the sine-Gordon equation and systems 
of discrete approximations to it which are respectively a model of the 
Josephson junction and models of coupled-point Josephson junctions. The 
so-called current-drive case is considered. The voltage response of these 
devices is the average of the time derivative of the solution of these 
equations and the authors demonstrate the existence of running periodic 
solutions for which the average exists. Static solutions are also studied. 
These together with the running solutions explain the multiple-valuedness 
in the response of a Josephson junction in several cases. The stability of 
the various solutions is described in some of the cases. Numerical results 
are displayed which give details of structure of solution types in the case 
of a single point junction and of the one-dimensional distributed junction. 
10 refs. 
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Abstract: It is shown here that certain systems of nonlinear (parabolic) 
reaction-diffusion equations have solutions which are approximated by 
oscillatory functions in the form R( XI MINUS c TAU )P(t*) where P(t) 
represents a sinusoidal oscillation on a fast time scale t and R( XI MINUS 
c TAU ) represents a slowly-varying modulating amplitude on slow space ( XI 
) and slow time ( TAU ) scales. Such solutions describe phenomena in 
chemical reactors, chemical and biological reactions, and in other media 
where a stable oscillation at each point (or site) undergoes a slow 
amplitude change due to diffusion. 19 refs. 
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Abstract: This is a study of the large time behavior of solutions x(t) to 
the nonlinear equation x(t) = GAMMA INTEGRAL x( ETA ) (1 MINUS x( ETA ))d 
ETA , where the integration is from ETA = t MINUS 1 to ETA = 1, which 
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known that such solutions approach constants as t YIELDS INFINITY . This 
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Abstract: Several examples are presented which illustrate some 
capabilities and some limitations of the method of matched asymptotic 
expansions for solving evolution equations. The results are listed 
according to^ spectral properties of the linear problem ^resulting . near a, 
known steady state of the system. When the linear problem is stable, it is 
shown that the solution can be written as a (finite) sum of terms, each 
responding on a different time scale. When the linear problem is unstable, 
it is shown that the method can be used to determine initial data which 
excite only decaying modes, and, in the case of bifurcation of new steady 
states, to construct the new states as well as the transients to them. 35 
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Abstract: The model presented here describes the spread of an infection, 
in a population by keeping track of the chronological ages of the 
participants as well as their " class ages 11 (i. e. the length of time 
since entering their present state) . The reasoning behind this model is 
similar to that used in the equation of age dependent population growth. 7 
ref s . 
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Abstract: Important components of neural networks are input synapses, 

action potential generators and output synapses. Rather than modeling a 
whole neuron in terms of a few ionic channels or as having 
Hodgkin-Huxley, Morris-Lecar or FitzHugh-Nagumo dynamics, we describe a 
neuron's action potential generator (APG) . An APG may be at the hillock 
region at the base of an axon or another specific region of a cell. We 
model it using bifurcation theory based on observations by A. F. 
Hodgkin about membrane excitability. The result is a simplified model 
that leads us to view a neural network as comprising input and output 
synapses (electrical or chemical) that network APGs . These centers of 



activity are coupled by transfer functions from input synapses to an 
APG and from an APG to output synapses. The transfer functions account 
for time delays and signal attenuation that result from input and 
output structures. While this falls far short of a complete biophysical 
model of specific neurons in a network, it is consistent with empirical 
data, it is easily formulated, it is analytically tractable, and 
computer simulations' based on it' are straightforward. One outcome is a 
precise description of the cumulative distribution function (CDF) of 
action potentials. Since records of cell firing amount to collections 
of CDFs, the model is for a variable that is accessible to experimental 
observation. This, methodology is applied here to describe bursting 
neural circuits and embedded loop networks similar to those occurring 
in basal ganglia. 
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Abstract: A neuronal network inspired by the anatomy of the cerebral cortex 
was simulated to study the self-organization of spiking neurons into 
neuronal groups. The network consisted of 100 000 reentrantly 
interconnected neurons exhibiting known types of cortical firing 
patterns, receptor kinetics, short-term plasticity and long-term 
spike-timing-dependent plasticity (STDP), as well as a distribution of 
axonal conduction delays. The dynamics of the network allowed^us to 
study the fine temporal structure o*f emerging "firing patterns with 



millisecond resolution. We found that the interplay between STDP and 
conduction delays gave rise to the spontaneous formation of neuronal 
groups - sets of strongly connected neurons capable of firing 
time-locked, although not necessarily synchronous, spikes. Despite the 
noise present in the model, such groups repeatedly generated patterns 
• of activity 'wrt'h millisecond s^ike-^iming precision ' Exploration of'the 
model allowed us to characterize various group properties, including 
spatial distribution, size, growth, rate of birth, lifespan, and 
persistence in the presence of synaptic turnover. Localized coherent 
input resulted in shifts of receptive and projective fields in the 
model similar to those observed in vivo. 
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Abstract: Networks of coupled periodic oscillators (similar to the Kuramoto 
model) have been proposed as models of associative memory. However, 
error-free retrieval states of such oscillatory networks are typically 
unstable, resulting in a near zero capacity. This puts the networks at 
disadvantage as compared with the classical Hopfield network. Here we 
propose a simple remedy for this undesirable property and show 
rigorously that the error-free capacity of our oscillatory, 
associative-memory networks can be made as high as that of the Hopfield 
network. They can thus not only provide insights into the origin of 
biological memory, but can also be potentially useful for applications 
in information science and engineering. 
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Abstract: A model is presented that reproduces spiking and bursting 

behavior of known types of cortical neurons. The model combines the 

biologically plausibility of Hodgkin-Huxley-type dynamics and the 
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model, one can simulate tens of thousands of spiking cortical neurons 

in real time (1 ms resolution) using a desktop PC. 
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Abstract: In this paper we extend the results of Frankel and Kiemel [SIAM 
J. Appl. Math, 53 ( 1993), pp. 1436 - 1446] to a network of slowly 
coupled oscillators. First, we use Malkin ! s theorem to derive a 
canonical phase model that describes synchronization properties of a 
slowly coupled network. Then, we illustrate the result using slowly 
coupledoscillators ( 1) near Andronov - Hopf bifurcations, ( 2) near 
saddle-node on invariant circle bifurcations, and ( 3) near relaxation 
oscillations. We compare and contrast synchronization properties of 
slowly and weakly coupledoscillators. 
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Abstract: Small-world and scale-free networks are known to be more easily 
synchronized than regular lattices, which is usually attributed to the 
smaller network distance between oscillators. Surprisingly, we find 
that networks with a homogeneous distribution of connectivity are more 
synchronizable than heterogeneous ones, even though the average network 
distance is larger. We present numerical computations and analytical 
estimates on synchronizability of the network in terms of its 
heterogeneity parameters. Our results suggest that some degree of . 
homogeneity is expected in naturally evolved structures, such as neural 
networks, where synchronizability is desirable. 

Identifiers — KeyWord Plus(R): COMPLEX NETWORKS; FRAGILITY; RESONANCE; 
EVOLUTION; STABILITY; DYNAMICS; INTERNET; TOPOLOGY; WEB 

Cited References: 

ALBERT R, 2002, V74, P47, REV MOD PHYS 
ALBERT R, 1999, V401, P130, NATURE 

AMARAL LAN, 2000, V97, P11149, P NATL ACAD SCI USA 

BARABASI AL, 1999, V286, P509,^ SCIENCE g 

' BARAHONA M, '2002, V89, 054101," PHYS REV LETT*' 1 " 

BOLLOBAS B, 2002, HDB GRAPHS NETWORKS 

BRAITENBERG V, 1998, P93, CORTEX STAT GEOMETRY 

COHEN R, 2003, V90, 058701, PHYS REV LETT 

DOROGOVTSEV SN, 2000, V50, PI, EUROPHYS LETT 

DOROGOVTSEV SN, 2002, V51, P1079, ADV PHYS 

DOROGOVTSEV SN, 2000, V62, P1842, PHYS REV E A 

FALOUTSOS M, 1999, V29, P251, COMP COMM R 

FINK KS, 2000, V61, P5080, PHYS REV E A' 

GADE PM, 2000, V62, P6409, PHYS REV E A 

GOH KI, 2001, V87, 278701, PHYS REV LETT 

HONG H, 2002, V65, 026139, PHYS REV E 2 

IZHIKEVICH EM, 2003, V26, P161, TRENDS NEUROSCI 

JOST J, 2001, V65, PHYS REV E 

KANDEL ER, 2000, PRINCIPLES NEURAL SC 

KARBOWSKI J, 2001, V86, P3674, PHYS REV LETT 

KLEMM K, 2002, V65, 036123, PHYS REV E 2A 

KLEMM K, 2002, V65, 057102, PHYS REV E 2 

KWON 0, 2002, V298, P319, PHYS LETT A 

LAGOFERNANDEZ LF, 2000, V84, P2758^ PHYS REV^LETT, 
' LATORA' V, 2001,' V87,' 198701, PHYS REV ' LETT 

LILJEROS F, 2001, V411, P907, NATURE 

MONTOYA JM, 2002, V214, P4 05, J THEOR BIOL 

MOTTER AE, 2002, V65, 065102, PHYS REV E 2 

MOTTER AE, 2002, V66, 065103, PHYS REV E 2 



NEWMAN ME J, 2000, V84, P3201, PHYS REV LETT 
NEWMAN ME J, '2001, V64, 026118/ PHYS* REV E 2 
NEWMAN MEJ, 2001, V98, P404, P NATL ACAD SCI USA 
NEWMAN MEJ, 2001, V64, 016132, PHYS REV E 2 
NISHIKAWA T, 2002, V66, 046139, PHYS REV E 2 
PECORA LM, 1998, V80, P2109, PHYS REV LETT 
SOLE RV, 2001, V268, P2039, P ROY SOC LOND B BIO 
STEPHAN KE, 2000, V355, Pill, PHILOS T ROY SOC B 
STROGATZ SH, 2001, V410, P268, NATURE 
WANG XF, 2002, V49, P54, IEEE T CIRCUITS-I 
WANG XF, 2002, V12, P885, INT J BIFURCAT CHAOS 
WANG XF, 2002, V12, P87, INT J BIFURCAT CHAOS 
WATTS DJ, 1999, SMALL WORLDS 
WATTS DJ, 1998, V393, P440, NATURE 



1/5/87 (Item 10 from file: 34) 

DIALOG ( R) File 34 : SciSearch ( R) Cited Ref Sci 
(c) 2004 Inst for Sci Info. All rts. reserv. 

11727667 Genuine Article#: 688EW Number of References: 19 
Title: Relating' STDP to BCM * /t . . . „ . 

Author (s): Izhikevich EM (REPRINT) ; Desai NS 

Corporate Source: Inst Neurosci,San Diego//CA/92121 (REPRINT); Inst 

Neurosci,San Diego//CA/ 92121 
Journal: NEURAL COMPUTATION, 2003, V15, N7 (JUL) , P1511-1523 
ISSN: 0899-7667 Publication date: 20030700 

Publisher: MIT PRESS, FIVE CAMBRIDGE CENTER, CAMBRIDGE, MA 02142 USA 
Language: English Document Type: ARTICLE 
Geographic Location: USA 

Journal Subject Category: COMPUTER SCIENCE, ARTIFICIAL INTELLIGENCE 
Abstract: We demonstrate that the BCM learning rule follows directly from 

STDP when pre- and postsynaptic neurons fire uncorrelated or weakly 

correlated Poisson spike trains, and only nearest-neighbor spike 

interactions are taken into account. 
Identifiers — KeyWord Plus (R) : TERM SYNAPTIC PLASTICITY; VISUAL- CORTEX; 

PYRAMIDAL CELLS; NMDA RECEPTORS; DEPENDENCE; EXPERIENCE; NEURONS; 

MODEL; LTP 
Cited References: 

ABARBANEL HDI, 2002, V99, P10132, P NATL ACAD SCI USA 

BEAR MF, 1994, V4, P389, CURR OPIN NEUROBIOL 

BEAR MF,- 1987/ V23-7, P42/ SCIENCE ' ■ . ■ , . 

BI GQ, 1998, V18, P10464, J NEUROSCI 

BIENENSTOCK EL, 1982, V2 , P32, J NEUROSCI 

CASTELLANI GC, 2001, V98, P12772, P NATL ACAD SCI USA 

DEBANNE D, 1998, V507, P237, J PHYSIOL-LONDON 

FELDMAN DE, 2000, V27, P45, NEURON 

FROEMKE RC, 2002, V416, P433, NATURE 

KEMPTER R, 1999, V59, P4498, PHYS REV E 

KEMPTER R, 2001, VI 3, P27 09, NEURAL COMPUT 

KIRKWOOD A, 1996, V381, P526, NATURE 

KIRKWOOD A, 1993, V260, P1518, SCIENCE 

MARKRAM H, 1997, V275, P213, SCIENCE 

PHILPOT BD, 2001, V29, P157, NEURON 

SENN W, 2001, V13, P35, NEURAL COMPUT 

SJOSTROM PJ, 2001, V32, P1149, NEURON 

SONG S, 2000, V3, P919, NAT NEUROSCI 

VANROSSUM MCW, 2000, V20, P8812, J NEUROSCI 



1/5/88 (Item 11 from file: 34) 

DIALOG ( R) Fil'e 34 rSciffearch ( R) Cited' RSf Sci .. .. 

(c) 2004 Inst for Sci Info. All rts. reserv. 

11439623 Genuine Article#: 651YQ Number of References: 34 
Title: Bursts as a unit of neural information:, selective communication via 
resonance 



Author (s ):■ -izhikevich EM (REPRINT) « ; **Desai NS; 'Walcotrt EC; -Hoppensteadt 
FC 

Corporate Source: Inst Neurosci, 1064 0 John Jay Hopkins Dr/San 

Diego//CA/92121 (REPRINT); Inst Neurosci, San Diego//CA/92121 ; Arizona 
State Univ,Ctr Syst Sci, Tempe//AZ/85287 

Journal: TRENDS IN NEUROSCIENCES, 2003, V26, N3 (MAR), P161-167 

ISSN: 0166-2236 Publication date: 20030300 

Publisher: ELSEVIER SCIENCE LONDON, 84 THEOBALDS RD, LONDON WC1X 8RR, 
ENGLAND 

Language: English Document Type: REVIEW 
Geographic Location: USA 

Journal Subject Category: NEUROSCIENCES 

Abstract: What is the functional significance of generating a burst of 

spikes, as opposed to a single spike? A dominant point of view is that 
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neurons. Here, we discuss the alternative, but complementary, 
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Abstract: Motivated by the practical consideration of the measurement of 
chaotic signals in experiments or the transmission of these signals 
through a physical medium, we investigate the effect of filtering on 
chaotic symbolic dynamics. We focus on the linear, time-invariant 
filters that are used frequently in many applications, and on the two 
quantities characterizing chaotic symbolic dynamics: topological 
entropy and bit-error rate. Theoretical consideration suggests that the 
topological entropy is invariant under filtering. Since computation of 
this entropy requires that the generating partition for defining the 
symbolic dynamics be known, in practical situations the computed 
entropy may change as a filtering parameter is changed. We find, 
through numerical computations and experiments with a chaotic 
electronic circuit, that with reasonable care the computed or measured 
entropy values can be preserved for a wide range of the filtering 
parameter. (C) 2003 American Institute of Physics. 
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paper we present an alternative but complementary hypothesis. We 
consider the effect of a short burst on a model postsynaptic cell 
having damped oscillation of its membrane potential. The oscillation 
frequency (eigenf requency ) plays a crucial role. Due to the 
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variable is complex. The model provides geometric illustrations to many 
interesting phenomena occurring in biological neurons having 
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such neurons prefer a certain resonant frequency of the input that is 
nearly equal to their eigenf requency, they can be excited or inhibited 
by a doublet (two pulses) depending on its interspike interval, and 
they can fire in response to an inhibitory input. All these properties 
could be observed in Hodgkin-Huxley-type models. We use the 
resonate-and-f ire model to illustrate possible sensitivity of 
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train. Being an analogue of the integrate-and-f ire model, the 
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Abstract: Periodic bursting behavior in neuron is a recurrent transition 
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bifurcation and the transition to the quiescent state occurs via fold 
limit cycle bifurcation, the burster is said to be of elliptic type 
(also know as a "subHopf / f old cycle" burster) . Here we study the 
synchronization dynamics of weakly connected networks of such bursters. 
We find that the behavior of such networks is quite different from the 
behavior of weakly connected phase oscillators and resembles that of 
strongly connected relaxation oscillators. As a result, such weakly 
connected bursters need few (usually one) bursts to synchronize, and 
synchronization is possible for bursters having quite different 
quantitative features. We also find that interactions between bursters 
depend crucially on the spiking frequencies. Namely, the interaction 
are most effective when the presynaptic interspike frequency matches 
the frequency of postsynaptic oscillations. Finally, we use the 
FitzHugh Rinzel, Morris-Lecar , and Hodgkin-Huxley models to illustrate 
our major results. 
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Abstract: The diffusion gradient chamber (DGC) is a novel device developed 
to study the response of chemotactic bacteria to combinations of 
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given. An ADI (Alternating Direction Implicit) method is constructed 
for finding numerical solutions of the model and carrying out computer 
simulations. The numerical results of the model successfully reproduced 
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continuous (possibly noninvertible) change of variables. We also derive 
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Abstract: Many sci.enti.sts believe ..that, *all pulseycpupled . neural networks 
are toy models that are far away from the biological reality, We show 
here, however, that a huge class of biophysically detailed and 
biologically plausible neural-network models can be transformed into a 
canonical pulse-coupled form by a piece-wise continuous, possibly 
noninvertible, change of variables, Such transformations exist when a 
network satisfies a number of conditions; e.g., it is weakly connected; 
the neurons are Class 1 excitable (i.e., they can generate action 
potentials with an arbitrary small frequency) ; and the synapses between 
neurons are conventional (i.e., axo-dendritic and axe-somatic). Thus, 
the difference between studying the pulse-coupled model and 
Hodgkin-Huxley-type neural networks is just a matter of a coordinate 
change. Therefore, any piece of information about the pulse-coupled 
model is valuable since it tells something about all weakly connected 
networks of Class 1 neurons. For example, we show that the 
pulse-coupled network of identical neurons does not synchronize 
in-phase. This confirms Ermentrout ! s result that weakly connected Class 
1 neurons are difficult to synchronize, regardless of the equations 
that describe dynamics of each cell. 
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Abstract: Our study of thalamo-cortical systems suggests a new architecture 
for a neurocomputer that consists of oscillators having different 
frequencies and that are connected weakly via a common medium forced by 



an external input. Even though such oscillators are all interconnected 

homogeneously, the external input imposes a dynamic connectivity. We 

use Kuramoto's model to illustrate the idea and to prove that such a 

neurocomputer has oscillatory associative properties. Then we discuss a 

general case. The advantage of such a neurocomputer is that it can be 

built using voltage controlled oscillators, optical oscillators, 

lasers, microelectromechanical systems, Josephson junctions, 

macromolecules, or oscillators of other kinds. (Provisional patent 

60/108,353) [S0031-9007 (99) 08813-4] . 
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Abstract: We consider all models of the thalamo-cortical system that 

satisfy the following two assumptions: (I) each cortical column is an 
autonomous oscillator; (2) connections between cortical columns and the 
thalamus" are weak.' Our goal is to deduce from these assumptions general 
principles of thalamo-cortical interactions that are independent of the 
equations describing the system. We find that the existence of synaptic 
. connections between any two cortical columns does not guarantee that 
the columns interact: They interact only when there is a certain nearly 
resonant relation between their frequencies, which implies that the 
interactions are frequency modulated (FM) . When the resonance relation 
holds, the cortical columns interact through phase modulations. Thus, 
communications between weakly connected cortical oscillators employ a 
principle similar to that in FM radio: The frequency of oscillation 
encodes the channel of communication, while the information is 
transmitted via phase modulations. If the thalamic input has an 
appropriate frequency, then it can dynamically link any two cortical 
columns, even those that have non-resonant frequencies and would 
otherwise be unlinked. Thus, by adjusting its temporal activity, the 
thalamus has control over information processing taking place in the 
cortex. Our results suggest that the mean firing rate (frequency) of 



periodically spiking neuron does not carry any information other than 

identifying a channel of communication. Information (i.e. neural code) 

is carried through modulations of interspike intervals. (C) 1998 

Elsevier Science Ireland Ltd. All rights reserved. 
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Abstract: We describe the model of the hippocampus consisting of 
interactive oscillators with input from the entorhinal cortex 
(modulating the main information dow by a theta rhythm) and the septum 
(a theta rhythm generator) . When interconnections between oscillators 
are allowed to strengthen in an adaptive way, the network can be 
trained using a series of lessons. This results in a connection matrix 
that memorizes the temporal sequence of inputs. Presenting one of the 
lessons to the trained network results in reproduction of the remainder 
of the sequence. In this paper, we create such a connection matrix, 
derive from it an appropriate Markov chain and simulate the chain to 
illustrate its dynamics. (C) 1998 Elsevier Science Ireland Ltd. All 
rights reserved. 
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Abstract: Studying weakly connected oscillators leads to phase models. It 

has been proven recently that weakly connected oscillators with delayed 

interactions do not lead to phase models with time delays even when the 

delay is of the same order of magnitude as the period of oscillation. 

This has resulted in a fading interest in such models. We prove here 

that if the interaction delay between weakly connected oscillators is 

much longer than the period of oscillation, then the corresponding 

phase model does have an explicit time delay. 
Identifiers--KeyWord Plus (R) : OSCILLATORS 
Cited References: 

BLECHMAN II, 1971, SYNCRHONIZATION DYAN 

BUCK J, 1976, V234, P74, SCI AM 
. ERMENTROUT GB,. .198.1,. V12 r - P327 v J MATH BIOL 

ERMENTROUT GB, 1991, V29, P195, J MATH BIOL 

ERMENTROUT GB, 1994, P79, NEURAL MODELING NEUR 

FENICHEL N, 1971, V21, P193, INDIANA U MATH J 

GRAY CM, 1994, VI, Pll, J COMPUT NEUROSCI 

HIRSCH MW, 1977, INVARIANT MANIFOLDS 

HOPPENSTEADT F, 1997, WEAKLY CONNECTED NEU 

KIM S, 1997, V79, P2911, PHYS REV LETT 

KURAMOTO Y, 1984, CHEM OSCILLATIONS WA 

LUZYANINA TB, 1995, V6, P43, NETWORK-COMP NEURAL 

MALKIN IG, 1949, METHODS POINCARE LIA 

MALKIN IG, 1956, SOME PROBLEMS NONLIN 

NIEBUR E, 1991, V67, P2753, PHYS REV LETT 

PESKIN CS, 1975, MATH ASPECTS HEART P 

SHEPHERD GM, 1983, NEUROBIOLOGY 



1/5/113 (Item 36 from file: 34) 

DIALOG (R) File 34 : SciSearch { R) Cited Ref Sci 



(c) 2004 Inst for Sci Info. All rts . reserv. 

06881363 Genuine Article#: ZZ053 Number of References: 30 
Title: Multiple cusp bifurcations 
Author (s): Izhikevich EM (REPRINT) 

Corporate Source: ARIZONA STATE UNIV, CTR SYST SCI & ENGN/TEMPE//AZ/85287 

(REPRINT); MICHIGAN STATE UNIV, DEPT MATH/E LANSING/ /MI/ 4 882 4 
Journal: NEURAL t NETWORKS, 1998, Vl^N^ (APR), Pf9 t 5-50.8_ . . , 

ISSN: 0893-6080' Publication date: 19980400 

Publisher: PERGAMON- ELSEVIER SCIENCE LTD, THE BOULEVARD, LANGFORD LANE, 

KIDLINGTON, OXFORD OX5 1GB, ENGLAND 
Language: English Document Type: ARTICLE 
Geographic Location: USA 

Subfile: CC ENGI — Current Contents, Engineering, Computing & Technology 
Journal Subject Category: COMPUTER SCIENCE, ARTIFICIAL INTELLIGENCE 
Abstract: The cusp bifurcation provides one of the simplest routes leading 
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dynamics. First, we prove that any such network can be transformed into 
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Abstract: This is the second of two articles devoted to analyzing the 
relationship between synaptic organizations (anatomy) and dynamical 
properties (function) of networks of neural oscillators near multiple 
supercritical Andronov-Hopf bifurcation points. Here we analyze 
learning processes in such networks. Regarding learning dynamics, we 
assume (1) learning is local (i.e. synaptic modification depends on 
pre- and postsynaptic neurons but not on others), (2) synapses modify 
slowly relative to characteristic neuron response times, (3) in the 
absence of either pre- or postsynaptic activity, the synapse weakens 
(forgets) . Our major goal is to analyze all synaptic organizations of 
oscillatory neural networks that can memorize and retrieve phase 
information or time delays. We show that such networks have the 
following attributes: (1) the rate of synaptic plasticity connected 
with learning is determined locally by the presynaptic neurons, (2) the 
excitatory neurons must be long-axon relay neurons capable of ^forming 
distant "connections with other excitatory anci inhibitory neurons, (3) 
if inhibitory neurons have long axons, then the network can learn, 



passively forget and actively unlearn information by adjusting synaptic 
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Abstract: We study weakly connected networks of neural oscillators near 
multiple Andronov-Hopf bifurcation points. We analyze relationships 
between synaptic organizations (anatomy) of the networks arid their 
dynamical properties (function) . Our principal assumptions are: (1) 
Each neural oscillator comprises two populations of neurons: excitatory 
and inhibitory ones; (2) activity of each population of neurons is 
described by a scalar (one-dimensional) variable; (3) each neural 
oscillator is near a nondegenerate supercritical Andronov-Hopf 
bifurcation point; (4) the synaptic connections between the neural 
oscillators are weak. 
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All neural networks satisfying these hypotheses are governed by the 
same dynamical system, which we call the canonical model. Studying the 
canonical model shows that: (1) A neural oscillator can communicate 
only with those oscillators which have roughly the same natural 
frequency. That is, synaptic connections between a pair of oscillators 
having different natural frequencies are functionally insignificant. 
(2) Two neural oscillators having the same natural frequencies might 
not communicate if the connections between them are from among a class 
of pathological synaptic configurations. In both cases the anatomical 
presence of synaptic connections between neural oscillators does not 
necessarily guarantee that the connections are functionally 
significant. (3) There can be substantial phase differences (time 
delays) between the neural oscillators, which result from the synaptic 
organization of the network, not from the transmission delays. Using 
the canonical model we can illustrate self-ignition and autonomous 
quiescence (oscillator death) phenomena. That is, a network of passive 
elements can exhibit active properties and vice versa. We also study 
how Dale's principle affects dynamics of the networks, in particular, 
the phase differences that the network can reproduce. We present a 
complete, class.if" iq.at.ion of alL possj-ble synaptic, organizations from, 
this point of view. The theory developed here casts some light on 
relations between synaptic organization acid functional properties of 
oscillatory networks, The major advantage of our approach is that we 



obtain results about all networks of neural oscillators, including the 

real brain. The major drawback is that our findings^are valid^only when 

the brain operates' hear a critical "regime, viz*, for a multiple 

Andronov-Hopf bifurcation. 
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01330798 INSIDE CONFERENCE ITEM ID: CN013173034 

Canonical Models for Bifurcations from Equilibrium in Weakly Connected 
Neural Networks 

Hoppensteadt, F. C. ; Izhikevich, E. 

CONFERENCE: WCNN ' 95-World congress on neural networks 

WORLD CONGRESS ON NEURAL NETWORKS , 1995; VOL 1 P: I-80-I-83 

Lawrence Erlbaum Associates, 1995 

ISBN: 0805821252 

LANGUAGE: English DOCUMENT TYPE: Conference Papers 

CONFERENCE SPONSOR: International Neural Network Society 
CONFERENCE LOCATION: Washington, DC 
CONFERENCE DATE: Jul 1995 (19950) (19950) 

BRITISH LIBRARY ITEM LOCATION :. 9353 . 44 64 80 
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Held as the annual meeting of the INNS 
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01184449 INSIDE CONFERENCE ITEM ID: CN011619198 

Persistence of Singular Perturbation Solutions in Noisy Environments 
Hoppensteadt, F. C. 

CONFERENCE: Trends and developments in ordinary differential equations - - 
International symposium on ordinary differential equations and 
applications 
P: 141-142 

World Scientific, 1994 

ISBN: 9810215304 

LANGUAGE: English DOCUMENT TYPE: Conference Papers 
CONFERENCE EDITOR(S): Alavi, Y. ; Hsieh, P.-F. 

CONFERENCE LOCATION: Kalamazoo,, MI_ . . ... 

CONFERENCE DATE: May 1993 (199305) (199305) 

BRITISH LIBRARY ITEM LOCATION: 95/31713 Trends 
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Simulation of Tonotopic Neural Circuit 
Hoppensteadt, F. C. 

CONFERENCE: Mathematics applied to biology and medicine-lst European 
congress 
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Pattern recognition via synchronization in phase -locked loop neural 
networks 

Hoppensteadt, FC ; Izhikevich, EM 
Center for Syst. Sci. & Eng., Arizona State Univ., Tempe, AZ, USA 
IEEE Transactions on Neural Networks, vll, n3, pp734-738, 2000 
Document type: journal article Language: English 
Record type: Abstract 
ISSN: 1045-9227 

ABSTRACT : . 

We propose a novel architecture of an oscillatory neural network that 
consists of phase-locked loop (PLL) circuits. It stores and retrieves 
complex oscillatory patterns as synchronized states with appropriate phase 
relations between r\eurpns . ... . . . , 

DESCRIPTORS: NEURAL CHIPS; IMAGE RECOGNITION; PHASE LOCKED LOOPS; 
SYNCHRONIZATION; VCO- -VOLTAGE CONTROLLED OSCILLATORS 
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Oscillatory model of the hippocampal memory 
Borisyuk, R; Hoppensteadt, F 
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IJCNN'99. International Joint Conference on Neural Networks. Proceedings 
(Cat. No.99CH36339) , 10-16 July 1999, Washington, DC, USA1999 
Document type: Conference paper Language: English 
Record type: Abstract 
ISBN: 0-7803-5529-6 

ABSTRACT: 

We describe a biologically inspired oscillatory neural network for 
memorizing temporal sequences of neural activity patterns. The neural 
network consists of interactive neural oscillators with all-to-all 
excitatory connections forced by a slow T-periodic signal. The dynamics of 
the network are viewed through a time window with duration T. The network 
memorizes binary patterns in terms of low and high activity of the 
corresponding oscillators. The learning rule is temporally asymmetric, and 
it takes into account the activity level of pre- and post- ' synaptic 1 
oscillators in two contiguous time windows. Recall by the network is fast: 
all memorized patterns of sequences are reproduced in the correct order 
during the same time window, but with a short time delay. The applicability 
of these results to studies of the hippocampus is discussed. 

DESCRIPTORS: ARTIFICIAL NEURAL NETWORKS; NEUROPHYSIOLOGY; DELAY TIME; 
SYSTEM THEORY 
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Slowly Coupled Oscillators: Phase Dynamics and Synchronization 

Izhikevich, Eugene M; Hoppensteadt , Frank C 

SIAM Journal on Applied Mathematics v. 63 no6 (Aug. /Sept. 2003) p. 1935-53 
DOCUMENT TYPE: Feature Article ISSN: 0036-1399 LANGUAGE: English 
RECORD STATUS: New record 

ABSTRACT: In this paper we extend the results of Frankel and Kiemel [SIAM 
J. Appl. Math, 53 (1993), pp. 1436-1446] to a network of slowly coupled 
oscillators. First, we use Malkin ! s theorem to derive a canonical phase 
model that describes synchronization properties of a slowly coupled 
network. Then, we illustrate the result using slowly coupled oscillators 
(1) near Andronov-Hopf bifurcations, (2) near saddle-node on invariant 
circle bifurcations, and (3) near relaxation oscillations. We compare and 
contrast synchronization properties of, js lowly ancl ^weakly . coupled , „ . 
oscillators. Reprinted by permission of the publisher. 
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System of Phase Oscillators with Diagonalizable Interaction 

Nishikawa, Takashi; Hoppensteadt, Frank C 

SIAM Journal on Applied Mathematics v. 63 no5 (June/Aug. 2003) p. 1615-26 
DOCUMENT TYPE: Feature Article ISSN: 0036-1399 LANGUAGE: English 
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ABSTRACT: We consider a system of N phase oscillators having randomly 
distributed natural frequencies and diagonalizable interactions among the 
oscillators. We show that, in the limit of N --> [infinity], all solutions 
of such a system are incoherent wi£h # probability one. far . any .strength of 
coupling, which implies that there is no sharp transition from incoherence 
to coherence as the coupling strength is increased, in striking contrast to 
Kuramoto 1 s (special) oscillator system. Reprinted by permission of the 
publisher . 
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Numerical and experimental investigation of the effect of filtering on 

chaotic symbolic dynamics 
Zhu, Liqiang; Lai, Ying-Cheng; Hoppensteadt, Frank C 
Chaos v. 13 nol (Mar. 2003) p. 410-19 

DOCUMENT TYPE: Feature Article ISSN: 1054-1500 LANGUAGE: English 
RECORD STATUS: Corrected or revised record 

ABSTRACT: Part of a special issue on the control and synchronization of 
chaotic dynamical systems. An investigation of the effect of filtering on 
chaotic symbolic dynamics is presented. Special attention is given to the 
linear, time-invariant filters that are frequently employed in many 
applications and to the 2 quantities representing chaotic symbolic 
dynamics: topological entropy and bit-error rate. Theoretical studies 
indicate that the topological entropy does not vary under filtering. As 
computation of this entropy requires that the generating partition for 
defining the symbolic dynamics be known, in practical situations the 
computed entropy may vary as a filtering parameter is varied. It was 



found, via numerical computations and experiments with a chaotic electronic 
circuit, that with reasonable care the computed or measured entropy values 
can be preserved for a wide range of the filtering parameter. 
DESCRIPTORS: Chaos (Science); Linear time varying systems; Bit error rates 



1/5/148 (Item 4 from file: 99) 

DIALOG (R) File 99:Wilson Appl . Sci & Tech Abs 
(c) 2004 The HW Wilson Co. All rts. reserv. 

2575332 H.W. WILSON RECORD NUMBER: BAST01024174 

Synchronization of MEMS resonators and mechanical neurocomputing 

Hoppensteadt, Frank C ; Izhikevich, Eugene M 
IEEE Transactions on Circuits and Systems. Part I, Fundamental Theory and 
Application's v. '48* 'noZ (Feb . "2001 T p*. 1*33-8 

DOCUMENT TYPE: Feature Article ISSN: 1057-7122 LANGUAGE: English 
RECORD STATUS: Corrected or revised record 

ABSTRACT: The authors developed a network of coupled 
microelectromechanical system (MEMS) oscillators that function as a 
neurocomputer with an oscillatory autocorrelative associative memory. The 
system is based on the fact that networks of arbitrary oscillators have 
associated memory when coupled appropriately. A canonical model describing 
the nonlinear dynamics of a single MEMS oscillator is presented and a 
network of such oscillators is considered. The theory is applied to a 
pattern recognition problem. 

DESCRIPTORS: MEMS resonators; Neural networks; Associative memories; 
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Pattern recognition via synchronization in phase-locked loop neural 
networks 

Hoppensteadt, Frank C ; Izhikevich, Eugene M 
IEEE Transactions on Neural Networks v. 11 no3 (May 2000) p. 734-8 
DOCUMENT TYPE: Feature Article ISSN: 1045-9227 LANGUAGE: English 
RECORD STATUS: Corrected or revised record 

ABSTRACT: The authors developed a novel architecture for an oscillatory 
neural network that consisted of phase-locked loop circuits. The circuit 
stores and retrieves complex oscillatory patterns as synchronized states 
with appropriate phase relations between neurons. 

DESCRIPTORS: Phase locked loops--Design ; Pattern recognition — Neural 
network models; Synchronism; 
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Weakly connected neural networks [book review] 

Hoppensteadt, F. C (Frank Charles; Izhikevich, Eugene M; Ermentrout, Bard 
reviewer 

SIAM Review v. 41 nol (Mar. '99) p. 178-9 
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Singular perturbation solutions of noisy systems 
Hoppensteadt , Frank C ; 

SIAM Journal on Applied Mathematics v. .55 (Apr. '95) p. 544-51 
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Mathematics in medicine and the life sciences [book review] 

Hoppensteadt, F. C (Frank Charles; Peskin, Charles. S; Milton, John 
reviewer 

SIAM Review v. 36 (Mar. '94) p. 134-5 
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A particle method for population waves 
Chiu, Chichia; Hoppensteadt, Frank C 

SIAM Journal on Applied Mathematics v. 54 (Apr. '94) p. 466-77 
DOCUMENT TYPE: Feature Article ISSN: 0036-1399 LANGUAGE: English 
RECORD STATUS: New .record . . ... 

ABSTRACT: The authors derive a particle technique for solving the weak 
formulation of a phase model that permits the growth rate of cells to 
change at different phases of cell cycle. Phase models are useful for 
analyzing synchronization of bacterial cell culture growth and other 
biological occurrences associated with cell cycles. Particle techniques 
are numerical techniques that are especially useful in computational fluid 
dynamics, where they are sometimes referred to as vortex methods. The 
particle technique can be a natural way of explaining how synchronized cell 
cultures can be induced by starvation-nutrition cycles. The convergence of 
the proposed particle method is proved for the linear version of the model. 
Error estimates for the approximation are derived. 

DESCRIPTORS: Bacteriology — Cultures and culture media; Differential 
equations — Numerical solutions ; Cell division (Biology) — Mathematical 
models; 
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Signal processing by model neural networks 
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SIAM Review v. 34 (Sept. '92) p. 426-44 
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Title : Dynamics of the Josephson junction 

Author(s): Levi, M. ; Hoppensteadt , F.C. ; Miranker, W.L 

Affiliation: Courant Institute, New York University 

Source: Q. Appl . Math. (United States) v 37:3. Coden: QAMAA 

Publication Date: Jul 1978 

p 167-198 

Document Type: Journal Article 

Language: English ^ t tt ... . . . „ 

Journal Announcement: EDB7810 

Subfile: INS (US Atomindex input); AIP (SPIN). 
Country of Origin: United States 

Abstract: We study the sine-Gordon equation and systems of discrete 

approximations to it which are respectively a model of the Josephson 
junction and models of coupled-point Josephson junctions. We do this in 
the so-called current-driven case. The voltage response of these 
devices is the average of the time derivative of the solution of these 
equations and we demonstrate the existence of running periodic 
solutions for which the average exists. Static solutions are also 
studied. These together with the running solutions explain the 
multiple-valuedness in the response of a Josephson junction in several 
cases. The stability of the various solutions is described in some of 
the cases . Numerical results are displayed with give details of 
structure of solution types in the case of a single point junction and 
of the one-dimensional distributed junction.; 
Major Descriptors: * JOSEPHSON JUNCTIONS — SINE-GORDON EQUATION 
Descriptors: BOUNDARY CONDITIONS; COUPLING; ONE- DIMENSIONAL CALCULATIONS 
Broader Terms: EQUATIONS; FIELD EQUATIONS; SUPERCONDUCTING JUNCTIONS 
Subject Categories^ 656101* Solid State Physics. .7- .Superconductivity -- 

' General' Theory ' — '(-1987)' 
INIS Subject Categories: A17* -- Low Temperature Physics 
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Oscillatory associative memory network with perfect retrieval 

NISHIKAWA Takashi; HOPPENSTEADT Frank C ; LAI Ying-Cheng 

Department of Mathematics and Center for Systems Science and Engineering 
Research, Arizona State University, Tempe, Arizona 85287, United States; 
Courant Institute of Mathematical Sciences, New York University, New York 
NY 10012, United States; Department of Electrical Engineering, Arizona 
State University, Tempe, Arizona 85287, United States 

Journal: Physica. D, 2004, 197 (1-2) 134-148 

ISSN: 0167-2789 CODEN: PDNPDT Availability: INIST-145D; 
354000120345920080 

No. of Refs.^22 ref. t ...... 

Document Type: P (Se'rial) ; A (Analytic) 

Country of Publication: Netherlands 

Language: English 

Inspired by the discovery of possible roles of synchronization of 
oscillations in the brain, networks of coupled phase oscillators have been 
proposed before as models of associative memory based on the concept of 
temporal coding of information. Here we show, however, that error-free 
retrieval states of such networks turn out to be typically unstable 
regardless of the network size, in contrast to the classical Hopfield 
model. We propose a remedy for this undesirable property, and provide a 
systematic study of the improved model. In particular, we show that the 
error-free capacity of the network is at least 2 epsilon SUP 2 / log n 
patterns per neuron, where n is the number of oscillators (neurons) and 



epsilon the strength of the second-order mode in the coupling function. 

English Descriptors: Synchronization; Coupled oscillator; Hopfield model; 
Second order; Mode coupling; Neural networks; Random matrix; Non linear 
phenomenon 

French Descriptors: Synchronisation; Oscillateur couple; Modele Hopfield; 
Ordre 2; Couplage mode; Reseau neuronal; Matrice aleatoire; Phenomene non 
lineaire 
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Bursts as a unit of neural information: selective communication via 
resonance 

IZHIKEVICH Eugene M ; DESAI Niraj S; WALCOTT Elisabeth C; HOPPENSTEADT 
Frank C 

The Neurosciences Institute, 10640 John Jay Hopkins Drive, San Diego, CA 
92121, United States; Center for Systems Science, Arizona State University, 
Tempe,- AZ 85287, United States 

Journal: Trends in neurosciences : (Regular edition), 2003, 26 (3) 
161-167 

ISSN: 0166-2236 Availability: INIST-18018B; 354000114 892 680100 
No. of Refs . : 38 ref . . 

Document Type: P (Serial) ; A (Analytic) 

Country o*f Publication: United Kingdom " i- - ■ <• - 

Language: English 

What is the functional significance of generating a burst of spikes, as 
opposed to a single spike? A dominant point of view is that bursts are 
needed to increase the reliability of communication between neurons. Here, 
we discuss the alternative, but complementary, hypothesis: bursts with 
specific resonant interspike frequencies are more likely to cause a 
postsynaptic cell to fire than are bursts with higher or lower frequencies. 
Such a frequency preference might occur at the level of individual synapses 
because of the interplay between short-term synaptic depression and 
facilitation, or at the postsynaptic cell level because of subthreshold 
membrane potential oscillations and resonance. As a result, the same burst 
could resonate for some synapses or cells and not resonate for others, 
depending on their natural resonance frequencies. This observation suggests 
that, in addition to increasing reliability of synaptic transmission, 
bursts of action potentials might provide effective mechanisms for 
selective communication between neurons. 

English Descriptors: Review; Neuron; Discharge pattern; Resonance; 
Communication; Synapse; Synaptic plasticity; Facilitation; Membrane 
potential';" Oscill'atron; Resonance frequency; Synaptic "transmission ; 
Synaptic potential; Action potential 

Broad Descriptors: Electrophysiology ; Electrophysiologie ; Electrof isiologia 

French Descriptors: Article synthese; Neurone; Mode decharge; Resonance; 
Communication; Synapse; Plasticite synaptique; Facilitation; Potentiel 
membrane; Oscillation; Frequence resonance; Transmission synaptique; 
Potentiel synaptique; Potentiel action 
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Capacity of Oscillatory Associative-Memory Networks with Error-Free 
Retrieval 

NISHIKAWA Takashi; LAI Ying-Cheng; HOPPENSTEADT Frank C 

Department of Mathematics, Center for Systems Science and Engineering 
Research, Arizona State University, Tempe, Arizona 85287, USA; Department 
of Electrical Engineering, Arizona State University, Tempe, Arizona 85287, 
USA 

Journal:. Physical review letters*, . 2004-03-12 , < -92 ■■( 10* • 108 101-198101-4 
ISSN: 0031-9007 CODEN: PRLTAO Availability: INIST-8895 
Document Type: P (Serial) ; A (Analytic) 
Country of Publication: United States 
Language: English 

Networks of coupled periodic oscillators (similar to the Kuramoto model) 
have been proposed as models of associative memory. However, error-free 
retrieval states of such oscillatory networks are typically unstable, 
resulting in a near zero capacity. This puts the networks at disadvantage 
as compared with the classical Hopfield network. Here we propose a simple 
remedy for this undesirable property and show rigorously that the 
error-free capacity of our oscillatory, associative-memory networks can be 
made as high as that of the Hopfield network. They can thus not only 
provide insights into the origin of biological memory, but can also be 
potentially useful for applications in information science and engineering. 

English Descriptors: Theoretical study; Neurophysiology; Hopfield neural 
nets; Brain models; Asymptotic stability; Patterning 
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Slowly coupled oscillators: Phase dynamics and synchronization 
IZHIKEVICH E M ; HOPPENSTEADT F C 
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Journal: SIAM Journal on Applied Mathematics, 2003, 63 (6) 1935-1953 
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In this paper we extend the results of Frankel and Kiemel Degree SIAM J. 
Appl. Math, 53 (1993), pp. 1436-14 4 6<pilcrow> to a network of slowly 
coupled oscillators. First, we use Malkin's theorem to derive a canonical 
phase model that describes synchronization properties of a slowly coupled 
network. Then, we illustrate the result using slowly coupled oscillators 
(1) near Andronov-Hopf bifurcations, (2) near saddle-node on invariant 
circle bifurcations, and (3) near relaxation oscillations. We compare and 
contrast synchronization properties of slowly and weakly coupled 
oscillators. 

English Descriptors: Coupled oscillators; Phase dynamics; Phase model; 
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Theory; Experiments 
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Simple model of spiking neurons 
IZHIKEVICH E M 
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A model is presented that reproduces spiking and bursting behavior of 
known types of cortical neurons. The model combines the biologically 
plausibility. of Hodgkin-yuxleyrtype* . dynamics^ - .an.d the . .computational 
efficiency of integrate-and-f ire neurons. Using this model, one can 
simulate tens of thousands of spiking cortical neurons in real time (1 ms 
resolution) using a desktop PC. 

English Descriptors: Spiking neurons; Hodgkin-Huxley model; Cortical 
neurons; Theory; Mathematical models; Computational methods; Computer 
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We demonstrate that the BCM learning rule follows directly from STDP when 
pre- and postsynaptic neurons fire uncorrelated or weakly correlated 
Poisson spike trains, and only nearest-neighbor spike interactions are 
taken into account. 
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Small-world 'an-d s-cale-free networks* are known *• "to be* more' easily 
synchronized than regular lattices, which is usually attributed to the 
smaller network distance between oscillators. Surprisingly, we find that 
networks with a homogeneous distribution of connectivity are more 
synchronizable than heterogeneous ones, even though the average network 
distance is larger. We present numerical computations and analytical 
estimates on synchronizability of the network in terms of its heterogeneity 
parameters. Our results suggest that some degree of homogeneity is expected 
in naturally evolved structures, such as neural networks, where 
synchronizability is desirable. 

English Descriptors: Theoretical study; Computerized simulation; 
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Motivated by the practical consideration of the measurement of chaotic 
signals in experiments or the transmission of these signals through a 
physical medium, we investigate the effect of filtering on chaotic symbolic 
dynamics. We focus on the linear, time-invariant filters that are used 
frequently in many applications, and on the two quantities characterizing 
chaotic symbolic dynamics: topological entropy and bit-error rate. 
Theoretical consideration suggests that the topological entropy is 
invariant under filtering. Since computation of this entropy requires that 
the generating partition for defining the symbolic dynamics be known, in 
practical situations the computed entropy may change as a filtering 
parameter is changed. We find, through numerical computations and 
experiments with a chaotic electronic circuit, that with reasonable care 
the computed or measured entropy values can be preserved for a wide range 
of the filtering parameter. (c) 2003 American Institute of Physics. 
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Efficiency in passage times is an important issue in designing networks, 
such as transportation or computer networks. The small-world networks have 
structures that yield high efficiency, while keeping the network highly 
clustered. We show that among all networks with the small-world structure, 
the most efficient ones have a single center node, from which all shortcuts 
are connected to uniformly distributed nodes over the network. The networks 
with several centers and a connected subnetwork of shortcuts are shown to 
be almost as efficient. Genetic-algorithm simulations further support our 
results . 
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We describe here how a network of oscillators can perform neural 
computations. In particular, it shown how the connectivity within the 
network can be created to memorize data in terms of phase relations between 
synchronized states. The memorized states are extracted through correlation 
calculations. The influence of noise on the system is discussed. 
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We suggest a simple spiking t mode^l-resonate-axKi-fAre neuron* which is 
similar to the ' integrate-and- fire neuron except that the state variable is 
complex. The model provides geometric illustrations to many interesting 
phenomena occurring in biological neurons having subthreshold damped 
oscillations of membrane potential. For example, such neurons prefer a 
certain resonant frequency of the input that is nearly equal to their 
eigenf requency, they can be excited or inhibited by a doublet (two pulses) 
depending on its interspike interval, and they can fire in response to an 
inhibitory input. All these properties could be observed in 
Hodgkin-Huxley-type models. We use the resonate-and- f ire model to 
illustrate possible sensitivity of biological neurons to the fine temporal 
structure of the input spike train. Being an analogue of the 
integrate-and-f ire model, the resonate-and-f ire model is computationally 
efficient and suitable for simulations of large networks of spiking 
neurons . 
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Periodic bursting behavior in neurons is a recurrent transition between a 
quiescent state and repetitive spiking When the transition to repetitive 
spiking occurs via a subcritical Andronov-Hopf bifurcation and the 
transition to the quiescent state occurs via fold limit cycle bifurcation, 
the burster is said to be of elliptic type (also known as a "subHopf / f old 
cycle" burster) . Here we study the synchronization dynamics of weakly 
connected networks of such bursters. We find that the behavior of such 
networks is quite different from the behavior of weakly connected phase 
oscillators and resembles that of strongly connected relaxation 
oscillators. As a result, such weakly connected bursters need few (usually 
one) bursts to synchronize, and synchronization is possible for bursters 
having quite different quantitative features. We also find that 
interactions between bursters depend crucially on the spiking frequencies, 
Namely the interactions are most effective when the presynaptic interspike 
frequency matches the frequency of postsynaptic oscillations. Finally, we 
use the FitzHugh-Rinzel , Morris-Lecar , and Hodgkin-Huxley models to 
illustrate our major results v , » . - . * 
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The transition to phase synchronization in systems consisting of a large 
number' (N)" of c'dupl'ed norilinear ' oscillators ' Via v the route" of phase 
clustering (phase synchronization among subsets of oscillators) is 
investigated. We elucidate the mechanism for the merger of phase clusters 
and find an algebraic scaling between the critical coupling parameter 
required for phase synchronization and N. Our result implies that, in 
realistic situations, phase clustering may be more prevalent than full 
phase synchronization. 
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We combine here two well-known and established concepts: 
microelectromechanical systems (MEMS) and neurocomputing. First, we 
consider MEMS ' oscillators "having low amplitude' activity arid we derive a 
simple mathematical model that describes nonlinear phase-locking dynamics 
in them. Then, we investigate a theoretical possibility of using MEMS 
oscillators to build an oscillatory neurocomputer having autocorrelative 
associative memory. The neurocomputer stores and retrieves complex 
oscillatory patterns in the form of synchronized states with appropriate 
phase relations between the oscillators. Thus, we show that MEMS alone can 
be used to build a sophisticated information processing system (U.S. 
provisional patent 60/178,654). 
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We use the Malkin theorem to derive phase equations for networks of 
weakly connected relaxation oscillators. We find an explicit formula for 
the connection* functions when the oscillators have one-dimensional slow 
variables. 'The functions are discontinuous in the relaxation limit mu 
rightarrow 0, which provides a simple alternative illustration to the major 
conclusion of the fast threshold modulation (FTM) theory by Somers and 
Kopell that synchronization of relaxation oscillators has properties that 
are quite different from those of smooth (nonrelaxation) oscillators. We 
use Bonhoef f er-Van Der Pol relaxation oscillators to illustrate the theory 
numerically. 
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We investigate here possible neurocomputational features of networks of 
laser oscillators. .Our approach is similar to classical optical 
neurocomputing where artificial neurons are lasers and connection matrices 
are holographic media. However, we consider oscillatory neurons 
communicating via phases rather than amplitudes. Memorized patterns 
correspond to synchronized states where the neurons oscillate with equal 
frequencies and with prescribed phase relations. The mechanism of 
recognition is related to phase locking. 
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We propose a .novej. . architecture of an oscillatory -neural network that 
consists of phase-locked loop (PLL) circuits. It stores and retrieves 
complex oscillatory patterns as synchronized states with appropriate phase 
relations between neurons. 
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Bursting behavior in neurons is a recurrent transition between a 
quiescent state and repetitive spiking. When the transition to repetitive 
spiking occurs via a subcritical Andronov-Hopf bifurcation and the 
transition to the quiescent state occurs via double limit cycle 
bifurcation, the burster is said to be of subcritical elliptic type. When 
the fast subsystem is near a Bautin (generalized Hopf) point, both 
bifurcations occur for nearby values of the slow variable, and the 
repetitive spiking has small amplitude. We refer to such an elliptic 
burster as being of local Bautin type. First, we prove that any such 
burster can be converted into a canonical model by a suitable continuous 
(possibly noninvertible ) change of variables. We also derive a canonical 
model for weakly connected networks of such bursters. We find that behavior 
of such networks is quite different from the behavior of weakly connected 
phase oscillators, and it resembles that of strongly connected relaxation 
oscillators. As a result, such weakly connected bursters need few (usually 
one) bursts to synchronize. In-phase synchronization is possible for 
burst ers having quite different quantitative featpur^s^ , whereas out^— of-phase 
synchronization may be difficult to achieve. We also find that interactions 
between bursters depend crucially on the spiking frequencies. Namely, the 
interactions are most effective when the presynaptic interspike frequency 
matches the frequency of postsynaptic oscillations. Finally, we use the 
FitzHugh-Rinzel model to evaluate how studying local Bautin bursters can 
contribute to our understanding of the phenomena of subcritical elliptic 
bursting. 
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We prove that weakly connected networks of quasi-periodic 
(multif requency) oscillators can be transformed into a phase model by a 
continuous change of variables. The phase model has the same form as the 
one for periodic oscillators with the exception that each phase variable is 
a vector. When the oscillators have mutually nonresonant frequency 
(rotation) vectors, the phase model uncouples. This implies that such 
oscillators do not interact even though there might be physical connections 
between them. When the frequency vectors have mutual , low- order resonances, 
the oscillators interact via phase deviations. This mechanism resembles 
that of the FM radio, with a shared feature - multiplexing of signals. 
Possible applications to neuroscience are discussed. 
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We study pulse-coupled neural networks that satisfy only two assumptions: 
each isolated neuron fires periodically, and the neurons are weakly 
connected. Each such network can be transformed by a piece-wise continuous 
change of variables into a phase model, whose synchronization behavior and 
oscillatory associative properties are easier to analyze and understand. 
Using the phase model, we can predict whether a given pulse-coupled network 
has oscillatory associative memory, or what minimal adjustments should be 
made so that it can acquire memory. In the search for such minimal 
adjustments we obtain a large class of ^simple pulse-coupled neural networks 
that can memorize and reproduce synchronized temporal patterns the same way 
a Hopfield network does with static patterns. The learning occurs via 
modification of synaptic weights and/or synaptic transmission delays. 
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Many scientists believe that all pulse-coupled neural networks are toy 
models that are far away from the biological reality. We show here, 
however, that a huge class of biophysically detailed and biologically 
plausible neural-network models can be transformed into a canonical 
pulse-coupled form by a piece-wise continuous, possibly noninvertible, 
change of variables. Such transformations exist when a network satisfies a 
number of conditions; e.g., it is weakly connected; the neurons are Class 1 
excitable (i.e., they can generate action potentials with an arbitrary 
small frequency); and the synapses between neurons are conventional (i.e., 
axo-dendritic and axo-somatic) . Thus, the difference between studying the 
pulse-coupled model and Hodgkin-Huxley-type neural networks is just a 
matter of a coordinate change. Therefore, any piece of information about 
the pulse-coupled model is valuable since it tells something about all 
weakly connected networks of Class 1 neurons. For example, we show that the 
pulse-coupled network "of identical 4 neurons does hot 'synchronize " ih-phase . 
This confirms Ermentrout's result that weakly connected Class 1 neurons are 
difficult to synchronize, regardless of the equations that describe 
dynamics of each cell. 

English Descriptors: Neural network; Synapse; Desynchronization ; Biophysics 
; Axodendritic synapse; Axosomatic synapse; Canonical model; Hodgkin 
neural network 

French Descriptors: Reseau neuronal; Synapse; Desynchronisation; 
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Modele phase; Modele canonique; Reseau neuronal Hodgkin 
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Our study of thalamo-cortical systems suggests a new architecture for a 
neurocomputer that consists of oscillators having different frequencies and 
that are connected weakly via a common medium forced by an external input. 
Even though such oscillators are all interconnected homogeneously, the 
external input imposes a dynamic connectivity. We use Kuramoto<right single 



quotation mark>s model to illustrate the idea and to prove that such a 
neurocomputer has oscillatory associative properties. Then we discuss a 
general case. The advantage of such a neurocomputer is that it can be built 
using voltage controlled oscillators, optical oscillators, lasers, 
microelectromechanical systems, Josephson junctions, macromolecules, or 
oscillators of other kinds. (Provisional patent 60/108,353) 
English Descriptors: Instrumentation; Theoretical study; Brain; 
Biocomputers ; Neurophysiology; Oscillators; Neural networks 
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The cusp bifurcation provides one of the simplest routes leading to 
bistability and hysteresis in neuron dynamics. We show that weakly 
connected networks of neurons near cusp bifurcations that satisfy a certain 
adaptation condition have quite interesting and complicated dynamics. 
First, we prove that any such network can be transformed into a canonical 
model by an appropriate continuous change of variables. Then we show that 
the canonical model can operate as a multiple attractor neural network or 
as a globally asymptotically stable neural network depending on the choice 
of parameters . 

English Descriptors: Neural network; Connectionism; Multiple system; Cusp 
configuration; Bifurcation; Pitching; Canonical ensemble; Learning; 
Bistability; Perception; Theoretical study 
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Studying weakly connected oscillators leads to phase models. It has been 
proven recently that weakly connected oscillators with delayed interactions 
do not lead to phase models with time delays even when the delay is of the 
same order of magnitude as the period of oscillation. This has resulted in 
a fading interes/t iri . such models. Ke prove here -tha*t • if the interaction 
delay between weakly connected oscillators is much longer than the period 
of oscillation, then the corresponding phase model does have an explicit 
time delay. 

English Descriptors: Theoretical study; Oscillations; Nonlinear dynamical 
systems 

French Descriptors: 8710; 0545; 0705M; 4279T; Etude theorique; Oscillation; 
Systeme dynamique non lineaire 



Classification Codes: 002A03A; 001B00E45; 001B00G05M; 001B40B79T 
Copyright (c) 1998 American Institute of Physics. All rights reserved. 



1/5/180 (Item 25 from file: 144) 

DIALOG (R) File 144: Pascal 

(c) 2004 INIST/CNRS. All rts . reserv. 

12520919 PASCAL No.: 96-0194943 

An averaging principle for dynamical systems in Hilbert space with Markov 
random perturbations 
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Journal: Stochastic processes and their applications, 1996, 61 (1) 
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Copyright (c) 1996 Elsevier Science B.V. All rights reserved. We study 
the asymptotic behavior of solutions of differential equations dx SUB 
epsilon (t) /dt=A (y (t/ epsilon ))x SUB epsilon (t) , x SUB epsilon (0)=x SUB 
0 , where A(y) , for y in a space Y, is a family of operators forming the 
generators of semigroups of bounded linear operators in a Hilbert space H, 
and y(t) is an ergodic jump Markov process in Y. Let A= integral A(y) rho 
.(dy) where rho (dy). is. the er.godia distribution o.f, y.( t) — We show that under 
appropriate conditions as epsilon rightarrow 0 the process x SUB epsilon 

(t) converges uniformly in probability to the nonrandom function &xmacr; (t) 
which is the solution of the equation d&xmacr ; ( t ) /dt=A&xmacr ; ( t ) , 
&xmacr;(0)=x SUB 0 and that epsilon SUP - SUP 1 SUP / SUP 2 (x SUB epsilon 

( t ) -&xmacr ; (t ) ) converges weakly to a Gaussian random function &xtilde; (t) 
for which a representation is obtained. Application to randomly perturbed 
partial differential equations with nonrandom initial and boundary 
conditions are included. 

English Descriptors: Dynamical system; Stochastic system; Averaging method; 
Asymptotic expansion; Partial differential equation; Hilbert space; 
Bounded operator; Markov process; Jump process; Random function 
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English Descriptors: Stochastic equation; Volterra integral equations; 
Volterra equation 

French Descriptors: Equation stochastique; Equation integrale Volterra; 
Equation Volterra 
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Recent work on singular perturbation solutions that persist in the 
presence of noise is described. Two different settings are considered: 
small deviation theory in quasi-static problems, where there are small 



amplitude but highly irregular perturbations, and averaging problems where 
there are ergodic stochastic perturbations. In the first case, it is shown 
that quasi-static approximations can be valid when the underlying problem 
experiences small deviation perturbations in problems that are stable under 
persistent disturbances. In the second, averaging principles are described 
for certain dynamical systems in Hilbert spaces that include applications 
to a wide variety of initial-boundary value problems for partial 
differential equations and for Volterra integral equations. These methods 
are applied here to four problems arising in applications 
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Reponses en frequence des neurones, modele du type bscillateur commande 
par la tension, reseaux de neurones 
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On demontre 1 ! existence de solutions periodiques pour des equations 
integrodif f erentielles de type Volterra 
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On montre comment la forme de 1 1 approximation perturbation pour la 
resolution de systemes d 1 equations dif f erentielles a petit parametre 
identifiable peut s'utiliser pour generer des equations non stiff ou 
•relaxees . • - > . . • - • - - 
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On deduit un algorithme de perturbation multitemps pour etudier des 
systemes de comrnande synchrones et des processus de renouvellement non 
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Le modele de Fitzttugh-Nagumo permet de decrire de facon simple et 
satis faisante les horloges radiales isochrones (RIC) qui sont etudiees ici 
lorsqu'elles sont soit isolees soit couplees. Modelisation des ruptures de 



rythmes observees chez les poissons, les oiseaux et les mammiferes et 
consequents a des modifications de l'intensite lumineuse 
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Au laboratoire, en modifiant la photoperiode, on peut provoquer, chez 
M.m. , le passage, d'une activite nocturne a une t .actlvite d4urne..Dans la 
nature, il s'agit peut etre d'un mecanisme pour eviter les micro-climats 
trop rigoureux 
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English Descriptors: DIFFUSION EQUATION; DIFFUSION PROCESS; EQUATION SYSTEM 

; NON LINEAR THEORY; TRANSPORT THEORY 
English Generic Descriptors: THEORETICAL PHYSICS 

French Descriptors: THEORIE TRANSPORT; PROCESSUS DIFFUSION; THEORIE NON 

LINEAIRE; EQUATION DIFFUSION; SYSTEME EQUATION 
French . Generic Descriptors : PHYSIQUE. THEORIQUE . . . - . 



Classification Codes: 130A05E 



1/5/204 (Item 49 from file: 144) 

DIALOG (R) File 144: Pascal 

(c) 2004 INIST/CNRS. All rts. reserv. 

01777096 PASCAL No.: 78-0002124 

MULTITIME METHODS FOR SYSTEMS OF DIFFERENCE EQUATIONS. 

HOPPENSTEADT F C ; MI RANKER W L 
Journal: STUD. APPL. MATTR. , 1977, 56 (3) 273-289 
Availability: CNRS-546 
No. of Refs . : 5 REF. 

Document Type: P (SERIAL) ; A (ANALYTIC) 
Country of Publication: USA 
Language: ENGLISH 

ON ETUDIE DES SYSTEMES D 1 EQUATIONS AUX DIFFERENCES CONTENANT DE PETITS 
.PARAMETRES. .A L'ALDE . D 1 UN ...SCHEMA .DEL PERTURBATION- CON STRUCT I F ANALOGUE A 
CELUI DEVELOPPE POUR L 1 ETUDE DES EQUATIONS DI FFERENTIELLES . 

English Descriptors: EQUATION; DIFFERENCE EQUATION; EQUATION SYSTEM; 

PERTURBATION THEORY 
English Generic Descriptors: MATHEMATICS 

French Descriptors: EQUATION; EQUATION DIFFERENCES; SYSTEME EQUATION; 

THEORIE PERTURBATION 
French Generic Descriptors: MATHEMATIQUES 



Classification Codes: 110A02A 



1/5/205 (Item 50 from file: 144) 

DIALOG (R) File 144: Pascal 

(c) 2004 INIST/CNRS. All rts. reserv. 

01512853 PASCAL No.: 77-0245496 
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LES ECLOSIONS SYNCHRONISEES SE MANIFESTENT POUR DES INSECTES DONT LE 
CYCLE EST DE 10 ANS ET PLUS ET LES ECLOSIONS NON SYNCHRONISEES POUR LES 
INSECTES DONT LE CYCLE EST INFERI EUR A 10 ANS. 
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ON ETUDIE LES SOLUTIONS DE L 1 EQUATION AUX DIFFERENCES X(N+1)=MX(N) X 
(l-X(N)) POUR M COMPRIS ENTRE 0 ET 4 . { ON ETUDIE LE COMPORTEMENT DES 
SOLUTIONS EN REGIME CHAOTIQUE. ^ FINALEMENT ON CALCULE NUMERIQUEMENT t LES 
FONCTIONS " DE DENSITE POUR DECRIRE LE^COMPORTEMENT DYNAMIQUE *DES * SOLUTIONS 
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The goal of this book chapter is to present several interesting, 
canonical types of bifurcations (including cusp, fold and Andronov-Hopf ) 
that occur in continuous and discrete time models of population biology. 
These include models of epidemics, fisheries and predator-prey systems 
(without and with random perturbations) . As stated by the authors, the 
selection of the models reflects their interests. They give references to 
the literature in which the models are described in more detail, but do not 
attempt to -give' pro'ofs* of results br' a "comprehensive" preference" list'. The 
authors have succeeded admirably in their goal. The chapter is recommended 
reading for anyone wanting to find out about these interesting bifurcations 
in the context of population biology. 

\{For the entire collection see MR 2004c : 00011 . \ } 
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Y. Kuramoto [Chemical oscillations, waves, and turbulence, Springer, 
Berlin; 198'4"; MR 8*7e: 9*2054] showed* that" in' his mean- field mod£l consisting ' 
of $N$ coupled oscillators, in the limit as $N \to \infty$ there exists a 
critical coupling strength $\epsilon\sb c$ such that for $\epsilon < 
\epsilon\sb c$ the solution is incoherent but for $\epsilon > \epsilon\sb 
c$ partially coherent solutions appear with a nonzero fraction of locked 
pairs of oscillators. In the present paper the authors consider a class of 
$N$ phase oscillators with randomly distributed natural frequencies and 
diagonalizable interactions between the oscillators. For this class of 
system complete separation of variables through an appropriate change of 
variable is possible. Some properties of diagonalizable systems are 
established. The behaviour of a generic diagonalizable system in the limit 
as $N \to \infty$ is described. It is shown that in the limit as $N \to 
\infty$, all solutions of above types of systems are incoherent with 
probability one for any strength of coupling. The implication of this 
result is that there is no "sharp transition from incoherence to coherence 
as the coupling strength is increased in diagonalizable systems. 

Reviewer: Rajasekar, S. (Tirunelveli ) 
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Random perturbation methods with applications in science and 
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Abstract Length: LONG (77 lines) 
In the first half of their book (Chapters 1 — 7) the authors introduce 
the reader to the mathematical theory of randomly perturbed dynamical 
systems. They use the tools of this theory to study in detail the dynamics 
of real-world systems subject to random impact which are introduced in the 
second ■ half' -of the 'book (Chapters •*8-*- 12*) .* The systems v c"dnside red ' here 
appear in mechanics (conservative systems with two degrees of freedom, 
dynamical systems on the torus (Chapters 8 and 9)), in quantum mechanics 
and electronics (the phase-locked loop, a standard electronic circuit 
(Chapter 10) ) , in population biology (predator-prey type systems (Chapter 
11)) and in genetics (gene pool, plasmid stability for bacteria, evolution 
of the genome (Chapter 12)). 

The common feature of these systems is that they can be modeled as 
Volterra integral equations which, mostly, can be reduced to ordinary 
differential equations, or as difference equations with randomly varying 
parameters. These parameters in turn are modeled as ergodic stationary or 
ergodic Markov processes $y ( t/\epsilon) $ which vary at a much faster 
time-scale than the unperturbed system. A randomly perturbed system of thi 
kind is close, in a certain sense, to a deterministic system that is 
obtained by averaging out the fast moving noise (first order approximation 
Chapter 3: Averaging). Typically, if $\epsilon > 0$ is very small, the 
solution $x\sb { \epsilon} ( t ) $ of the differential equation $\dot x\sb 
{ \epsilon} (t ) = a(x\sb { \epsilon ) ( t ) , y ( t/\epsilon) ) $ is close to the 
solution $\dot { \overline{x) } (t) $ of the deterministic ODE 
$\dot{\overline{x) } (t) =\overline { a ) (\dot { \overline { x} } (t) ) $, where 
$\overline{a-J (x) = *\irrt { a (x, y) d\rh*o (y) • $\rho$ ' the' ergodic measure 
associated with the stochastic process $y ( t/\epsilon) $ . 



The second order approximation is due to the fact that the suitably 
scaled differences, e.g. $\tilde { x} \sb { \epsilon} ( t ) = ( l/\epsilon) [x\sb 
{ \epsilon} ( t ) - \overline {x} (t) ] $, converge to a Wiener process in the weak 
sense either with respect to only all finite-dimensional distributions, 
called weak convergence, or, stronger, with respect to the distribution of 
the. random function $ Vtilde { x } \sb >{ \epsdlon } ( \cdot*) $• wi-fch values in • the - 
infinite-dimensional space of continuous functions, called weak convergence 
in $C$ (Chapter 2 and Chapter 4: Normal Deviations). These ideas are 
analogous to the convergence properties of real-valued random variables 
(Law of Large Numbers and Ergodic Theorems, Chapter 1, and the Central 
Limit Theorem (CLT) ) and they include CLTs for strongly mixing Markov 
processes as well as for stationary and ergodic processes. 

In Chapter 5 (Diffusion Approximation) the authors give conditions under 
which the speeded up state $x\sb { \epsilon} ( t/\epsilon) $ of the randomly 
perturbed system can be approximated by the solution of a suitable 
diffusion equation, while in the comparatively long Chapter 6 (Stability) 
they investigate the long term stability properties (exponential growth 
rates or Lyapunov exponents) of randomly perturbed difference equations, 
ordinary differential equations and convolution integral equations as well 
as large deviations of systems starting near a stable steady state and 
noise induced oscillations between two steady states (stochastic 
resonance) . 

The theoretical part closes with a short introduction to the theory of 
Markov chains in random environments (Chapter 7) which provides the tools 
for studying the evolution of the genome in a random environment (Chapter 
12, . Section, 3 ) . , . . - . * ■ • - ■ * . 

The book, published in the Springer series of Applied Mathematical 
Sciences, is written in mathematical terms. But the authors make the 
laudable attempt to introduce others, e.g. engineers and biologists (and 
non-expert mathematicians) , to the area of randomly perturbed dynamical 
systems by starting the book with an Introduction that in less technical 
terms illustrates the main facts of the theory and how these facts help one 
better understand the dynamical systems considered by scientists and 
engineers. This, together with many computer-simulated visualizations, may 
make the book accessible to scientists and engineers with a general 
mathematical background. On the other hand, the chapters in the applied 
part of the book provide a knowledge of biological facts that are needed by 
a mathematician to understand the mathematical model of the real world. 

The comments referring to the sources are sufficient, though very short. 
For a second edition a larger index and correction of the misprints would 
be appreciated. 

Reviewer: Wihstutz, Volker (1-NC3) 
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This is an introductory book on mathematical modeling in the 
bio-sciences. It is written for mathematicians as well as for life 
scientists. Simple models are presented, and previous knowledge of biology 
is not required for understanding the book. All the essential biological 
background. . i.s given in .the text, .jwhile. basic mathematical knowledge is. 
sufficient for reading a large part of the book. 

In each chapter, the material is organized in increasing order of 
complexity followed by exercises. Some of the exercises deal with the 
material of that chapter, while others are projects that extend the 
preceding material. Many chapters contain sections with suggestions for 
computing projects. Simulations are done in Matlab and computer code is 
included in the text. 

There are two major parts of the text. The first part deals with models 
from physiology, while the second part deals with models from population 
biology. In Chapter One, the authors derive models of blood flow and 
pressure, and mechanisms for controlling them. They describe a model of 
uncontrolled circulation, the baroreceptor loop, as neural control of 
circulation, and autoregulation . At the end of the chapter, they consider 
the dynamics of the arterial pulse, with nonconstant pressure, flow and 
volume in time. 

Gas exchange in the lungs is described in Chapter Two. Equations of gas 
transport in one alveolus are derived from the ideal gas law. They are used 
to describe gas transport in the lung. Chapter Three deals with control 
of cell volume. Based on osmotic effects and neglecting all electrical 
effects, a. sample modej. .of cell volume .^control is. derived. By .introducing 
electrical effects, a more complex model is obtained. A special section 
covers the Hodgkin-Huxley equation for the nerve action potential. Chapter 
Four describes the renal mechanism. Here the dynamics of sodium ions and 
water in nephrons are modeled. Muscle mechanics is explained by a 
force-velocity curve and crossbridge dynamics in Chapter Five. Chapter Six 
deals with neural systems and the main topic is a model of a neural 
network. Population dynamics is considered in Chapter Seven. Bacterial 
cultures are used to illustrate a simple model of exponential growth. In a 
more complex example of age structures, Euler ! s renewal theory is 
presented, and in an example of microbial ecology, more complex continuous 
models of population growth are presented. Nonlinear reproduction curves 
and controlling populations are other topics in this chapter. There are 
two topics from genetics in Chapter Eight: population genetics and an 
application of mathematics in biotechnology. 

Epidemic models are described in Chapter Nine. The authors present the 
spread of an infection within a family, the threshold of an epidemic and 
prediction of the severity of an epidemic. In the last chapter, they 
describe modeling of patterns of population growth. 
Reviewer: Marusic, Miljenko (CT-ZAGR) 
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Summary: ' 'The diffusion gradient chamber (DGC) is a novel device 
developed to study the response of chemotactic bacteria to combinations of 
nutrients and attractants [D. Emerson, R. M. Worden and J. A. Breznak, 
Appl . Environ. Microbiol. 60 (1994), no. 4, 1269]. Its purpose is to 
characterize genetic variants that occur in many biological experiments. In 
this paper, a mathematical model which describes the spatial distribution 
of a bacterial population within the DGC is developed. We give a 
mathematical analysis of the model concerning positivity and boundedness of 
the solutions. An alternating direction implicit method is constructed for 
finding numerical solutions of the model and carrying out computer 
simulations. The numerical results of the model successfully reproduce the 
patterns that are observed in the experiments using the DGC. 1 r 
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Abstract Length: MEDIUM (13 lines) 

Summary: N s We combine here two well-known and established concepts: 
microelectromechanical systems (MEMS) and neurocomputing. First, we 



consider MEMS oscillators having low amplitude activity and we derive a 
simple mathematical model that describes nonlinear phase-locking dynamics 
in them. Then we investigate the theoretical possibility of using MEMS 
oscillators to build an oscillatory neurocomputer having autocorrelative 
associative memory. The neurocomputer stores and retrieves complex 
oscillatory patterns in the form of synchronized states with appropriate 
phase relations between the oscillators. Thus, we show that MEMS alone can 
be used to build a sophisticated information processing system (U.S. 
provisional patent 60/178 , 654 ). 1 ' 
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The first edition of this book has been reviewed [MR 94a:34003]. 
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The subtitle of this book, v 'Modeling in the frequency domain'', reflects 
its focus on a neuronal model called ' x VCON'', or voltage controlled 
oscillator neuron' based on an analogy between neurons and phase-locked 
loops of voltage-controlled oscillators in electronic circuits. Informally, 
.the idea is, .that th.e neuronal axon, hillock, generates, auction potentials of 
relatively constant waveform; what varies depending on inputs to the neuron 
is the timing of the action potential. Therefore the output can be written 
$V (\phi ( t ) ) $ , where $V(\cdot)$ is a fixed function (the waveform), and 
$\phi(t)$ is the phase. Modeling of a neuronal network is said to be in the 
N ' frequency domain' 1 if the dynamics of the network can be expressed in 
terms of differential equations with the phases as the variables . 
Phase-locking behavior is sought, i.e. conditions under which the neurons 
in a network tend toward constant phase differences in the course of time. 
The analogy is an interesting and fruitful one. 

The book begins with a review of electrical circuits, focusing on 
voltage-controlled oscillators (VCOs), where the time derivative of the 
output phase is a linear function of the input voltage. A phase-locked 
loop consists of a phase detector, a low-pass filter and a VCO in series, 
with feedback from the output of the VCO to the phase detector. After a 
review of other models of neuronal firing, including the Hodgkin-Huxley ion 
channel model, the VCO-VCON analogy is developed. A chemical synapse can 
be modeled as a low-pass filter, obeying the relationship $\tau 
\dot {X}+X=S\sb X$, where $X$ is the postsynaptic potential and $\tau$ is 
the reciprocal of the sum of the diffusion and postsynaptic binding 
constants fp.r tr;e .neurotransmitter,, replacing the. time.. , const ant $RC$ of, a 
resistive-capacitative filter. $S\sb X$ represents the source of the 
neurotransmitter, and is proportional to $\dot { \phi } $ , where $\phi$ is the 
input phase, since the time derivative of the input phase is related to the 
rate of arrival of presynaptic impulses. Electrical synapses can be 
similarly modeled. The phase detector of the VCO is replaced by a mixer 
that combines the postsynaptic potentials with the cell's current state. 
The mixing function is typically sigmoidal. VCONs may have multiple stable 
states, and can also undergo a saddle-node bifurcation in which the 
relative output phase no longer tends to a constant, but increases without 
bound, as the neuron fires repetitively. This bifurcation process is 
discussed in more mathematical detail in another recent book by the author 
and E. M. Izhikevich [Weakly connected neural networks, Springer, New York, 
1997; MR 98k:92004] . 

In later chapters many applications of the VCON model to small and large 
neuronal networks are given. Small networks include, for example, the 
' x atoll model 11 , in which an excitatory and an inhibitory cell interact in 
such a way that a burst in one cell is followed by a slow pulse in the 
other. The atoll oscillator model is later applied on a larger scale, to 
control of attention, through a model involving excitatory cells in the 
thalamus interacting with inl^ibito.ry, cej-ls. in the. thalamic reticular 
complex. Applications to central pattern generators such as the controller 
of the gastric mill in Crustacea, flight in moths, respiration, and 
binaural sound location are given. Larger neuronal networks are analyzed 
through continuum approximations, and also discussed in terms of v "mnemonic 
surfaces' 1 , energy-like functions whose minima correspond to stable choices 
for the relative phases of neurons in the network, and are capable of 
encoding information. 

Reviewer: Matthysse, Steven (Belmont, MA) 
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The asymptotic behavior of Markov chains with small perturbations of the 
$n$-step transition probabilities as follows is considered: $P\sp {*}\sb n 
(\epsilon) =P+ \epsilon Q\sp {*}\sb n$, where $P$ is a transition matrix 
and $Q\sp {*}\sb n$ is an ergodic stationary matrix-valued process, and 
$\epsilon$ is a small parameter. It is proved that the $n$-steps 
transition matrix of the Markov chain above converges to a nonrandom matrix 
with probability 1 under some general conditions as $n$ goes to infinitity 
and $\epsilon$ goes to 0. The result above is applied to the evolutionary 
path of bacteria. Similar considerations in more detail can be found in 
another paper by the authors [Random Oper. Stochastic Equations 4 (1996), 
no. 3, 205--227; MR 97j:60127], 

\{For the entire collection see MR 98b:62002.\} 
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The book under review is devoted to an analysis of general continuous, 
weakly connected neural networks of the form (1) $\dot x\sb i=f\sb i(x\sb 
i)+\epsilon g\sb i (x\sb 1, \cdots, x\sb n, \epsilon) $ . Here $x\sb i$ are 
"neural ■ activity 'vectors,' $x\Sb i\iVi \bo*ld' R\sp m$' of* $'n$ neurons , * 
$i=l, \cdots , n$ ; $f\sb i$ describes the dynamics of the $i$th neuron, and 
$g\sb i$ describes interactions between neurons. The small parameter 
$\epsilon$ indicates the strength of synaptic connections. The method used 
is that of treating (1) as an $\epsilon$-perturbation of the uncoupled 
system with $\epsilon=0$ . Then a bifurcation analysis is performed in 
detail for two appropriate models which fulfill the conditions of (1), 
namely the additive neural network model with a sigmoidal transfer function 
and the Wilson-Cowan model for the pair of excitatory and inhibitory 
neurons, creating a neural oscillator. In the book only limit cycles of 

weakly connected networks of resting neurons equilibrium points, and 

periodically spiking neurons are analyzed. As the authors say in the 

nicely written preface: N ^While it is feasible to study weakly connected 
networks of chaotic neurons, this problem goes beyond the scope of this 
book. 11 

After the introduction, written according to the authors in ordinary 
language, and well readable even for laymen, follows a nicely written 
Chapter 2 on bifurcations in neuron dynamics which must be read. Here also 
spiking and bursting phenomena are clearly described. Chapter 3 contains a 
short sketch of nonhyperbolic (when the Jacobian matrix of (1) has at least 
one- eigenvalue with zero real part*) netfrai networks'.' " The ' remaining part 'of' 
the book is mainly devoted to canonical models (Chapter 4), their 
derivation (Chapters 6--9) , and their analysis (Chapters 10--12). The term 
canonical model is not precisely defined here. The authors say that a model 
is canonical if there is a continuous change of variables that transforms 
any other model from a given class into this one. As the method of deriving 
the canonical models, the authors exploit the normal form theory. Canonical 
models treated in the book have only restricted value: They provide 
information about local behavior of (1) when there is an exponentially 
stable limit cycle but they say nothing about global behavior of (1), 
including the transients. The last Chapter 13 describes the relationship 
between synaptic organizations and dynamical properties of networks of 
neural oscillators. In other words, the problem of learning and 
memorization of phase information in the weakly connected network of 
oscillators corresponding to multiple Andronov-Hopf bifurcation is treated 
analytically . 

Surprisingly the book ends without any conclusions. Also there are no 
appendices to the book. The references are representative and sufficiently 
cover the problematics treated in the book. 
Reviewer: Andrey, Ladislav (Prague) 
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In this paper the authors consider a system in a linear phase space $X$ 
with discrete time $n\in\bold Z\sb +$, which is perturbed by a random 
process $\{y\sb n,\ n\in\bold Z\sb +\}$ defined on a measurable space 
$(Y,\scr C)$. The system depends on a small positive parameter $\epsilon$. 
The. state of- the s.ystem .at time $n$,.$»\sp \epsilon\sb *n\in X$, is - - - 
determined by the recurrence relations $x\sp \epsilon\sb {n+l}=f (x\sp 
\epsilon\sb n) +\epsilon\phi ( x\sp \epsilon\sb n,y\sb {n+l})$, for $n\geq 
0$, $x\sp \epsilon\sb 0=x\sb 0$, where $x\sb 0$ is given, $f\colon X\to X$ 
is a continuous function and $\phi\colon X\times Y\to X$ is a function 
which is continuous in $x$ and measurable in $y$. The system is studied in 
the following phase spaces: (1) $X$ a separable Banach space; (2) $X$ a 
separable Hilbert space; (3) $X=\bold R\sp d$ . The authors investigate the 
asymptotic behavior of the system as $\epsilon\to 0$ and $n\to\inf ty$ , and 
give applications of the results obtained to epidemics, to genetics and to 
demography. 
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The authors consider the equations $$x\sb \epsilon ( t ) =\phi ( t ) +\int\sp 
t\sb 0k\bigg (t, s, y\bigg (\f rac s\epsilon\bigg) (s ) \bigg) ds\quad{ \rm and}\tag 
1$$ $$\overline x (t ) =\phi ( t) +\int\sp t\sb 0\overline k (s , t, \overline 
x(t))ds,\tag 2$$ where $y(s)$ is an ergodic stationary process, or a 
stationary ergodic Markov process. Let $\rho(dy)$ denote the ergodic 
distribution of $y(t)$ and let $\overline k ( s , t, x) =\int 
k (s, t, y, x) \rho (dy) $ . First they find sufficient conditions for the 
convergence of $x\sb \epsilon(t)$ to $\overline x(t)$ when $\epsilon\to 0$, 
more precisely, for the relation $P\ { \lim\sb {\epsilon\to 0}\sup\sb {t\leq 
T(\Vert x\sb \epsilon ( t ) -\overline x(t)\Vert =0\}=1$ to be true for all 
$T>0$. They also show that, under some smoothness conditions on the kernel 
with respect to $x$, the solution $x\sb \epsilon(t)$ of (1), with 
$\phi ( t ) =x\sb 0$, converges weakly, when $\epsilon\to 0$, to a process 
$x(t)$ which satisfies the stochastic integral equation (3) $x(t)=x\sb 
0+\int\sp tVsb Oxfs')dw(s ) $ . In (1)", $w*ft)$ is an'^tfAbold R\sp d) $-valued 
Wiener process for which $Ew(t)=0$ and the covariance is expressed by means 
of $\rho(dy)$ and $k ( y) =\lim\sb { s , t\to\inf ty, t-s\to\inf ty } k ( s , t , y ) $ . As 
techniques, the authors show that $x(t)$ may be represented by a series of 
multiple stochastic integrals, $$x ( t ) =\phi ( t ) +\sum\sp \infty\sb 
{n=l}\int\sb {0<s\sb K\cdots<s\sb n<t}d\sb {s\sb n}w(s\sb n,t), 
\cdots,d\sb {s\sb l}w(s\sb l,s\sb 2)\phi(s\sb 1),$$ where $\phi(t)$ and 
$x(t)$ are $\bold R\sp d$-valued processes for all $t\in\bold R\sb +$, and 
$w(s,t)$ is an $L(\bold R\sp d) $-valued Wiener process in $s$, on the 
interval $[0,t]$. The authors treat particularly interesting problems that 
arise for integral equations $$x\sb \epsilon ( t) =\phi ( t) +\int\sp t\sb 
0M\bigg ( t-s, y\bigg ( { s\over\epsilon } \bigg) , x\sb \epsilon (s) \bigg) ds\tag4$$ 
of convolution type. These equations are treated by introducing Laplace 
transform techniques. The authors also discuss the applications of these 
perturbed Volterra equations in a number of domains: epidemics, 
demographics and electrical engineering. This is an interesting and 
valuable paper in both the theoretical and applied fields. 
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The convergence of the Cesaro averages of $n$-step transition 
probabilities is proved for finite Markov chains in random environments 
with the following $n$-step transition probability matrix: $P\sp {*}\sb n 
(\epsilon ) =P+\epsilon Q\sp {*}\sb n$, where $P$ is a transition matrix 
and the entries of $P\sp {*}\sb n(\epsilon )$ are stationary ergodic 
sequences and, as $\epsilon$ goes to zero, the limit distribution goes to 
the ergodic limit distribution of the homogeneous Markov chain with the 
transition matrix of the the mean of $P\sp {*}\sb n(\epsilon )$. The result 
above is applied to the evolutionary path of bacteria, i.e., to define and 
to express the fitness of a chromosome evolving in a random enviroment in 
terms that are observable in experiments. About the products of random 
matrices readers may also refer to J. E. Cohen [Bull. Amer. Math. Soc. 
(N.S.) 1 (1979), no. 2, 275—295; MR 81j:92029]. 
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A typical result of this paper is the theorem on the averaging principle 
for stochastic differential equations of jump type in Hilbert space. 
Consider the solution $x\sb \epsilon(t)$ of the differential equations 
$dx\sb \epsilon (t) /dt=A(y (tAepsilon) ) x\sb \epsilon ( t ) $ , $x\sb 
\epsilon (0) =x\sb 0$, $0<\epsilon\ll 1$, where $A(y)$, for $y$ in a space 
$Y$, is a family of operators forming the generators of semigroups of 
bounded linear operators in a Hilbert space $H$, and $y(t)$ is an ergodic 
jump Markov process in $Y$ with its ergodic probability distribution 
$\rho(dy)$. Further, let $\overline x(t)$ be the solution of the equation 
$d\overline{x} ( t ) /dt=\overline { A} \ over line { x } (t)$, $\overline { x } (0)=x\sb 
0$, where $\overline A=\int A ( y ) \rho (dy ) $ . Then the authors find conditions 
under which $x\sb \epsilon ( t ) \to\overline { x } ( t ) $ in probability uniformly 
on all finite intervals as $\epsilon\to 0$, and they also prove that, under 
some additional conditions, $\epsilon\sp {-l/2}(x\sb_ 
\epsilon (t) -\overl'ine fx}' (t) ) £ converges weakly in distribution as 
$\epsilon\to 0$ to a Gaussian process $\tilde x(t)$ with independent 
increments in $t$ for which the limiting distribution is obtained. They 
give applications of these results to partial differential equations with 
random perturbations and take examples arising in various fields, such as 
bacterial growth, random advection and diffusion in random media. 
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In this paper the author gives a brief survey of recently obtained 
results of his and coauthors [F. C. Hoppensteadt , Analysis and simulation 
of chaotic systems, Springer, New York, 1993; MR 94a: 34003; F. C. 
Hoppensteadt, R. Z. Khasminskii and H. Salehi, Random Oper. Stochastic 
Equations 2 (1994), no. 1, 61 — 78; MR 95d: 35189; F. C. Hoppensteadt, H. 
Salehi and A. V. Skorokhod, N 'An averaging principle for dynamical systems 
in Hilbert space with Markov random perturbations 1 ' , Stochastic Process. 
Appl., to appear; F. C. Hoppensteadt, H. Salehi and A. V. Skorokhod, 
N 'Randomly perturbed Volterra integral equations and some applications 1 1 , 
Stochastics Stochastics Rep., to appear]. He describes recent work on 
singular perturbation solutions that persist in the presence of noise in 
two different settings: one is the small deviation theory in quasi-static 
problems where there are small-amplitude but highly irregular 
perturbations, and the other is that of averaging problems where there are 
ergodic stochastic perturbations. He states that his new methods can be 
applicable to topics, such as Lorenz's equations, controls, diffusions in 
randomly fluctuating media and epidemic models. 
\{For the entire collection see MR 95j:00015\}. 
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also 34K50, 60H10, 93E03) ; 35R60 -Partial differential equations- 
Miscellaneous topics involving partial differential equations (For 
equations on manifolds, see 58Gxx; for manifolds of solutions, see 58Bxx; 
for stochastic PDEs, see also 60H15 ) -Partial differential equations with 



randomness (See also 60H15) ; 60H15 -Probability theory and stochastic 
processes (For additional applications, see llKxx, 62-XX, 90-XX, 92-XX, 
93-XX, 94-XX. For numerical results, see 65U05 ) -Stochastic analysis (See 
also 58G32) -Stochastic partial differential equations (See also 35R60) ; 
60H20 -Probability theory and stochastic processes (For additional 
applications, see llKxx, 62-XX, 90-XX, 92-XX, 93-XX, 94-XX. For numerical 
results, see 65U05) -Stochastic analysis (See also 58G32 ) -Stochastic 
integral equations 
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The authors investigate initial value problems for equations of the form 
•$ (\partial/\parfcia*l t)-u\sb \ epsil on (x, t>, \ omega ) +b\str ■ * - 

\epsilon (x, t, \omega) \cdot\nabla u\sb \epsilon (x, t, \omega) =F\sb 
\epsilon (x, t, \omega) $ $(t>0$, $x\in\bold R\sp k)$, $u\sb 
\epsilon (x, 0, \omega) =\phi\sb \epsilon (x, \omega) $ $ (x\in\bold R\sp k)$, 
where $\epsilon>0$ is a small parameter and $\omega$ denotes the randomness 
of the coefficients and consequently, of the solutions. Two cases are 
considered in detail: (1) If the randomness in the coefficients has order 
$0 (\epsilon) $ , then the perturbed problem is shown to have a unique 
solution $u\sb \epsilon$ for each small value of $\epsilon$. Furthermore, 
if $$\bold E\vert b\sb \epsilon (x, t ) -b\sb 0(x,t)\vert \sp p+\bold E\vert 
F\sb \epsilon (x, t ) -F\sb 0(x,t)\vert \sp p\leq c\epsilon\sp p$$ holds for 



some $p\geq 1$, then it follows that $\bold E\vert u\sb \epsilon (x, t ) -u\sb 
0(x,t)\vert \sp p\leq c\sb l\epsilon\sp p$, where, specifically in the 
case $\epsilon=0$ , all coefficients are nonrandom. (2) If the coefficients 
are highly oscillatory, in particular, if there exist mappings $b$ and $F$ 
being Lipschitz continuous with respect to $x$ and such that $b\sb 
\epsilon (x, t, \omega) =b (x, t/\epsilon, \omega) $ and $F\sb 
\epsilon (x, t, \omega) =F (x, t/\epsilon, \omega) $, the authors show that 
$\lim\sb {\epsilon\to 0}\bold E\vert u\sb \epsilon (x, t ) -u\sb 0(x,t)\vert 
=0$ uniformly on any finite time interval, providing the existence of a 
function $\overline b$ with $$\lim\sb { T\to\inf ty } \sup\sb {t\sb 0>0}\bold 
E\Big\vert T\sp { -1 } \cdot\int\sp {t\sb 0+T}\sb {t\sb 
0 }b (x, s) \, ds-\overline b (x) \Big\vert =0,$$ and similarly for $F$. In 
addition, under a mixing condition on the coefficients the authors 
approximate the distribution of the normalized solution $\epsilon\sp 
{-l/2}(u\sb \epsilon (x, t, \omega) -u\sb 0(x,t))$ in terms of a Gaussian 
Markov process that is found as the solution of an associated linear 
.stochastic &if f ere.ntia.1 equation. ..Finally,, the results, .are applied to 
approximate the moments of the solution to the case where $b\sb 0\equiv 0$, 
and to the case where the coefficients $b$ and $F$ have a special form of 
almost periodic random processes. 
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Abstract Length: MEDIUM (17 lines) 
Particle methods, sometimes referred to as vortex methods, are 
numerical methods that have been developed and used in computational 
physics, especially in computational fluid dynamics. The idea is to 
simulate a fluid flow by a finite number of fluid particles. The solution 
.can be .descxibecl ma,thejnaticaJ,ly by, a. sum pf Dirac. fielta functions .. In this 
way, the solution of equations governing the fluid flow can be obtained by 
tracking a finite number of fluid particles and evaluating • velocities by 
discretizing certain singular integrals. Similar ideas can be applied to 
solving equations arising in cell biology. 

Phase models are useful for studying synchronization of bacterial cell 
culture growth and other biological events associated with cell cycles. 
This paper considers a model that allows the growth rates of cells to 
change at different phases of the cell cycle. A particle method is derived 
for solving the weak formulation of this model, and convergence of this 
particle method is proved. 
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This book is a somewhat disorganized introduction to dynamical systems 
and their perturbations and bifurcations. Although a great many theorems 
are stated (almost always with no proof or at best a very sketchy proof) , 
these theorems are often hidden under obscure or misleading section titles. 
(This confusion of titles begins with the title of the book itself, the 
great majority of -which deals" with' nonchaotic systems 'Under the section 
heading Stroboscopic methods 11 in the chapter on Free oscillations'' 
one finds brief discussions of chaotic interval maps, entropy, Markov 
chains, circle maps, annulus maps, and homoclinic points, all of which are 
primarily of importance for forced rather than free oscillations.- Under 
''Implicit function theorems'' one finds the Fredholm alternative and the 
Lyapunov- Schmidt method. The preface claims that the book makes 
significant use of the concept of stability under persistent disturbances 
(SPD) . In fact, this concept is defined on page 91 (it is a somewhat odd 
combination of certain aspects of orbital and structural stability) , one 
theorem about it is stated without proof on page 104, and on page 191 it is 
claimed that the ideas of the method of averaging are related to those of 
SPD; these are the only references to SPD that the reviewer could find. 
Occasionally the statement of a theorem is misleading; the section on 
Floquet ' s theorem claims (as usual without proof) that the logarithm of a 
nonsingular matrix always exists, without pointing out that such a 
logarithm may be complex and that in these cases a real Floquet reduction 
can only be obtained by using double the least period of the system. 

The chapters of this book cover the following topics. (The descriptions 
are the reviewer's own, not those of the chapter and section titles.) 
Chapter 1,- "linear *s*ystrems. Chapter 2 , "introductory 'ideas about ridnlinear 
systems and maps. Chapter 3, stability of various kinds. Chapter 4, 
bifurcation of rest points and periodic solutions. This much constitutes 
Part I of the book; the rest, Part II, is titled ''Perturbation methods 11 . 
Chapter 5, regular perturbation theory. Chapter 6, an introduction to 
forced oscillations and resonance, with the standard example of Duffing' s 
equation, and a somewhat disconnected section about fractals, Newton's 
method, Julia sets, and fractal basin boundaries. (Recall that Part II of 
the book is supposed to be about perturbation methods; the text gives no 
clue that these fractal basin boundaries do not occur within the range of 
parameters for which perturbation methods are suitable.) Chapter 7, 
methods of averaging, including a brief reference to the KAM theorem but 
none to higher-order averaging or the Nekhoroshev theorem. Chapter 8, 



initial layer and boundary layer problems and relaxation oscillations. 
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From the preface: VN The techniques presented here range in mathematical 
difficulty up to calculus and matrix theory. The material is presented in 
general order of increasing mathematical difficulty. Some exercises deal 
with material in preceding sections, others are projects that extend 
preceding material . 

% "Our purpose in this book is not the systematic presentation of 
mathematical material, although there are important threads that run 
through several chapters . Instead, we hope to illustrate how mathematics 
can be used. In particular, our goal is to make available to students, 
having at least one term of calculus, topics in the. .life. . sciences ^ and 
medicine that have* benefited from mathematical modeling and analysis. In 
addition to exposing students to current ideas, the material is intended to 
reinforce their mathematics education by presenting familiar mathematical 
topics from novel points of view. Finally, enabling students to think in 
terms of models early in their academic experience should motivate them to 
develop and apply modeling skills further. 

""The mix of topics, taken largely from population biology and from 
physiology, includes important phenomena that are within reach of the 
students described above. The population part of the book draws its 
material from the areas of demographics, genetics, epidemics, and 
biogeography, while the physiological part surveys cardiovascular, 
pulmonary, renal, and muscle physiology. The final chapter is intended to 
introduce students to models of nerve cells and some neural circuits as a 
basis for studying how the brain works. 11 

The chapter headings are: Introduction; 1. The mathematics of 
populations: demographics; 2. Inheritance; 3. A theory of epidemics; 4. 
Biogeography; 5. The heart and circulation; 6. Gas exchange in the lungs; 

7. Control of cell volume and the electrical properties of cell membranes; 

8. The renal countercurrent mechanism; 9. Muscle mechanics; 10. Biological 
clocks and mechanisms of neural control. 
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A voltage-controlled oscillator neuron model (VCON) of a single neuron 
with cell body potential $\cos x(t)$, stimulated through a synapse having 
"presynaptic* potential '$\cos y(t)$,' is Modeled by 1 the' Equation * 
$dx/dt=\omega+ C\cos\sb +x(t)\cos\sb +y(t)$; where $\cos\sb +(u)$ denotes 
$\cos(u)$ if $\cos(u)\geq 0$ and 0 otherwise, and $C$ denotes connection 
strength. It is shown numerically that phase locking to input frequencies 
occurs in this model. A plot is made of the output/input frequency ratio 
$R=\lim\sb { t\rightarrow\infty} x(t)/y(t)$, in the case $y(t)=\mu t$ . A 
graph of $R$, as a function of $\omega/\mu$ , is shown to have visible 
plateaus at 1, $\tfrac 12$, etc. Learning is demonstrated by a second model 
$dx/dt=\omega+u$, $\tau\ , du/dt+u=C\cos\sb +x\cos\sb +y$, with time constant 
$\tau$. For this model, the plateau at $R=1$ becomes noticeably wider. In a 
third example, a network of VCONs is modeled by the system $dx/dt=w+f (x) $ . 
It is shown that if $w$^is nearly proportional to a vector of integers $n$, 
then $x\sb i/x\sb k$ can tend to $n\sb i/n\sb k$ as $t\rightarrow\inf ty$ , 
and again phase locking can occur. A circular array of VCONs, connected in 
rotationally symmetric fashion, is used to demonstrate self-organization in 
both the space and time domains. For this model the connection strengths 
are also modified. 
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Summary: Synchronization of cell doubling times due to alternating 
starvation-nutrition cycles is studied here using a method based on 
nonlinear Lexis diagrams and the assumption that the cell cycle has three 

phases pre-replication, replication and post-replication the middle of 

which is always of fixed length once started. 1 ' 
{For the entire collection see MR 90j: 92032}. 
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A mathematical analysis of small mammal populations . 
Hoppensteadt , F. C. (Department of Mathematics, Michigan State 
University, East Lansing, 48824, Michigan) 
Murphy, L. (Department of Mathematics, Oregon State University, 

Corvallis, 97331, Oregon 
Corporate Source Codes: 1-MIS; 1-ORS 
J. Math. Biol. 

Journal of Mathematical Biology, 1987, 25, no. 3, 263 — 274. ISSN: 
0303-6812 CODEN: JMBLAJ 
Language: English 
Document Type: Journal 
Journal Announcement: 1916 
Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (12 lines) 

Summary: v s Populations of Microtus montanus, the montane vole, have been 
extensively studied. It is known that their reproductive activity is 
closely linked to the availability of the chemicals in growing plants. We 
use a mathematical model here to study how the length of the vegetative 
season and the natural reproduction rhythm of voles are involved in the 
long term dynamics* 'of the population' ntimbers . In' particular,' *we ' use data 
obtained from Timpie Springs, Utah, and from Jackson Hole, Wyoming, to 
formulate a model. The novelty of this model is its use of littering 
curves that highlight the temporally discrete nature of vole reproduction. 
The model shows how the timing of the vegetative season can influence vole 
population sizes. ' ' 
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The book aims to be a textbook on neuronal modeling for an undergraduate 
course in mathematical biology. In fact it outlines mainly the studies of 
the author concerning the oscillatory aspects of the neuronal activity, by 
using his own model of the neuron and methods of electronic circuits. In 
the first part, the model of axonic propagation of Hodgkin and Huxley and 
its simplified form by FitzHugh and Nagumo are reviewed briefly. Later, a 
neuromime model is introduced and its characteristics of phase-resetting 
are analyzed. The neuromime is based on an integrated circuit which is 
able to simulate the activity of a neuron only for a variable related to 
the phase of firing. The properties of phase-blocking in coupled 
neuromimes' a're studied by using the 4 rotation ve6tbr* method . A section is 
devoted to the simulation of the control system of the respiratory 
activity in humans and to the reproduction of the activity-splitting 
patterns in rodents. The study of neural networks is undertaken in the 



last two sections, in which the model of Hopfield is briefly described and 
a continuous network, composed of neuromimes of the type considered by the 
author, is introduced. Here, the properties of phase-blocking in a 
discrete network,* -wit-h eight elements p are also 1 st'udi-e'd ' by means'of a ' 
numerical simulation. A short appendix on the mathematical background, 
with a model of the harmonic oscillator, and the Fourier and Laplace 
transforms, ends the book. The bibliography does not give complete 
information about the modeling of (rhythmic) activity either of single 
neurons or of neural networks. All the sections are supplied with 
exercises, the solutions for which are given. 

Reviewer: Ventriglia, F. (Naples) 

Review Type: Signed review 

Descriptors: *92A09 -Biology and behavioral sciences-Physiology, 
biochemistry (See also 76-XX, in particular 76Zxx, and 78A70, 80A30, 80A32, 
92A27, 92A4 0) 



1/5/250 (Item 40 from file: 239) 

DIALOG (R) File 239:Mathsci 

(c) 2004 American Mathematical Society. All rts . reserv. 
02000309 CIS 8801217 

A mathematical analysis of small mammal populations 
Hoppensteadt; C* * i * •■ *. . - - - 

Murphy, L. 

J. Math. Biol. JMathBio (CIS abbrev) 

Journal of Mathematical Biology, 1987, 25, 263-274 ISSN: 0303-6812 
CODEN: jmblaj 

Language: English 
Document Type: Journal 

Subfile: CIS (Current Index to Statistics) ASA/IMS 
Identifiers: Renewal theory 



1/5/251 (Item 41 from file: 239) 

DIALOG (R) File 239:Mathsci 

(c) 2004 American Mathematical Society. All rts. reserv. 
01978068 MR 87h#92024 

Synchronized oscillations in networks of neuron analogue circuits. 
Hoppensteadt, F. C. (Department of Mathematics, University of Utah, 
Salt Lake City, 84112, Utah 
Corporate Source Codes : 1-UT 

IMA J. Math. Appl-. Med. Biol. k • * * - 

IMA Journal of Mathematics Applied in Medicine and Biology, 1984, 1, 
no. 2, 135--148. ISSN: 0265-0746 
Language: English 
Document Type: Journal 
Journal Announcement: 1816 
Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: SHORT (6 lines) 

Summary: * N We introduce a simplified circuit analogue of a nerve cell. 
This circuit is based on modulation of a voltage-controlled oscillator by 
signals entering through a circuit analogue of a chemical synapse. We 
describe phase locking of frequency-encoded information, and we show that 
large networks can sustain stable spatial patterns of phase-locked 
behaviour. ' 1 

Reviewer: Summary 

Review Type: Abstract 

Descriptors: *92A09 -Biology and behavioral sciences-Physiology, 
biochemistry (See also 76-XX, in particular 76Zxx, and 78A70, 80A30, 80A32, 
92A27, 92A40) 
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Stable oscillations of weakly nonlinear Volterra integro -differential 

equations . . ■ - ■ - 

Hoppensteadt, F. C. (Department of Mathematics, University of Utah, 
Salt Lake City, 84112, Utah) 
Schiaffino, A 

Corporate Source Codes; 1-UT 
J. Reine Angew. Math. 

Journal fur die Reine und Angewandte Mathematik , 1984, 353, 1--13. 
ISSN: 0075-4102 CODEN: JRMAA8 
Language: English 
Document Type: Journal 
Journal Announcement: 1704 
Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: SHORT (10 lines) 

For a class of Volterra integro-di ff erential equations results about the 
existence and stability of periodic solutions are derived which are 
analogous to known results for ordinary differential equations. The 
equations are small (as measured by a parameter $\varepsilon$ ) nonlinear 
perturbations of linear ones. There are two cases. In the first the 
perturbation is periodic and the result is about forced oscillations. In 
the second the unperturbed linear problem describes uncoupled harmonic 
oscillators. -and , the result is about phase-locked* self-.s-ustained - 
oscillations. An equation describing an electronic circuit is discussed as 
an example. 

Reviewer: Diekmann, O. (Amsterdam) 
"Review Type: Signed review 

Descriptors: *45J05 -Integral equations-Integro-ordinary differential 
equations ; 70K20 -Mechanics of particles and systems (For relativistic 
mechanics, see 83A05 and 83C10; for statistical mechanics, see 82-XX)- 
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An algorithm for approximate solutions to weakly filtered synchronous 
control systems and nonlinear renewal processes . 

Hoppensteadt, F. C. (Department of Mathematics, University of Utah, 
.Salt Late City,, 84112, Utah . . - . 

Corporate Source Codes : 1-UT 

SIAM J. Appl. Math. 

SIAM Journal on Applied Mathematics, 1983, 43, no. 4, 834 — 843. 
ISSN: 0036-1399 CODEN: SMJMAP 
Language: English 
Document Type: Journal 
Journal Announcement: 1524 
Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (18 lines) 

The author considers singular perturbation problems for nonlinear 
Volterra integro-diff erential equations whose kernels have components 
near the Dirac delta functional. Although in general the equations do not 
reduce to equivalent systems of differential equations, the author finds a 
way to attack the problem by natural generalizations of methods used to 
study singularly perturbed initial and initial-boundary problems for 
ordinary and partial differential equations. Using earlier results 
concerning differential equations [Comm. Pure Appl. Math. 24 (1971), 
807 — 840; MR 44#5576; Arch. Rational Mech. Anal. 35 (1969), 284—298; MR 
40#1694] the author applies the method of matched asymptotic expansions to 
the. Volterra- equations directly. Ln a separate, section- -the case ®f - linear ■ 
convolution equations is studied in great detail. The problems considered 
are well motivated by two instructive examples, one from the study of 
weakly filtered phase-locked loop electrical circuits, the other from 



nonlinear renewal equations describing the population dynamics of a 
demographic process. 

Reviewer: Groh, Jurgen (Jena) 

Review Type: Signed review 

Descriptors: *45D05 -Integral equations-Volterra integral equations (See 
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An averaging method for Vol terra integral equations with applications to 
phase-locked feedback systems. 

Equadiff 82 (Wurzburg, 1982) 
Hoppensteadt , F. C. (Department of Mathematics, University of Utah, 
Salt Lake City, 84112, Utah 
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1983, 

Springer, Berlin-New York, ; 256--26S, , 

Series: Lecture Notes in Math., 1017, 

Language: English 

Document Type: Proceedings Paper 

Journal Announcement: 1605 

Subfile: MR (Mathematical Reviews) AMS 

Abs tract Length: SHORT (5 lines) 

An averaging scheme is justified for the Volterra integro-dif f erential 
equation $\dot x(t)-f (t, \varepsilon) +\varepsilon\int\sb 0\sp 
tk(t-s)F(x(s) , \varepsilon) \,ds$, where $x(t)\in{\bf R}\sp N$ and 
$\varepsilon>0$ is a small parameter. A number of practical examples 
containing phase-locked loops are considered. 

{For the entire collection see MR 84j: 00010}. 

Reviewer:' Bainov," D.* "(Sofia) 

Review Type: Signed review 
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An extrapolation method for the numerical solution of singular 
perturbation problems . 

Hoppensteadt , F. C. (Department of Mathematics, University of Utah, 
Salt Lake City, 84112, Utah) 
Miranker, W. L. (IBM Thomas J. Watson Research Center,, Yorktown Heights, 
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SIAM J. Sci. Statist. Comput . 

Society for Industrial and Applied Mathematics. Journal on Scientific 
and Statistical Computing, 1983, 4, no. 4, 612 — 625. ISSN: 0196-5204 
CODEN: SIJCD4 

Language: English 

Document Type: Journal 

Journal Announcement: 1605 

Subfile: MR (Mathematical Reviews) AMS 

Abstract Length: SHORT, 9 lines 

Authors 1 summary: N % We show how the form of the perturbation 
approximation for the solution of stiff systems of ordinary differential 



equations with an identifiable small parameter can be used to generate 
associated nonstiff or relaxed equations. Solutions of these relaxed 
equations are easily calculated, and appropriate combinations of these 
solutions furnish numerical approximations to the original stiff problem. 
Variations of this method are applied to two classes of initial value 
problems: those with highly oscillatory solutions and those with rapidly 
equilibrating solutions. 11 

Reviewer: Summary 

Review Type: Abstract 

Descriptors: *65L99 -Numerical analysis-Ordinary differential equations- 
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Document Type: Technical Report 
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0303-6812 CODEN: JMBLAJ 
Language: English 
Document Type: Journal 
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Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (17 lines) 

Under certain conditions, the FitzHugh - Nagumo model of nerve 
membranes can be averaged and rewritten as r= mu r(l - r) , theta =1, from 
which one sees that all solutions starting with r > 0, theta = theta 
{sub}0, have the same phase for all time. This is the content of the 
notion of radial isochron clock (RIC) . The authors study the application 
of a stimulus A at regular time intervals when the phase is determined by 
theta =1+ delta (t/T - l/T)f(A, theta), where f is the strength of a 
phase resetting change, which occurs when t is an integer multiple of T, 
and delta is the Dirac measure. The coupling between an external clock 
(phase theta) and an endogenous clock (phase psi) is described by theta = 
mu, psi = omega +F(psi - theta), where F is the coupling, mu is the 
external clock's free frequency and finally F(phi +2 pi)=F(phi). RIC 
systems can be shown to exhibit phase locking; this phenomenon is 
described by a system theta = omega + epsilon F (theta, epsilon) , theta 
and F being n-vectors. 

Reviewer: Haimovici, A. (Iasi) 

Review Type: Signed review 

Descripto-r-s : *92A09 -Biology and* othe*r natural * sciences , behavioral 
sciences-Physiology, biochemistry (See also 76-XX, in particular 762xx, and 
78A70, 80A30, 80A32, 92A27, 92A40) ; 34C15 -Ordinary differential equations 
-Qualitative theory (See also 58Fxx) -Nonlinear oscillations 
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Integrate-and-f ire models of nerve membrane response to oscillatory 
input . 
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SIAM J. Appl. Math. 

SIAM Journal on Applied Mathematics, 1981, 41, no. 3, 503 — 517. 
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Journal Announcement: 1408 
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Abstract Length: LONG (30 lines) 

In the integrate-and- f ire model for the activity of nerve cells 
it is assumed that the nerve membrane potential at the locus where nerve 
impulses are generated summarizes or integrates the inputs (applied 
currents or post-synaptic potentials). At each time when the integral 
reaches a threshold value, an impulse is produced and the potential reset 
to its resting value. The integration is not perfect in time: the 
membrane potential steadily undergoes a linear decay process with a 
constant rate. Mathematically the model is described by the differential 
equation dv/dt= - av+f (t) in connection with the threshold condition 
v(t{sup)+)=0 if v(t)=v {sub}T, where v is the membrane potential, f the 



input function, a the decay rate, and v{sub}T the threshold. The 
authors give a complete typology of what happens if the input is 
sinusoidal: f (t ) =S { sub} 0+S { sub }m cos (omega t+ phi). Namely, the behavior of 
the model is characterized in terms of a unique, parameter-dependent 
rotation number, describing the asymptotic relationship between impulse 
times and phases of forcing. Three types of behavior are associated with 
regions in parameter space having positive Lebesgue measure: either the 
sequence of impulses is periodic andjphase locked t to .the ingut ,( rational 
rotation number)', or the impulse times are ergodically distributed across 
the input phases (irrational rotation number) , or the impulse sequence is • 
finite. Other types of behavior for parameter sets of measure zero are 
discussed as well. Besides these analytical results, some numerical 

examples of typical cases are presented (phase densities and interspike 
time histograms) . Finally, a close relationship is revealed to another 
biological oscillator model studied by L. Glass and M. C. Mackey [J. 
Math. Biol. 7 (1979), 339 - 352; MR 83b:92014], 
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Mathematical aspects of physiology (Proc. Summer Sem. , Univ. Utah, Salt 
Lake City, Utah, 1980) 
Hoppensteadt , F . C . 
1981, 

Amer. Math. Soc, Providence, R.I.,; pp. 327 — 344,, 

Series: Lectures in Appl . Math., 19, 

Language: English 

Document Type: Proceedings Paper 

Journal Announcement: 1323 
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Abstract Length: SHORT (5 lines) 

Several tunnel diode oscillator circuits are described which 
model FitzHugh - Nagumo type neurons. The phase-locking behavior of these 
neural models is described in detail using Fourier methods, stroboscopic 
methods and rotation numbers. (For the entire collection see MR 

82f :92003. ) 

Reviewer: Scott, Alwyn C. (Los Alamos, N.M.) 
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The response of an individual to a drug is modeled in the case 
when free and bound sites equilibrate rapidly and also in the case when the 
drug binds permanently to the"' receptbr '"sites . Then 1 the" population 'dynamics 
of drug use is modeled by using an infectious disease model where 
nonusers are susceptibles and active drug users are infectives. Nonusers 
are assumed to become users through contact with active users. A 
threshold parameter is identified so that below the critical threshold 
value, drug use dies out. Above the threshold value, there is an epidemic 
of drug use. No drugs or populations to which the model could be applied 
are mentioned. 

Reviewer: Hethcote, H. W. (Iowa City, Iowa) 
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Computation by extrapolation of solutions of singular perturbation 
problems 

Analytical and numerical approaches to asymptotic problems in analysis 
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Authors* summary: x N We show how the asymptotic form of the 
solution of singular perturbation problems can be used to generate 
associated unperturbed or relaxed equations. Solutions of these relaxed 
equations are easily calculated, and appropriate combinations of them 
furnish numerical approximations to the original problem. Variations of 
this method are applied to two classes of initial problems: those with 
rapidly equilibrating solutions and those with highly oscillatory 
solutions. ,,M ** ' '(Tor* the eriti re collection see 1 '8 im765005 .* ) ' " 
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This book gives a survey at an advanced undergraduate or graduate 
level of some of the more prominent mathematical methods which have proved 
useful for studying biological systems (i.e., population and genetic 
phenomena). The first two chapters deal with population dynamics (i.e., 
maximal sustained yield, maximal current revenue, optimal discounted future 
revenue), renewal theory, and reproduction matrices (i.e., honest and 
dishonest matrices), using difference equations in models of total 
population and population age structure. The third chapter surveys models 
of random population events using Markov chains (i.e., the two-state, 
hypergeometric, Polya, Fisher-Wright, and Reed - Frost chains). The last 
two chapters look at some of the mathematical methods used to examine the 
qualitative behavior of more complicated difference equations covering 
perturbation methods and dispersal processes. 
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Some Mathematical Questions in Biology (Proc. 13th Sympos . Math. Biol., 
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Authors' summary: N ^ A model of respiration control is formulated 
and analyzed. The model describes the blood concentration of a single 
solute (CO{sub}2) and the velocity of blood in brainstem circulation. The 
systemic and pulmonary circulations, the heart and the lung are lumped 
into one device (called the pump-filter) and the brainstem is represented 
by a single device (called the sensor-controller) . The model exhibits a 
range of behavior from regular to chaotic breathing patterns, reminiscent 
of Cheyne-Stokes breathing, as parameters are changed. 1 1 
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Bacteria can grow in spatial patterns in response to diffusion of 
a needed nutrient. In this paper the authors present a model which 
describes the histidine concentration, the concentration of the growth 
medium's buffer and the size of the bacterial population. The 
mathematical- problem Involves a •■coupl-ed- system* o*f -three • - - 

reaction-diffusion equations in a one-dimensional cylindrical domain. The 
reaction functions in these equations are based on Michaelis - Menten 
kinetics. A heuristic argument and a numerical example are given. 
(For the entire collection see MR 82c:92004.) 
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Proceedings of the Twelfth Summer Seminar on Applied Mathematics held at 
the University of Utah, Salt Lake City, Utah, June 15 — 27, 1980. Edited by 
Frank C. Hoppensteadt. 

Contributors: Hoppensteadt, Frank C. 
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Contents: Frank C. Hoppensteadt, Foreword (pp. v - vi) ; Charles 
S. Peskin, Lectures on mathematical aspects of physiology (pp. 1 - 107); 
Richard Skalak, Blood rheology (pp. 109 - 139); Stephen Childress, Aspects 
of physiological fluid mechanics (pp. 141 - 163); Raul Mendez, Numerical 



study of incompressible flow in a region bounded by elastic walls (p. 
165); T. W. Secomb, Kinematics of close-packed red blood cells in shear 
flow (pp. 167 - 170); John L. Stephenson, Case studies in renal and 
epithelial physiology (pp. 171 - 212); J. E. Wood, A 

statistical-mechanical model of the molecular dynamics of striated muscle 
during mechanical transients (pp. 213 - 259); Dieter Schenzle, On 
neuroendocrine control of ovarian function (pp. 261 - 264); A. T. 
Winfree, Pec.uliari.ties, in the impulse response of pacemaker ..neurons (pp. 
265 - 279); John Rinzel, Models in neurobiology (pp. 281 - 297); James P. 
Keener, Chaotic cardiac dynamics (pp. 299 - 325); F. C. Hoppensteadt , 
Electrical models of neurons (pp. 327 - 344) ; Jonathan Bell, Threshold 
and conduction in prototype models of myelinated nerve axons: a 
preliminary report (pp. 345 - 347); David L. Barrow, Threshold for a 
reaction-diffusion equation related to nerve conduction (pp. 349 - 353) ; U. 
an der Heiden, M. C. Mackey and H. 0. Walther, Complex oscillations in a 
simple deterministic neuronal network (pp. 355 - 360) ; H. T. Banks, 
Parameter identification techniques for physiological control systems 
(pp. 361 - 383); Walter T. Kyner and Gary A. Rosenberg, Parameter 
estimation techniques used in the determination of the bulk flow of brain 
interstitial fluid (pp. 385 - 388); James Wiskin, Modelling of 
stimulation evoked acetylcholine release from myenteric neurons and 
estimation of the parameters (pp. 389 - 394). (The papers of 

mathematical interest that appear to be in final form are being reviewed 
individually. ) 
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Computer studies of nonlinear oscillators . 

Nonlinear oscillations in biology (Proc. Tenth Summer Sem. Appl . Math., 
Univ. Utah, Salt Lake City, Utah, 1978) 
Hoppensteadt , Frank C . 
1979, 

Amer. Math. Soc, Providence, R.I.,; pp. 131 — 139,, 
Series: Lectures in Appl. Math., 17, 
Language: English 
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Journal Announcement: 1215 

Subfile: MR (Mathematical Reviews) AMS 

Abstract Length: MEDIUM (25 lines) 

A sample of nonlinear oscillation problems is presented to 
illustrate several numerical methods which have been successfully applied 
to the study of oscillatory phenomena. Results are summarized in this 
survey article. First, studies of Mathieu's equation suggest computer 
experiments which can be performed to study the response of a pendulum to 
oscillation of its support point. A pendulum having a vertically 
oscillating support point can have three (at least) stable responses 
coexisting for the same parameter values of forcing and tuning; the 
straight up and straight down positions as well as a continual rotation 
are all stable solutions of this problem. Computer experiments can be 
used in these cases to determine the domains of attraction of the various 
modes of response. This example also demonstrates the interesting fact 
that otherwise unstable static states can be stabilized by external 
oscillatory forcing. Next, stable oscillatory responses of the van der 

Pol equation to periodic forcing are described. The concept of rotation 
number can be used to summarize numerical experiments on subharmonic 



solutions. Then it is described how computation of the power spectrum can 
be used to study the presence of higher harmonics. Finally, computer 

studies of certain chaotic dynamical systems are discussed. Numerical 
simulations are used to describe density functions characterizing random 
behavior of deterministic oscillators. (For the entire collection see 

MR 81a:92002. ) 1 
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Nonlinear oscillations in biology. 

Proceedings of the Tenth Summer Seminar on Applied Mathematics held at 
the University of Utah, Salt Lake City, Utah, June 12 — 23, 1978. Edited by 
Frank C. Hoppensteadt . 

Contributors: Hoppensteadt, Frank C. 

Publ: American Mathematical Society, ^Providence, % .B.. J . 
1'979, x+253 pp." ISBN: 0-8218-1117-7 
Series: Lectures in Applied Mathematics, 17. 
Price: $29.20. 
Language: English 

Document Type: Book; Proceedings 
Journal Announcement: 1215 

Nonlinear oscillations in biology; Summer Seminar: Applied Mathematics,; 
Biology; Salt Lake City, Utah, 10th 1978 
Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (12 lines) 

Contents: Frank C. Hoppensteadt, Foreword (pp. ix - x) ; Louis N. 
Howard, Nonlinear oscillations (pp. 1 - 67); Charles R. Steele, Studies of 
the ear (pp. 69 - 91); Arthur Winfree, 24 hard problems about the 
mathematics of 24 hour rhythms (pp. 93 - 126) ; Donald Ludwig, Stochastic 
modelling and nonlinear oscillations (pp. 127 - 129); Frank C. 
Hoppensteadt, Computer studies of nonlinear oscillators (pp. 131 - 139) ; 
O. E. Rossler, Chaotic oscillations: an example of hyperchaos (pp. 141 - 
156); Jack K. Hale, Nonlinear oscillations in equations with delays (pp. 
157 - 185); John Guckenheimer , A brief introduction to dynamical systems 
(pp. 187 - 253). t JMost of the papers are be|ng .reviewed indiyidually . ). 
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Threshold analysis of a drug use epidemic model 
Hoppensteadt, F. C. 

Murray, J. D. 
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Mathematical Biosciences. An International Journal, 1981, 53, 79- 
87 ISSN: 0025-5564 CODEN: mabiar 
Language: English 
Document Type: Journal 
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Slow selection analysis of genetic traits in synchronized populations. 

Conference on Deterministic Differential Equations and Stochastic 
Processes Models for Biological Systems (San Cristobal, N.M., 1977). 
Hoppensteadt , F. C. 
Rocky Mountain J. Math. 

The Rocky Mountain Journal of Mathematics, 1979, 9, no. 1, 93 — 97. 
ISSN: 0035-7596 CODEN: RMJMAE 
Language: English 
Document Type: Journal 
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Abstract Length: MEDIUM (15 lines) / 

The behaviour over time is considered for the frequencies of the gamete 
types for a genetic structure taken to be two autosomal loci each having 
two alleles in a diploid population. The population is taken to be 
synchronized in the sense that all reproductions occur at once followed by 
the removal of parents. Selection is assumed to act over a much longer time 
scale than reproduction. The first result is tha£ w.fyen .selection 4 s . slow 
and linkage right, Fisher's fundamental theorem of natural selection 
holds, namely the population's genetic structure changes in such a way as 
to increase its fitness. The second result is that the effects of slow 
selection and loose linkage on a population initially at linkage 
equilibrium are such that the fundamental theorem can only be verified 
after an initial transient period in which the gamete frequencies 
equilibrate to a certain functional form. 
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A nonlinear renewal equation with periodic and chaotic solutions. 

. Asymptotic methods and singular perturbations J Pr.pc ( ,$IAM-AMS Sympos . # 
Appl. Math'."/ New York/ 1976) *' 
Hoppensteadt, Frank C. 
1976, 

Amer. Math. Soc. , Providence, R.I.,; pp. 51 — 60. SIAM-AMS Proceedings, 
Vol. X,, 

Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (14 lines) 

Author's summary: ' % A nonlinear renewal equation which arises in several 
areas of mathematical population theory is studied by a combination of 
mathematical and numerical analysis. The model is characterized by two 
parameters: $m$, a measure of the population's viability and fertility and 
$\mu$, the (normalized) length of the population's reproductive window. 
Solutions are described for all values of these parameters, $0\leq m\leq 
4$, $0\leq\mu\leq 1$, by a combination of multi-time perturbation analysis 
and numerical solution. In various regions, the solutions are shown to be 



described by Burgers' equation, a Korteweg-de Vries equation and by a 
nonlinear difference equation. Numerical methods are used to investigate 
the remaining regions . 1 ' 

\{For the entire collection see MR 58\#15769.\} 
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Mathematical theories of populations: demographics, genetics and 
epidemics . 

Regional Conference Series in Applied Mathematics. 

Hoppens teadt , Frank 
Publ: Society for Industrial and Applied Mathematics, Philadelphia, Pa., 
1975, vii+72 pp. 
Price: $7.05. 
Language: English 
Document Type: Book 
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From the preface: % "The monograph begins with a study of population age 
structure. A basic model is derived first, and it reappears frequently 
throughout the remainder. Various extensions and modifications of the basic 
model are th"en appTied to several '^popuration phenomena*/ such cis stable age 
distributions, self-limiting effects and two-sex populations. The second 
part is" devoted to population genetics, and it contains a summary of some 
of the most successful applications of mathematics in the biological 
sciences. Attention is focused on the derivation and analysis of a model 
for a one-locus, two-allele trait in a large randomly mating population. 
Then extensions of the system are considered which account for more 
complicated social structure (assortative mating and migration) and for age 
structure. This part ends with a description of Fisher's model for the 
propagation of a gene in a spatially distributed population, and stable 
gene waves are shown to exist. A reason for the success of mathematical 
theories in genetics has been the abundance of precise collectable data. 
Unfortunately, this is not the case in the topics discussed in Parts I and 
III. The final part, Part III, is concerned with the dynamics of contagious 
phenomena in a population. These are studied in the context of epidemic 
diseases, but the same methods can be used to describe other phenomena such 
as rumors, fads and information as well as models for two interacting 
systems. Several classic examples are discussed first, then a general 
age-dependent theory is formulated. However, the emphasis in Part III is 
placed on studies of qualitative properties of several typical models. 
First, a threshold 'theorem is" deriVed for' an age- dependent epidemic; and 
then the long time behavior of solutions to a relapse-recovery model is 
determined. Finally, models for the spatial spread of contagion are derived 
and extensively discussed. 1 1 
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Frequency entrainment of a forced van der Pol oscillator. 
Flaherty, J. E. 
Hoppensteadt , F. C. 



Studies in Appl . Math. 
1978, 58, no. 1, 5 — 15. 
Language: English 
Document Type: Journal 
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This is an interesting study of a van der Pol relaxation oscillator 
subject to external sinusoidal forcing: $$ (d\sp 2u/dt\sp 2)+k(u\sp 
2-1) (du/dt) +u=\mu kB\cos (\mu t+\alpha) . $$ The model is studied for the 
parameters s*atisfyi'ng -$\alpha=0$, "$\mu^l$> $0<B<0 . '*$*0<l/k<0 . 2$ 7 The ' 
numerical computing of the rotation number $(\rho)$ suggests that it 
defines a continuous but piecewise constant surface except in overlap 
regions where it is double-valued, having what looks like folds. The 
parameter ranges where $\rho$ is single-valued illustrate the phenomenon of 
locking phase, the successive bifurcation of stable subharmonic and almost 
periodic oscillations. The paper gives an indication of the nature of 
stable responses of this system, which poses difficult analytic problems. 
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Optimal exploitation of a spatially distributed fishery. 

New trends in systems analysis (Proc. Internat. Sympos . , Versailles, 
1976) 

Hoppensteadt , Frank C. 

1977, 

Springer, Berlin,; pp. 3--18. Lecture Notes in Control and Informat. 
Sci., Vol. 2,, 

La-nguage: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (22 lines) 

Consider a fish population of density $u(x,t)$ at time $t$, where $x$ is 
the distance to that boundary of the habitat bordering a breeding ground. 
In this one-dimensional model, the second boundary $x=X$ is assumed to 
adjoin an unfavorable environment. The dynamical boundary condition at 
$x=0$ is of the form $\partial uApartial x=-A(u)=0$ if $\int\sb 0\sp 
Xu (x, t) \,dx<T$, and $\partial uApartial x=-A(u)=-A$ if $\int\sb 0\sp 
Xu (x, t ) \, dx\geq T$, where $T$ is a threshold parameter. A zero flux 
boundary condition prevails at the hostile boundary. 

If a bang^bang prirrciple motivates fhe* harvest' effort, then a free 
boundary (stop-harvest) of the form $x=S(t)$ appears on either side of 
which distinct diffusion relations hold, and on which a zero revenue 
relation holds. This is, of course, reminiscent of a two-phase Stefan 
problem of heat conduction. 

The author discusses in detail equilibrium solutions of this dynamical 
problem and relates this analysis to various important qualitative 
questions such as harvesting quotas, competition and ultimate fishery / 
collapse. In particular, threshold parameters are derived for the onset of 
collapse. A number of other questions are also discussed in this brief but 
interesting paper. 

\{For the entire collection see MR 57\#5216.\) 
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Dynamics of the Josephson junction. 
Levi, M. 
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Quart. Appl. Math. 
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Language: English 
Document Type: Journal 
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Abstract Length: SHORT (8 lines) 

The jump in the electron wave function across the junction gap is studied 
by considering the sine-Gordon equation. This analysis is performed for the 
current-driven case. Replacement of the derivatives by divided differences 
gives the discrete version of this problem. A numerical analysis of the 
static solutions, including their stability, is described, and the 
existence of a running periodic solution for both the discrete and 
continuous cases is proved. 

Reviewer: Sips, Vladimir (Zagreb) 
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Periodic solutions of a logistic difference equation. 
Hoppensteadt, F, C. 

Hyman, J. M. 

SIAM J. Appl. Math. 

1977, 32, no. 1, 73—81. 

Language: English 

Document Type: Journal 
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Consider a logistic difference equation $(*)$ $x\sb {n+l}=mx\sb n(l-x\sb 
n)$, where $m$ is a nonnegative real number. Since the authors are 
interested in the solution of $(*)$ that lies in the unit interval, they 
examine equation $(*)$ only for values of $m$ satisfying $0\leq m\leq 4$. 
Under this condition, by employing simple geometric arguments and numerical 
computations-, they -study periodic •■solutions of equation* "$(*)$ They show 
that as $m$ increases from zero, solutions having successively higher 
periods branch from old ones until the value $m\doteq 3.57$ is reached, 
after which there is an infinity of periodic solutions. The authors also 
make an interesting study of the behavior of the solution of $(*)$ in the 
chaotic regime. In fact, it is shown how, as $m$ increases from 3.57, 
solutions having various other periods are added to the solution set until, 
at $m\doteq 3.83$, solutions of period three, and so solutions of all 
periods, are present. Furthermore, in order to study the dynamics of 
solutions in portions of the chaotic regime, the authors give numerical 
calculations for the density functions. In the present work, the authors 
have used an interesting result due to T. Li and J. Yorke [Amer. Math. 
Monthly 82 (1975), no. 10, 985—992; MR 52\#5898] . 
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An analysis of transient behavior in the onset of convection. 

Singular perturbations and boundary layer theory (Proc. Conf., Ecole 
Centrale, Lyon, 1976) 
Gordon, Noam 

Hoppensteadt , Frank C. 
1977, 

Springer, Berlin,; pp. 231--243. Lecture Notes in Math., Vol. 594,, 
Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
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The bifurcation of convective motion for a viscous, incompressible fluid 
contained between horizontal parallel plates is studied using matched 
asymptotic expansions for an associated ordinary differential equation in 
an abstract space. The perturbation series is rigorously justified on the 
basis of previous work of the authors [Comm. Pure Appl . Math. 28 (1975), 
no. 3, 355 — 373; MR 54\#5927] . The series obtained is similar to those 
derived by earlier- -workers and the novelty of this* work "lies in its* 
rigorous justification. Only the cases of % 'rolls' 1 or a singly periodic 
disturbance are dealt with, although the authors state that the other cases 
may be dealt with analogously. 

\{For the entire collection see MR 56\#810.\) 
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Iterated averaging methods for systems of ordinary differential 
equations with a small parameter. 
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The method of averaging in its simplest form applies to a vector 
differential equation $\dot x=\varepsilon f (x, t , \varepsilon) $ (where 
$\varepsilon$ is a small parameter and $f$ is periodic in $t$) and consists 
in replacing $f$ by its average over $t$ and asserting that the resulting 
solutions are a good approximation to the exact solutions. The first-order 
average is not always an adequate approximation, especially when some 
components of the average of $f$ vanish, and one may construct higher-order 
averaged equations whose solutions give better approximations -[see L. 
Perko, SIAM J. Appl. Math. 17 (1968), 698--724; MR 41\#2129] . The present 
authors establish asymptotic validity of what is essentially a second-order 



averaging method under various circumstances, sometimes on expanding 
intervals of time and sometimes on infinite intervals. The setting is $\dot 
w=\varepsilon E (w, x, y , z , t, \varepsilon) $ , $\dot x=\varepsilon 
F(w,x,y,z,t, \ va reps i Ion) $ , $\dot y=\varepsilon G (w, x, y, z, t , \varepsilon) $ , 
$\dot z=A( t) z+\varepsilon H(w,x,y,z,t, \ va reps i Ion) $ , where $w,x,y,z$ are 
vectors, $A(t)$ is a matrix such that all solutions of $\dot z=A{t) z$ 
approach zero exponentially, and $E,F,G$ are quasiperiodic in $t$ and 
decompose into sums of periodic functions. It is explicitly assumed that 
the time average of $E$ vanishes; a second-order average is calculated for 
the $w$ components only, which vary pn,J:he extray.s^ow v tame sc^le., . 
$\varepsilon\sp 2t$, the first order sufficing otherwise. Applications are 
given, including an approximation to transient solutions in a Hopf 
bifurcation problem. Simpler proofs (along the lines of Perko) could 
probably be given for some of the results of this paper, at least for 
expanding intervals of length $0 ( 1/ \varepsilon) $ ; however results are 
obtained here on intervals of length $0 ( l/\varepsilon\sp 2)$, and the 
proofs are designed to be applicable in a partial differential equations 
setting as well. 
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Multitime methods for systems of difference equations. 
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Studies in Appl . Math. 
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The authors develop a constructive perturbation scheme which reduces the 
problem of determining behavior of solutions of a perturbed difference 
system to that for an associated system of ordinary differential equations 
via an averaging procedure. The general system treated has the form 
$x (m+1 ) =Ax (m) +\varepsilon f (x (m) , \ varepsilon) $ . Assumptions are that $A$ is 
similar to a block diagonal matrix containing one oscillatory block and one 
stable block, that the part of the solution corresponding to the stable 
block has a smooth solution, and that the appropriate average exists 
smoothly. 

Particular results are obtained for the linear one-and two-parameter 
systems $x (m+1 ) = (A-f\ varepsilon B+\mu 4 C)^ (m) $ . The. me th.od , developed is 
applied to examples including a population genetics model, a training 
algorithm in pattern recognition, and the numerical analysis of a stiff 
differential equation. Results obtained are similar to those obtained by 
the authors for differential equations [J. Differential Equations 22 
(1976), no. 2, 237—249; MR 54\#10777] . 
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Mathematical methods of population biology. 



Hoppensteadt, F. C. 

Publ: Courant Institute of Mathematical Sciences, New York University, 
New York, 

1977, v+167 pp. 
Price: $5.25. 
Language: English 
Document Type : Book 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: LONG (33 lines) 

This 'is a survey of 'problems and* mathematical methods* "of population 
biology. The choice of topics and emphasis is, in this reviewer's opinion, 
rather idiosyncratic. If used as a text, these notes would require 
substantial amplification in places, and references are also scanty. 
However, mathematically inclined readers will find a wealth of ideas and 
insights . 

Chapter I .begins with a discussion of the dynamics of a single 
population, with emphasis on linear models. This theme is taken up again in 
Chapter V, where age structured models are presented, including Leslie 
matrices and stable age distributions. Chapter I closes with an 
introduction to the theory of exploitation of populations . 

Chapter II is entitled Mathematical Ecology. It is a brief account of 
the periodic cicada, and of the spruce budworm. This must be the first work 
on population-ecology which does not mention the Lotka-Volterra equations, 
not even to point out that they are structurally unstable. In fact, there 
is no discussion of interacting species at all. 

Chapter III is concerned with contagion, i.e., the spread of infectious 
disease. It discusses the Reed-Frost and Kermack-McKendrick models, and 
presents some threshold results. Genetic models are treated in Chapter IV, 
VI and VII. These are all of the one or two locus variety. The simplest 
case of bacterial p'las'mi'd inheritance is dealt with' ' Ii r st . Then 
deterministic theories for one locus with two or three alleles are 
presented, followed by two loci, with two alleles each. Finite populations, 
where random sampling effects appear, are dealt with in Chapter VII. 
Emphasis is placed upon branching process approximations, and diffusion 
approximations . 

The final chapter is devoted to spatial effects. Linear and nonlinear 
diffusion models are presented. Topics include critical patch size for 
phytoplankton, diffusive instability, travelling waves, and genetic clines . 
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As an example of a more general type of nonlinear diffusion equation, the 
authors discuss in detail the system $\partial uApartial t=\theta\sb 
1 (\partial\sp 2u/\partial x\sp 2 ) +\alpha ( \lambda ) u-\beta v-\lambda u\sp 
3$, $\partial v/\partial t=\theta\sb 2 ( \partial\sp 2v/\partial x\sp 
2)+\beta u+\ gamma ( \ lambda ) v-\mu\ lambda u\sp 2v$, where $u$ and $v$ are 
functions of $x$ and $t$, $\theta\sb 1$ and $\theta\sb 2$ are diffusion 
coefficients, $\beta$ and $\mu$ are positive constants and $\lambda$ is a 
parameter. For $\lambda=\lambda\sb 0$, $\alpha ( \lambda\sb 



0) =0=\gamma (\lambda\sb 0)$ and neglecting diffusion the system has a 

stable equilibrium point at $u=0$, $v=0$ . It is assumed that, for $\lambda$ 
close to $\lambda\sb 0$, $\alpha ( \lambda ) =\varepsilon\sp 2\alpha\sb 
2+0 (\varepsilon\sp 3)$ and $\gamma ( \lambda) =\varepsilon\sp 2\gamma\sb 
2+0 (\varepsilon\sp 3)$, where $\varepsilon=\lambda-\lambda\sb 0$ and 
$\alpha\sb 2,\gamma\sb 2$ are positive constants. This implies 

that neglecting diffusion the equilibrium point is unstable for 

$\varepsilon$ non-zero. 

Under these circumstances, the authors look for a finite amplitude 
solution for* the f*ull -equation, including- diffusion'/* of* 'the form* * • ' 
$u=R(\eta) \cos\{\beta t\sp \ast+\phi (\eta) \ } $ , $v=R (\eta) \sin\ { \beta t\sp 
\ast+\phi (\nu) \ } $, where $\eta=\varepsilon (x-\varepsilon ct)$ and $t\sp 
\ast/t=l+0 (\varepsilon) $, with $c$ constant. Thus $R$ gives a slowly 
varying amplitude and $\phi$ a slowly varying phase angle for the 
sinusoidal oscillation. The authors show that a solution of this kind is 
possible and indicate that similar solutions can occur in the more general 
case . 

Reviewer: Brown, Archibald (Canberra) 
Review Type: Signed review 
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Singular perturbations and boundary layer theory. 

Proceedings of the Conference held at the Ecole Centrale de Lyon, Lyon, 
December 8 — 10, 1976. Edited by C. M. Brauner, B. Gay and J. Mathieu. 
Lecture Notes in Mathematics, Vol. 594. 

Contributors: Brauner, C. M. ; Gay, B. ; Mathieu, Jean; Hoppensteadt , 
Frank C. ; Miranker, Willard L. 

Publ: Springer-Verlag, Berlin-New York, 

1977, vii+539 pp. 

Language: English 

Document Type: Book 

Singular perturbations and boundary layer theory; Conference: Singular 
Perturbations and Boundary Layer Theory; Boundary layer theory; Lyon, 
(Proc. Conf., Ecole Centrale, Lyon, 1976) (Ecole Centrale, Lyon, 1976) 1976 

Subfile: MR (Mathematical Reviews) AMS 

Abstract Length: LONG (63 lines) 

Table of Contents: Preface (pp. iii-iv) ; J. Baranger, Estimations 
d'erreur a l'interieur pour un probleme de couche limite (pp. 1--9) ; A. 
Bensoussan, J. L. Lions and G. Papanicolaou, Perturbations et 

augmentation ' 1 des conditions initiales (pp. 10 — 29); Yu. A. Berezin [Ju. 
A. Berezin]*,- G . - 1 Dudnikova/' V. A. Novikov and Nv N*. '-Ya'nenko • [N . : 
Janenko] , Analytical and numerical studies of equations with sign changing 
viscosity coefficient (pp. 30--38) ; A. Bourgeat and R. Tapiero, Calcul de 
correcteurs pour un probleme de couche limite provenant de la physique du 
plasma (pp. 39 — 49); C. M. Brauner and B. Nicolaenko, Singular 
perturbation, multiple solutions and hysteresis in a nonlinear problem (pp. 
50--76) ; J. S. Darrozes, Comportement singulier des ecoulements a grand 
nombre de Reynolds au voisinage du bord d ! attaque d'une aile mince (pp. 
77--107); A. S. Demidov, Sur les problemes elliptiques 

pseudo-dif f erenciels , a petit parametre dans l'operateur principal (pp. 
108--122) ; A. S. Demidov, Sur la perturbation vx singuliere 1 1 dans un 
probleme a frontiere libre (pp. 123--130) ; G. Duvaut, Comportement 
macroscopique d'une plaque perforee periodiquement (pp. 131--145) ; Wiktor 
Eckhaus, Matching principles and composite expansions (pp. 146 — 177); Ph. 
Gatignol, Theorie asymptotique des trains d'ondes lentement modules pour 
certaines classes d ! equations de conservation (pp. 178--200) ; J. Genet and 
M. Madaune, Perturbations singulieres pour une classe de problemes 
hyperboliques non lineaires (pp. 201--230) ; Noam Gordon and Frank C. 



Hoppensteadt, An analysis of transient behavior in the onset of convection 
(pp. 231 — 243) ; Jean-Pierre Guiraud, The dynamics of rolled vortex sheets 
tightly wound around slender ..vortex filaments iniinvisoid incompressible 
flow (French summary) (pp. 244 — 259); Maurice Holt and Mohammed Y. 
Hussaini, Laminar wakes in supersonic flow (pp. 260 — 274); Frank C. 
Hoppensteadt and Willard L. Miranker, Applications of multi-time methods to 
pattern recognition and other problems (pp. 275 — 287) ; Denise Huet, 
Perturbations singulieres de problemes elliptiques (pp. 288 — 300); Th. 
Levy, Equations et conditions d 1 interface pour des phenomenes acoustiques 
dans des milieux poreux (pp. 301 — 311); Pierre Louvet and Jean Durivault, 
Compressible countercurrent flow in a strongly rotating cylinder pp. 
312--333) ; John L. Lumley, Some applications of singular perturbations to 
problems in fluid mechanics (pp. 334 — 350); J. Mathieu and B. Gay, Quelques 
problemes asymptotiques en theorie de la couche limite (pp. 351 — 364); F. 
Mignot and J. P. Puel, Un resultat de perturbations singulieres dans les 
inequations variationnelles (passage du 2$\sp { \ text { V erne }} $ ordre au 
l$\sp {\text{er}}$ ordre) (pp. 365 — 399); F. Obermeier, The application of - 
singular perturbation methods to aerodynamic sound generation (pp. 
400 — 421); R. E. O'Malley, Jr. and J. E. Flaherty, Singular 
singular-perturbation problems (pp. 422 — 436) ;.E. Sanchez-Palencia, 
Perturbations spectrales liees a la vibration d'un corps elastique dans 
l ! air (pp. 437 — 455); D. Brian Spalding, Numerical computation, of steady 
boundary layers, (pp. 456 — 473); Luc Ta&tar, Homogeneisa-tion en ■ <• - 
hydrodynamique (pp. 474 — 481); V. A. Trenogin, On the theory of abstract 
elliptical equations containing a singular perturbation (pp. 482 — 490) ; V. 
A. Trenogin and N. A. Sidorov, Regularisation of computation of branching 
solutions of nonlinear equations (pp. 491 — 505); Milton Van Dyke, From zero 
to infinite Reynolds number by computer extension of Stokes series (pp. 
506--517); R. Kh. Zeytounian, Les modeles de couche limite en dynamique de 
l'atmosphere (pp. 518--539) . 

\{The papers of mathematical interest will be reviewed individually . \ } 
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Numerical- roe thods for stiff systems of differential -equations related - 
with transistors, tunnel diodes, etc. 

Computing methods in applied sciences and engineering (Proc. Internat. 
Sympos., Versailles, 1973), Part 1 
Miranker, Willard L. 

Hoppens teadt , Frank 
1974, 

Springer, Berlin,; pp. 416 — 432. Lecture Notes in Comput . Sci., Vol. 
10,, 
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Subfile: MR (Mathematical Reviews) AMS 
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The authors consider systems of ordinary differential equations whose 
solutions are composed of slowly varying components, highly damped 
components, highly oscillatory components and combinations of these. Such 
stiff systems arise, for example, in models of circuits containing 
transistors or tunnel diodes. A relationship between stiff differential 
equations and differential equations with singular perturbations,, 
discovered by Miranker, is used to develop numerical schemes for the 
solution of. .stiff .equations...- .» , .. t * ■ . 

First, the multitime technique of asymptotic expansions is combined with 
the method of averaging of Bogoljubov to produce a procedure for deriving 
the asymptotic form of solutions to singularly perturbed differential 
equations over the full indicated range of solution behavior. A new concept 



is then introduced for the numerical solution of differential equations. A 
quantity is accepted as an approximation to the solution at a point in time 
of the quantity t approximates ..any yalue.^which the , splut.lon assumes, on a , 
neighborhood if that point in time, the size of the neighborhood being 
arbitrary but positive. With this new concept of a numerical solution, a 
numerical method is constructed based on the asymptotic theory developed 
earlier. 

\{For the entire collection see MR 49\#4197.\} 
Reviewer: Katz, I. Norman 
Review Type: Signed review 
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Differential equations having rapidly, changing t spluti-ons : analytic 
methods for weakly nonlinear systems . 
Hoppensteadt , F . C . 

Miranker, Willard L. 

J. Differential Equations 

1976, 22, no. 2, 237 — 249. 

Language: English 
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Abstract Length: MEDIUM (13 lines) ' 

The authors derive an approximation to the solution of the initial value 
problem (1) $dx/d\tau=\varepsilon F ( \tau, t , x, y , \ varepsilon ) $ , 
$dy/d\tau=A\ { \tau) y+\ varepsilon G (\tau, t, x, y, \ varepsilon) $, 

$x (0) =\xi (\varepsilon) $, $y ( 0 ) =\eta ( \varepsilon) $ , where $\varepsilon$ is a 
small parameter $t=\varepsilon\tau, x, \xi$, $F\in E\sp m$, $y,\eta$, $G\in 
E\sp n$ and $A\in E\sp {n\times n}$. The solution of (1) is considered in 
the form $x=X\sb 0 ( \tau, t ) +\varepsilon X\sb 1 (\tau, t, \varepsilon) $, $y=Y\sb 
0 (\tau, t) +\varepsilon Y\sb 1 ( \tau, t , \varepsilon) $ , and it is shown that 
under suitable assumptions $X\sb 1$ and $Y\sb 1$ are smooth functions that 
remain bounded uniformly in $\varepsilon$ for $\tau\in [ 0 , T/\varepsilon] $ , 
.$t\in[0,T]$, ..... . 

Reviewer: Rab, M. 
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Analysis of some problems having matched asymptotic expansion solutions . 
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The author, discusses^ the capabilities^ and limi^tipos^of th§ method of 
matched asymptotic expansions for the solution of evolution equations. The 
method is well known to be a highly useful heuristic in the qualitative 
analysis of physical problems. The author presents a number of cases in 
which the method has been shown to be mathematically rigorous an area to 



which he has been a major contributor. The results are listed according to 
spectral properties of the linear part of the problem resulting from 
scaling near a known steady state. The author shows that, when the linear 
problem is- stable,- "the solution can be'*written as' a' *f ini'te sum of "terms; 
each responding on a different time scale. He also shows that, when the 
linear problem is unstable, the method can be used to determine initial 
data that excite only decaying modes, and, in the case of bifurcation of 
new steady states, to construct the new states as well as the transients to 
them. The bibliography lists thirty-five useful references. 

Reviewer: Greenlee, W. M. 
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Nonlinear stability analysis of static states which arise through 
bifurcation . 

Hoppensteadt , Frank 

Gordon, Noam 

Comm. Pure Appl . Math. 

1975, 28, no. 3, 355—373. 
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The authors consider the initial value problem 
$du/dt=Bu+F (u, \varepsilon) $, $u ( 0) =u ( \varepsilon) $, where $u$ is an element 
of some Banach space $E$, the linear operator $B$ is a densely defined 
(possibly unbounded) operator acting in $E$, and $F$ is some nonlinear 
operator. It is known that if the spectrum of $B$ lies in the stable 
half-plane, then the solution of the system can be written as the sum of a 
solution which is a smooth function of $\varepsilon$ at $\varepsilon=0$ 
(called the outer solution) and a function depending on $t$ and 
$\varepsilon$ (called the initial correction) which decays exponentially as 
•$t\rightarrrfw\inf ty$ . * ■ .. . , . 

In this paper the authors study the above nonlinear problems when $B$ 
has a finite dimensional null-space. They use an extension of the method of 
matched asymptotic expansions, derived in the second author's dissertation, 
and show that the solutions of the problem can be written as the sum of an 
outer solution and an initial correction. But there may be several outer 
solutions and the one pertinent to the above nonlinear problem is 
determined by the position of the initial data relative to the domains of 
attraction of these various choices'. 
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Solutions* Heat: bifurcated Steady stattes'. 
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The author uses an example to illustrate the application of two-timing 
methods to bifurcation near a steady state. He begins with the 
two-dimensional system $\dot x=T (x, y, \ lambda )$ , $\dot y=G (x, y, \lambda ) $ , 
and linearizes about the steady state $(\lambda\sb 0,x\sb 0,y\sb 0)$. He 
assumes that the coefficient matrix of the linear part is of the form 
$\text{diag} (0, a) $, so that the new system is $\dot u=F (u, v, \varepsilon ) $ , 
$\dot v-av+G (u, v, \varepsilon) $ , with steady state $(0,0,0)$. Nearby (in 
$\varepsilon$ ) steady states will satisfy the condition 

$F (u, v, \varepsilon) =0$, $av+G (u, v, \varepsilon ) =0$ , and the latter yields 
$v=v\sp \ast ( u, \ varepsilon) $ . This is substituted into the first equation, 
and- this equation d's s*olved by a symptofic procedure's*. ''The author ' discusses 
the details when the initial states $u ( 0 , \varepsilon) $ and 
$v(0, \varepsilon) $ are restricted. The result is an existence and 
uniqueness theorem, with asymptotic estimates. 
\{For the entire collection see MR 51\#966.\} 
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A slow selection analysis of two locus , two allele traits . 
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The author derives a deterministic, continuous-time model for the 
dynamics of two-locus, two-allele Mendelian traits in a large randomly 
mating diploid population. The model allows for frequency and time 
dependent birth and death rates. It is analyzed under the assumption that 
the selective forces acting in the population are small. Slow selection 
approximations to the system's solutions are then constructed. Two 
particular cases are studied. First, when linkage between Loci is tight, 
the population is shown to rapidly approach Hardy-Weinberg proportions, 
which then may vary on a (slow) time scale determined by differential 
fitness. In the case of constant birth and death rates, a measure of the 
population's fitness is shown to increase on the slow time scale after an 
initial rapid adjustment. The second case considered is for loose linkage; 
a population near linkage equilibrium is studied. It is shown that the 
epistatic parameters cancel and that the results agree with the tight 
linkage case up to the leading order. The linkage disequilibrium is 
described in both cases. 

Reviewer :' "Wilsony Susan * i . . . 
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Analysis of a stable polymorphism arising in a selection-migration model 
in population genetics . 
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Author's summary: % 4 We study the dynamic behaviour of a population in 
which there is a geographic variation in fitness of a one-locus, two-allele 
trait. The results derived here complement certain results by W. H. Fleming 
[J. Math. Biol. 2 (1975), no. 3, 219--233; MR 53\#7531 who described the 
existence of a stable polymorphism when the ratio of selection intensity to 
dispersion rate exceeds a critical value. Here an approximation to the 
polymorphic -state d-s obtained, and 1 the •'•evolution 1 of 'oth^r initial " states to 
it is studied. 1 1 

Reviewer: Bruter, C. P. 
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Thresholds for deterministic epidemics . 

Mathematical problems in biology (Conf., Univ. Victoria, Victoria, B.C., 
1973) 

Hoppensteadt , Frank 
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The author considers a model for a deterministic epidemic and obtains an 
equation for the limiting number of exposed susceptibles , involving the 
number $\gamma$ of susceptibles exposed to each infective. Not all 
assumptions are stated explicitly (e.g., the exponential life time for 
infectives), and the chart in fig. 1 does not seem to agree with the claim 
that the threshold value for $\gamma$ is 1. 

\{For the entire collection see MR 49\#10392.\} 

Reviewer: Fischler, R. 

Review Type: Signed review 

Descriptors: *92A15 -Biology and other natural sciences, behavioral 
sciences-Population dynamics, epidemiology 



1/5/292 (Item 82 from file: 239) 

DIALOG (R) File 239:Mathsci 

(c) 2004 American Mathematical Society. All rts. reserv. 
01385286 MR 52##6127 

A geometric approach to boundary value problems for nonlinear ordinary 
differential equations with a small parameter. 

Analytic theory of differential equations (Proc. Conf., Western Michigan 
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The author considers singularly perturbed boundary value problems for the 
vector system $\varepsilon\ , dy/dt=f ( t , y, \varepsilon) $ when $f(t,0,0)=0$ and 
$f\sb y(t,0,0)$ is conditionally stable. Appropriate boundary value 
problems are shown to have the trivial limiting solution with endpoint 
boundary layers determined by separate initial value problems with data 
asymptotically restricted to lower dimensional manifolds. Significant 
extensions have been given by the author [Comm. Pure Appl . Math. 24 (1971), 
807 — 840; MR 44\#5576] and by others. Related work appears in the book by 
A. B. Vasileva and V. F. Butuzov [Asymptotic expansions of the solutions of 
singularly perturbed equations (Russian), Izdat. ^Nauka' 1 , Moscow, 1973]. 
\{For the entire collection see MR 51\#6002.\} 
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Consider the nonlinear integral equation (E) 
$x(t)=\alpha(t)+\gamma\int\sb { \max (-1 , t-1 ) } \sp tx(l-x) ( \eta ) \ , d\eta\ 
(\gamma>0, t\geq-;l)$ t , with the conditio^ $0\leq x /t t ).\leq ,1$ . For $t\geq 0$, 
(E)' reduces to ( E$\sb '{ \text+} $ ) $x ( t ) =\gamma\int\sb {"t-l}\sp" 
tx ( 1-x) (\eta) \, d\eta$ . In order to estimate the rates at which the 
solution approaches equilibrium, consider moderate or large values of 
$\gamma$: (a) for $0<\gamma\leq 1$, $x=0$ is the stable equilibrium point 
of (E$\sb {\text+}$) and $0<x(t)<Y\sb n$, $n-l\leq t<n$, where $Y\sb 
l=\gamma/4$ and $Y\sb n=\gamma Y\sb {n-1} (1-Y\sb {n-l})\ (n=2, 3, \cdots ) $ , 
i.e., the convergence $x ( \cdot ) \rightarrow 0$ is majorized by the 
convergence to zero of iterates of the polynomial $\gamma x(l-x)$; (b) 
define $w$ by $x ( t ) = ( \gamma-l ) Agamma+w ( t) $ . For $\gamma>l$, the stable 
solution is $x= ( \gamma-l ) /\gamma$ , i.e., the convergence 

$w(\cdot) \rightarrow 0$ is majorized by the convergence to zero of iterates 
of $w (2-\gamma-\gamma w)$. For $\gamma\gg 1$, the solution has a linear and 
an oscillatory part and depends on time scales other than $t$. A multitime 
perturbation expansion is derived. 

Equation (E) describes the behavior of a deterministic epidemic, where 
$x(t)$ is the normalized number of infectives at time $t$ and $\alpha(t)$ 
is the number of initial infectives remaining at time $t$ (the function 
$\alpha (\cdot ) $ is continuous on $[-l,0)$ and satisfies the condition 
$0<\alpha (t)<l$, $d\alpha/dt\leq 0$ for $-l\leq t<0$, and $\alpha\equiv 0$ 
for $t>0$) ..The t existence and uniqueness of solutions v o„f . ( E$\sb {\text + ]$) 
have been proved by the second author and P. Waltman [Math. Biosci. 12 
(1971), 133—145; MR 46\#1373] . 
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Advances in differential and integral equations. 

A collection of papers presented at the Conference on Qualitative Theory 
of Nonlinear Differential and Integral Equations held at the University of 
Wisconsin, Madison, Wis., August 12--23, 1968. In memoriam: Rudolph E. 
Langer (March 8, 1894 to March 11, 1968). Edited by John A. Nohel . 

Contributors: Langer, Rudolph E.; Nohel, John A.; Grafton, R. B. ; Kaniel, 
S.; Berger, M. S . ; Newell, Alan C; Brunovsky, P.; Winston, Elliot; 
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From the editor's preface: v % This volume consists of invited lectures, 
together with abstracts of contributed papers, delivered at the conference. 
As indicated by the table of contents, the lectures divide naturally into 
five areas; some of the papers overlap several fields. The abstracts of 
contributed papers have been appended to the chapters in accordance with 
the topic. Several of the contributed papers are very closely related to 
particular longer papers discussed in the main part of the chapter. Two 
lectures were presented at the Conference in addition to the papers 
appearing here. Professor Jurgen Moser of the Courant Institute of 
Mathematical Sciences presented a new formulation of Sacker's result on 
continuation under perturbations of smooth invariant manifolds. Professor 
Stephen Smale, University of California at Berkeley, spoke on 'Global 
stability questions 1 , concerning dif f eomorphisms of a smooth, compact 
manifold; he discussed necessary and sufficient conditions for structural 
as well as $\Omega$-stability . His lecture will appear in the Proceedings 
of the 1968 AMS Summer Institute in Global Analysis. 1 1 

Table of Contents: John A. Nohel, Preface (pp. xi-xvi). 

Chapter 1: Problems in partial differential equations: Invited lectures: 
D. G. Aron'son, Local behavior of solutions" of nonlinear parabolic * equations 
(pp. 3--8); Avner Friedman, Nonlinear eigenvalue problems (pp. 9 — 13); 
Fritz John, Plane waves of finite amplitude (pp. 14 — 19); P. H. Rabinowitz, 
Existence of periodic solutions for some nonlinear hyperbolic partial 
differential equations (pp. 20--24); E. B. Fabes and N. M. Riviere, $L\sp 
p$-estimates $ ( Kp\leq\inf ty ) $ near the boundary for solutions of the 
Dirichlet problem (pp. 2S--34); James Serrin, Existence theorems for some 
compressible boundary layer problems (pp. 35--42). 

Abstracts of contributed papers: J. W. Dettman, Related well-posed and 
ill-posed problems in partial differential equations (p. 45); B. A. Fusaro, 
The Cauchy-Kowalewski theorem and a singular initial value problem (pp. 
46--47); S. P. Hastings, An existence proof for a boundary value problem 
arising in boundary layer theory (pp. 47--48) ; S. Kaniel, On the motion of 
a viscous incompressible fluid (p. 48); D. Mangeron and M. N. Oguztoreli, 
Boundary value and optimal problems concerned with various types of 
N% polyvibrating 1 1 equations (pp. 48--49); Alan C. Newell, A note on the 
multiple time scale method in perturbation problems (p. 49) ; B. J. 
Matkowsky, Asymptotic solution of a nonlinear stability problem (pp. 



49--51) ; D. Sattinger, Quantum scattering and asymptotic series (p. 51); R. 
E. Showalter and T. W. Ting, Pseudoparabolic partial differential equations 
(p. 51); Gabe Buis, Martin Eisen and W. G. Vogt, Lyapunov stability of 
partial differential equations (p. 52). 

Chapter 2: Delay and integral equations: Invited lectures: C. 
Corduneanu, Admissibility with respect to an integral operator and 
applications (pp. 55 — 63); Stephen Grossberg, A global prediction (or 
learning) theory for some nonlinear functional-differential equations (pp. 
'64-^70); Jac"k K.' Ha'le/ Behavior near' aperiodic 6'rbit v o*f ' functional * 
differential equations (pp. 71 — 75); J. J. Levin, Boundedness and 
oscillation of some Volterra and delay equations (pp. 76 — 81); R. K. 
Miller, The topological dynamics of Volterra integral equations (pp. 
82--87) . 

Abstracts of contributed papers: R. B. Grafton, A periodicity theorem 
for functional differential equations with applications (p. 91); Melvin L. 
Heard, On asymptotic behavior and periodic solutions of a certain Volterra 
integral equation (p. 92); A. P. Stokes, On the stability of integral 
manifolds for functional differential equations (pp. 92--93) ; Elliot 
Winston, Uniqueness of the zero solution for state dependent delay 
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second order equations (pp. 159--165) ; Taro Yoshizawa, Some remarks on the 
existence and the stability of almost periodic solutions (pp. 166--172) . 

Abstracts of contributed papers: M. S. Berger, On periodic solutions of 
Hamiltonian systems of ordinary differential equations (p. 175); P. 
Brunovsky, On stabilization of periodic solutions of nonlinear systems (p. 
175); Craig Comstock, On the limit cycles of $\ddot y-\mu\sin\dot y+y=0$ 
(pp. 175--177); A. M. Fink, Almost periodic solutions with module 
containment (p. 177); D. A. Sanchez, A note on periodic solutions of 
Riccati-type* equations* (pp. 177 — TIB).'* ' 

Chapter 5: Further topics in ordinary differential equations: Invited 
lectures: H. A. Antosiewicz, Set-valued mappings and differential equations 
(pp. 181--183) ; W. A. Harris, Holomorphic solutions of nonlinear 
differential equations at singular points (pp. 184 — 187); Henry Hermes, 
Existence and properties of solutions of $\dot x\in R(t,x)\sp \ast$ (pp. 
188--193) ; Aaron Strauss, Perturbing asymptotically stable differential 
equations (pp. 194 — 198). 

Abstracts of contributed papers: Arrigo Cellina, Multivalued 
differential equations and ordinary differential equations (pp. 201--202); 
J. W. Bebernes and Russell Wilhelmsen, A technique for solving 
two-dimensional boundary value problems (pp. 202--203); G. B. Gustafson, 
Examples in the oscillation theory of higher order equations (p. 203) ; 
Frank Hoppensteadt , Asymptotic expansions for singularly perturbed systems 
(pp. 203--204); J. D. Schuur, Saddle point type behavior for nonautonomous 
second order almost linear systems of differential equations (pp. 
204--205) ; James A. Yorke, Some extensions of Lyapunov 1 s second method (pp. 
206—207) . 
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following title: Advances in differential and integral equations (Conf. 
Qualitative Theory of Nonlinear Differential and Integral Equations, Univ. 
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Justification of matched asymptotic expansion solutions for some 
singular perturbation problems . 

Partial differential equations (Proc. Sympos . Pure Math., Vol. XXIII, 
Univ. California, Berkeley, Calif., 1971) 
Hoppensteadt , Frank 
1973, 

Airier. Math. Soc, Providence, R.I.,; pp. 337--341., 
Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: SHORT (10 linesj 

The author* gives* a Tucid arid concise discussion of methods for 
constructing asymptotic expansions for solutions of initial value problems 
of the form $\varepsilon u 1 (t) =F(t, u, \varepsilon) $, $u(0)=u\sb 
0 (\varepsilon) $, where $0\leq t\leq T\leq\inf ty$ , $u(t)$ belongs to a 
Banach space $E, \varepsilon$ is a small positive parameter, and $F$ is a 
mapping into $E$ . References are given to problems for which these methods 
have been verified. 

\{For more complete bibliographic information about the collection in 
which this article appears, including the table of contents, see MR 
48\#6620.\} 

Reviewer: Greenlee, W. M. 

Review Type: Signed review 

Descriptors: *35B25 -Partial differential equations-Qualitative 
properties of solutions-Singular perturbati ons ; 34G05 — Ordinarv 
differential equations-Differential equations in abstract spaces (See also 
58D25 ) -Differential equations in Banach and other abstract spaces (See also 
47Bxx.) 



1/5/296 (Item 86 from file: 239) 

DIALOG (R) File 239 :*Math'sci " 4 . % ...... . 

(c) 2004 American Mathematical Society. All rts. reserv. 
01338536 MR 49##3277 

Asymptotic stability in singular perturbation problems. II. Problems 
having matched asymptotic expansion solutions . 
Hoppensteadt, Frank 

J. Differential Equations 
1974, 15, 510 — 521. 
Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (24 lines) 

As in Part I [same J. 4 (1968), 350--358; MR 37\#1730] , the author 
considers the system (1) $dx/dt=f ( t , x, y , \varepsilon) $ , 

$\varepsilon\,dy/dt=g(t,x, y, \varepsilon) $ together with the reduced system 
(2) $dx/dt=f (t,x,y,0)$, $0=g ( t , x, y, 0 ) $ and assumes that (2) has a solution 
$(x\sb 0(t),y\sb 0(t))$ lying on a smooth branch of $g ( t , x, y, 0) =0$ . He 
makes two stability assumptions: (i) The linear system $dz/dt= ( f \sb x-f\sb 



yg\sb y{}\sp {-l}g\sb x)(t,x\sb 0(t),y\sb 0(t),0)z$ is exponentially 
asymptotically stable; (ii) the real parts of the eigenvalues of $g\sb 
y(t,x\sb 0(t),y\sb 0(t),0)$ are uniformly bounded above by a negative 
number. He shows that if $\xi ( \varepsilon) $ and $\eta ( \varepsilon) $ are 
smooth functions of $\varepsilon$ at $\varepsilon=0$ , and $\xi(0)$ and 
$\eta(0)$ are in certain 4 'domains of attraction 1 ' (independent of 
$\varepsilon$ ) , then (1) has a unique solution satisfying 
$x (0) =\xi (\varepsilon) $, $y ( 0) =\eta (\varepsilon) $ ; and if $ (x, y ) , ( \tilde 
x,\tilde y) $ are two such solutions (belonging to $(\xi,\eta)$ and 
$ ( \tilde\xi, \tilde\eta) $, respectively), then 

$ (x ( t, \ va reps i Ion) , y ( t, \varepsilon) ) - (\ tilde x ( t, \ va reps i Ion) , \ tilde 
y ( t, \varepsilon) ) \rightarrow 0$ as $t\rightarrow\inf ty$ . Under additional 
assumptions such solutions tend to finite limits as $t\rightarrow\inf ty$, 
and the author shows how the asymptotic expansions of these limits in 
powers of $\varepsilon$ can be computed. 
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A problem in the theory of epidemics. II. 
Hoppensteadt , Frank 

Waltman, Paul 

Math. Biosci. 

1971, 12, 133 — 145. 
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The authors modify earlier work JMatl\. Biosci . , .9, .( 197X1). , 71--91* MR . 
44\#7083] to obtain what they feel is an improved model of the spread of 
infection in an epidemic. In this new version an individual who is 
recovered from the infection is allowed to become reinfected. (By allowing 
for an infinitely long period before this reinfection is possible, they 
include the previous model.) Their scheme results in a set of five 
nonlinear integral equations, and they prove the basic results on 
existence, uniqueness and continuous dependence for these equations. They 
also discuss some numerical computations done on the model. 

Reviewer: Hastings, S. P. 
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Hoppensteadt, Frank 
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Math. Biosci. 
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Language: English 
Document Type: Journal 
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The authors consider a model, proposed by K. Cooke [Differential 
equations and dynamical systems (Proc. Internat. Sympos., Mayaguez, P.R., 
1965), pp. 167 — 183, Academic Press, New York, 1967; MR 36\#5461] , for th 



spread of an infection. The model leads to a functional differential 
equation where the delay depends on time, together with a coupled integral 
equation. The authors show existence, uniqueness, and continuous dependence 
on initial data. Then they discuss approximate solutions and numerical 
methods . 

Reviewer: Hastings, S. P. 
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Properties of solutions of ordinary differential equations with small 
parameters . 

Hoppensteadt , Frank 

Comm. Pure Appl . Math. 
1971, 24, 807--840 
Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: MEDIUM (22 lines) 

Initial value problems of the form $dx/dt=f ( t, x, y, \varepsilon) $ , 
$\varepsilon\, dy/dt=g ( t, x, y, \varepsilon) $, $x (0) =\xi ( \varepsilon) $, 
$y (0) =\eta (\varepsilon) $, where $x$ and $y$ are vector functions, are 
solved asymptotically for small $\varepsilon$ under assumptions more 
general than in the previous literature. As in the work of A. B. Vasileva 
[e.g., in Uspehi Mat. Nauk 18 (1963) no. 3 (111), 15--86; MR 28\#1363] and 
of others, the solution is represented as a sum of an ""outer 11 and an 
""inner 11 solution. The main new features of this work appear to be the 
following: (1) Let $x\sb 0(t)$, $y\sb 0(t)$ be the solution of the 
" "reduced* 1 system obtained for $\varepsilon=0$ , of which the solution to 
be found for" $\vare , psrlon>0$ "is a 'singiflar perturba'tidn." ' The Eigenvalues of 
the Jacobian matrix $\partial g(t,x\sb 0(t),y\sb 0 ( t ) ) Apartial y$, which 
are crucial in the theory, are subjected to milder conditions than in 
previous work. Some of them may even cross the imaginary axis, as'$t$ 
changes. (2) The " "matching' ' condition which determines the initial values 
at $t=0$ of the outer solution is put into a significantly simpler and more 
natural form. 

Short descriptions are given of extensions of the theory to two-point 
boundary value problems and to systems with several parameters. 
Reviewer: Wasow, W. 
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On quasilinear parabolic equations with a small parameter. 
Hoppensteadt ; Frank" * * % 

Comm. Pure Appl. Math. 
1971, 24, 17--38 
Language: English 
Document Type: Journal 
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The main purpose of this paper is to derive an expansion for the solution 
of the quasilinear boundary value problem $$ \multline \varepsilon ( \partial 
uApartial t)-\surn\sb {i,j = l}\sp na\sb { i j } (x, t, u, \nabla 
u, \varepsilon) \ , \partial\sp 2u/\partial x\sb i\partial x\sb j=\\ 
f (x, t, u, \nabla u, \varepsilon) , \endmultline $$ $u (x, t , \varepsilon ) =0$ , 



$x\in\partial\Omega$ / $0\leq t\leq T$, $u (x, 0, \varepsilon) =u\sb 
0 (x, \varepsilon) $, $x\in\Omega$ , in the smooth domain $\Omega\subset E\sp 
n$ which is valid as $\varepsilon\downarrow 0$. While the methods used here 
are formally similar to those used in the author's previous paper [Arch. 
Rational Mech. Anal. 35 (1969), 284--298; MR 40\#1694], where semilinear 
problems were treated and the resulting expansion was valid only in the 
$L\sp P(\Omega)$ norms, the current methods allow use of uniform norms 
everywhere . 

The author mentions that more general problems can be treated by 
variants of the present method, e.g., higher order parabolic equations with 
more general boundary conditions, systems of parabolic type, and parabolic 
equations coupled with integro-dif f erential equations. 
Reviewer: Sigillito, V. G. 
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Cauchy problems involving a small parameter. 
Hoppensteadt , Frank 

Bull. Amer. Math. Soc. 
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For the initial value problem 
$\varepsilon (dv/ dt) -A ( t, \ va reps i Ion) v=f ( t, v, \ va reps i Ion) $, $0\leq t\leq T$, 
$v(0)=v\sp 0 (\varepsilon) $ , where $v$ is an element of a Banach space $E$ 
and $\varepsilon>0$ is a small parameter, a method is given for finding an 
expansion for the solution which is valid as $\varepsilon\rightarrow 0$. 
The- expansion has -the -form $\sum\sb (r^^sp \inftyTvVsT5 ' r (t)"+V\sB ' 
r ( t/\varepsilon) -U\sb r ( t/\varepsilon) ] \varepsilon\sp r$, where the $v\sb 
r$ give the solution for $t$ away from zero, the $V\sb r$ give it for $t$ 
near zero, and the $U\sb r$ are matching terms to make the representation 
valid uniformly in $0\leq t\leq T$ . Three cases are treated: (i) Abstract 
parabolic case where, for each 

$ ( t , \varepsilon) \in [ 0, T] \times [ 0, \varepsilon\sb 0]$, $-A ( t , \varepsilon ) $ is 
the infinitesimal generator of an analytic semigroup of operators in $E$ . 
(ii) Abstract hyperbolic case where, for each $ ( t, \varepsilon) $, 
$-A( t, \varepsilon) $ is the infinitesimal generator of a semigroup of class 
$C\sb 0$. (iii) Parabolic case where $A$ is a positive definite elliptic 
operator in $E=L\sp \infty$. 
Reviewer: Hersh, R. 
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Asymptotic series solutions of some nonlinear parabolic equations with a 
small parameter. 
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Si considera il problema (1) 



$\varepsilon\, dv/dt+A(t, \varepsilon) v=f ( t, v, \varepsilon) $, $v ( t\sb 
O)=\overset 0\to v$, $v\in E$, ove $E$ e uno spazio di Banach, 
$\varepsilon$ e un parametro positivo, $A( t, \varepsilon) $ sono operatori 
lineari, chiusi, anche illimitati, con dominio di definizione 
$D=D (A( t, \varepsilon) ) $ denso in $E$ e indipendente da $t$ e da 
$\varepsilon$ per $0\leq t\leq T$, $0\!eq\varepsilon\leq\varepsilon\sb 0$; 
$f ( t, v, \varepsilon) =\sum\sb { i, j , k=0 } \sp \inf ty\varepsilon\sp it\sp jf\sb 
{ijk)v\sp k$ e $f\sb {ijk}$ sono operatori $k$-lineari, limitati in $E$ . Si 
suppone il problema (1) di tipo parabolico (si ammette cioe che per ogni 
$ ( t, \varepsilon) $, $-A ( t, \varepsilon) $ sia generatore di un semigruppo 
analitico di operatori limitati in $E$) . Al fine di pervenire a uno 
sviluppo asintotico della soluzione del problema considerato, l ! autore 
studia dapprima la soluzione fondamentale relativa alia parte lineare di 
(1), dandone una valutazione per $\varepsilon\rightarrow 0$. 
Successivamente dimostra, sotto opportune ipotesi per gli operatori $A$ e 
le f unzioni* "$f $ * e *$'v$,* che il' problema"^ 1) " ammette*, "pe*r "ogni $\varepsilon$ 
suf f icientemente piccolo, una sola soluzione $v ( t , \varepsilon) $ continua e 
dif f erenziabile per $0\leq t\leq T$ e assegna, per tale soluzione, supposto 
$\overset 0\to v\in D$, uno sviluppo asintotico valido per 

$\ varepsilon\rightarrow 0$. II metodo usato puo essere esteso per trattare 
problemi con condizioni iniziali piu generali di quelle considerate in (1). 
Reviewer: Vaghi, C. 
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On systems of ordinary differential equations with several parameters 
multiplying the derivatives . 
Hoppensteadt , Frank 

J. Differential Equations 

' 1969, 5y 106-'- lie ■ • " ..... , 

Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
Abstract Length: LONG (29 lines) 

The initial value problem $x 1 =f ( t , x, y\sb 1, \cdots , y\sb 
n, \varepsilon) , x ( t\sb 0)=x\sb 0 , \varepsilon\sb jy\sb j{}'=g\sb j(t,x,y\sb 
1, \cdots, y\sb n, \varepsilon) , y\sb j (t\sb 0)=y\sb { j 0 ) \ , ( j = l, \cdots , n) $ is 
studied asymptotically, as the $n$ positive parameters $\varepsilon\sb j$ 
tend to zero in such a way that $\varepsilon\sb { j +1 } /\varepsilon\sb 
j\rightarrow 0$ for $ j=l , \cdots , n$ . Here, $x,y\sb 1 , \cdots , y\sb n$ and 
$f,g\sb 1, \cdots, g\sb n$ are real vectors of dimensions $k,k\sb 
1, \cdots, k\sb n$, respectively, while $\varepsilon$ is the vector 
$ (\varepsilon\sb 1 , \varepsilon\sb 2 , \cdots , \varepsilon\sb n)$. A set of 
sufficient conditions is given under which the solution of the initial 
value problem tends, uniformly in unbounded $t$-intervals , to the solution 
of the reduced problem obtained formally by setting $\varepsilon=0$ and 
canceling the initial conditions on the $y\sb j$. The conditions involve 
the uniform smoothness of the coefficients on the infinite $t$-interval and 
the existence and stability of the solutions of certain auxiliary 
differential' equations' that may appropriately be ' call e*d v boundary layer 1 ' 
equations. The special case $n=l$ was the subject of the author ! s paper in 
Trans. Amer. Math. Soc. 123 (1966), 521--535 [MR 33\#2900] . The method of 
proof is similar to that of the earlier paper and involves the construction 
of Ljapunov functions for the boundary layer equations, but there are 
non-trivial differences, caused by the occurrence of a hierarchy of $n$ 
boundary layer equations. The author 1 s counterexamples in this as well as 
the earlier paper indicate that his hypotheses cannot be substantially 
weakened if uniform convergence on unbounded $t$-intervals is desired. 
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Asymptotic series solutions for nonlinear ordinary differential 
equations with a small parameter. 
Hoppensteadt , Frank 

J. Math. Anal. Appl . 
1969, 25, 521--536 
Language: English 
Document Type: Journal 

Subfile: MR (Mathematical Reviews) AMS 
■ Abstract Length:- 'SHORT (5 lines)" • • • * - 

In the author's previous work on singular perturbations [Trans. Amer. 
Math. Soc. 123 (1966), 521 — 535; MR 33\#2900] , he studied the dependence of 
solutions on a parameter for the infinite interval. Here he obtains an 
asymptotic expansion in powers of the parameter. 
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Asymptotic stability in singular perturbation problems. 
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The paper deals with real systems (1) $x 1 =f ( t, x, y, \varepsilon) $ , 
$\varepsilon y 1 =g ( t , x, y, \varepsilon ) $ , $x(t\sb 0)=x\sb 0$, $y(t\sb 0)=y\sb 
0$, where $x$ and $y$ are, respectively, $n$- and $m$-dimensional vectors, 
and $\varepsilon$ is a small positive parameter. For the sake of simplicity 
it is assumed that $ f ( t , 0, 0 , 0 ) =g ( t , x, 0, 0 ) =0$ . Two degenerate systems 
associated with (1) are (2) $x ' =f ( t, x, 0, 0) $, $x(t\sb 0)=x\sb 0$, and (3) 
$dy/ds=g(\alpna,\beta,y,0)$, $y(0)=y\sb 0$ ($\alpha$ and $\beta$ are 
parameters) . Appropriate differentiability and boundedness conditions are 
imposed on $f$ and $g$ and it is assumed that the solutions $x=0$ and $y=0$ 
of (2) and (3) are uniform-asymptotically stable, the latter uniformly in 
$\alpha, \beta$ . Theorem 1: Under the further assumption that $f$ and $g$ 
are periodic functions of $t$ (with a common period) , for each sufficiently 
small positive $\varepsilon$ there exists a uniformly asymptotically stable 
tube of solutions of (1). Theorem 2: Under certain further assumptions, 
including the assumption that $f$ and $g$ are independent of $\varepsilon$ , 
the solution $x=y=0$ for (1) is asymptotically stable (Ljapunov) for each 
sufficiently small positive $\ varepsilon$ . 
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$\varepsilon\rightarrow 0\sp +$ , of tine solutions of the initial value 
problem $x 1 =f ( t , x, y, \varepsilon) $ , $x(t\sb 0)=x\sb 0$; $\varepsilon 
y'=g(t,x,y,\varepsilon)$, $y(t\sb 0)=y\sb 0$, for $t\sb 0\leq t<\infty$. 
Here, $f$ and $g$ are real vectors of dimensions $k$ and $j$, 
respectively. Conditions are given under which the solution converges 
uniformly, as $\varepsilon\rightarrow 0\sp +$, not only on compact sets, 
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series of examples that uniform convergence need not take place if any of 
his conditions is omitted. One of these examples deserves particular 
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convergence on infinite intervals of linear singular perturbation problems. 
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functions, constructed with the help of a lemma of Massera [Ann. of Math. 
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